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Die Biorthonormierung ist bei linearen Matrizeneigenwertproblemen im Falle von Mehrfachwurzeln bei 
nichtsymmetrischen Matrizen von Bedeutung. Zur praktischen Durchführung wird eine Vorgehensweise an- 
gegeben, die unter Verwendung der allgemeinen Gramschen Matrix mit dem modernisierten Gaußschen Algo- 
rithmus arbeitet. Lineare und nichtlineare Elementarteiler, sowie reelle und komplexe Bigenwerte werden ge- 
trennt, aber nach einheitlichen Gesichtspunkten bearbeitet. 


The bi-orthonormal transformation is of importance for linear matrix problems involving multiple latent 
roots of non-symmetric matrices. For the numerical evaluation, the method given in this paper, uses the genera- 
lized Gram determinant and the modified Gauss’ algorithm. Linear and non-linear elementary divisors as 
well as real and complex latent roots are treated separately, but from a unique point of view. 


La transformation bi-orthonormale est d’importance & des probl&mes de matrices lineaires au cas de racines 
latentes multiples a des matrices non symmetriques. Pour l’Evaluation numerique une methode est indiquee, 
qui combine l’application de la matrice de Gram gin£rale avec l’algorithme de Gauss modifie. Les diviseurs 
lementaires lineaires et non lineaires et les racines latentes reelles ei complewes sont traites separement, mais 
suivant des points de vue uniformes. 


B ımHeäHnsIX IPoÖJeMaX COÖCTBEHHBIX 3HAYeHÜ MATPHI CYINeCTBEHHEIM ABAIAETCH ÖHOPTO- 
HOPMHFOBAHHE® B CIyyae KpATHBIX KOPHeÜ Y HeCHMMEeTPHYecKHX MaTtpnun. ua nparrayec- 
KOoü O0Pa0OTKH AaeTcA CIIOCOd, OCHOBBIBAWINMHÜCH HA YCOBePIIeHCTBOBAHHOM TAYCCOBOM AATO- 
prdMe Ipn HCNOANB30BAHHM O00ÖmMeHHOÄ Marpausı I'pama. JIumeürnste u HesImHeiiHgte 3TeMeH- 
TAPHbIe NEJIHTEJAH, a TAKIKEe NelicTBMTe/IbHBIe H KOMILIEKCHBIE COÖCTBEHHBIE 3HAYEeHHA, OÖPa- 
ÖATEIBAWTCH OTAENBHO, ONHAKO C EAHHOÄ TOYEH 3peHHn. 


1. Einleitung 
Gegeben sei das Eigenwertproblem 


er la tina Ka 21% 2 bau 6) 7013 


einer n-reihigen quadratischen Matrix U mit beliebigen reellen Elementen «a,;. Ihm zugeordnet 


ist das transponierte Problem 
ns. et (02) 


der transponierten Matrix W’, das die gleichen Eigenwerte besitzt wie (01). Die Eigenvektoren 
xt; von (Ol) und ); von (02) sind, wie man leicht zeigen kann, im Falle durchweg verschiedener 
reeller Eigenwerte A; zueinander orthogonal, und bei geeigneter Normierung gilt 
Our ti ek, 03 
ET Ward re eu nee (WI) 
Das System der 4, stellt dann zu dem der g; ein reziprokes Grundsystem im Sinne der Vektor- 
rechnung dar [3] und umgekehrt. Nur im Falle symmetrischer Matrizen, die uns hier nicht 
interessieren sollen, fallen beide Grundsysteme zusammen, die Eigenvektoren 7, sind unter- 
einander orthogonal. 
Im Falle durchweg verschiedener reeller Eigenwerte läßt sich ein beliebiger Vektor tv des 
n-dimensionalen linearen Vektorraumes R, nach den Eigenvektoren r, entwickeln: 
wear’ Fu SS euuud a  I aa, Ber Va Zn (04), 
und die Entwicklungskoeffizienten a; bestimmen sich mit Hilfe des reziproken Grundsystems 
der h; nach 
Hy 
—— —- = ;w ee ie re ae. (05). 
yii ü 
Dies folgt unmittelbar aus (04) durch Vormultiplikation mit 4; unter Beachtung von (03). 
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Sucht man demgegenüber einen Vektor 3, der eine Reihe von Eigenvektoren von (01) 
nicht enthält, sagen wir X, I», - . ., t, (r<n), so müssen die Entwicklungskoeffizienten a}, @,, 
., a, zu Null werden, der Vektor 53 muß also zu dem zugehörigen reziproken Teilgrundsystem 


orthogonal sein: 
v5 = Il el a 


Diese Aufgabe taucht bei der praktischen Behandlung des Eigenwertproblems auf. Hat man 
nämlich bereits einen oder mehrere Eigenvektoren von (01) ermittelt, so läßt sich das Problem 
mit Hilfe eines derartigen Vektors 3 reduzieren derart, daß nur noch an 3 beteiligte Eigenwerte 
und Eigenvektoren auftreten. Dies ist bekannt vom gewöhnlichen Iterationsverfahren her. 
Es ist aber von gleicher Bedeutung für die direkten Verfahren zur Auflösung der Eigenwert- 
aufgabe, was in einer weiteren Arbeit gezeigt worden ist?). 

So einfach wie geschildert sind die Verhältnisse indessen nur bei einfachen Eigenwerten. 
Im folgenden beschäftigen wir uns demgegenüber mit dem Auftreten von Mehrfachwurzeln 4; 
von der Vielfachheit 9, > 1, wobei wir zwischen linearen und nichtlinearen Elementarteilern 
der charakteristischen Determinante zu unterscheiden haben. Während sich im ersten Falle 
gegenüber dem der einfachen Wurzeln nur wenig ändert, hat man im zweiten das dann unvoll- 
ständige System der Eigenvektoren durch sogenannte Hauptvektoren (vgl. [1] S. 306, 
[12] S. 211) zu einem vollständigen Vektorsystem vom Range n zu ergänzen, und zwar sowohl 
das System der x, als auch das der 4,. Diese aus Eigen- und Hauptvektoren bestehenden Vek- 
torsysteme sind dann so zu biorthonormieren, daß auch im allgemeinen Falle die den Be- 
ziehungen (05) und (06) entsprechenden angegeben werden können. Das Ziel der Arbeit ist es, 
für die. Biorthonormierung eine in jedem Falle anwendbare praktische Vorgehensweise zu ent- 
wickeln. Es kann dann insbesondere ein Vektor 3 bestimmt werden, der keinen zu einem mehr- 
fachen Eigenwert gehörigen Eigen- oder Hauptvektor mehr enthält. 

Unabhängig vom Eigenwertproblem wird in 2. zuächst die Biorthonormierung zweier 
Vektorensysteme besprochen. Nach kurzem Hinweis auf die Bestimmung der Eigenvektoren tr; 
und dh, in einem Arbeitsgang in 3. werden in 4. einige Sätze über Hauptvektoren aufgeführt, so- 
weit sie für die Biorthonormierung gebraucht werden. Diese wird in 5. für lineare und in 6. für 
nichtlineare Elementarteiler durchgeführt. In 7. werden dann die Beziehungen (05) und (06) 
auf den Fall mehrfacher Eigenwerte verallgemeinert. Während bis dahin alle Eigenwerte als 
reell vorausgesetzt werden, gehen wir abschließend in 8. kurz auf den Fall komplexer Eigen- 
werte ein. 

Die Theorie der Biorthonormierung geht auf eine Arbeit von Enskog [2] zurück (vgl. 
auch [4] und [5] S. 278—286) und ist im wesentlichen vonSchmeidler ([6] S. 75—78) auf 
finite Matrizen übertragen worden. Hier soll ein praktisches Vorgehen zur zahlenmäßigen Be- 
stimmung des Biorthonormalsystems angegeben werden. Wir benutzen dabei die verallge- 
meinerte Gramsche Matrix in Verbindung mit dem Gau Bßschen Algorithmus, wobei von 
Beziehungen für Untermatrizen Gebrauch gemacht wird. 


2. Biorthonormierung zweier Vektorsysteme 


Gegeben seien im linearen Vektorraum R, zwei Vektorsysteme B und ® mit je o linear 
unabhängigen Vektoren (o<n): 


VB = (DD, RE, in, Wu a en ee: 
Die zugehörige verallgemeinerte Gr am sche Matrix ist 
TV BE ERW ET ZI IByDE Da er 


Wir haben dann zwei Fälle zu unterscheiden, je nachdem der Rang von X gleich oder kleiner 
ist als o. 
a))Gramsche Matrix ohne Rangabfall 
B o-reihige Matrix N habe den vollen Rang o. Dann gilt 


Satz 1: | Genau dann, wenn die o-reihige Gr am sche Matrix (2) der Vektorsysteme (1) den 
Rang o besitzt, können die beiden Systeme ®, W biorthonormiert werden zu ®*, 
W*, so daß 


Kr WB u... Da 


Der folgende Beweis gibt zugleich die praktische Durchführung der Biortkonormierung. Wegen 
det % + 0Oläßtsich N aufspalten in eine untere und eine obere Dreiecksmatrix Q und © nach 


RE) ara Mir ee 


Re 2) JuUbe direkte Verfahren bei Matrizeneigenwertproblemen. Wiss. Z. Techn. Hochsch. Dresden, 2 (1952/53) 
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(vgl. etwa [9] S. 9, [12] S. 251). Indem man diesen Eliminationsprozeß auf die neben und unter 
X stehenden Vektorsysteme W bzw, ® ausdehnt gemäß 


De N a DE (5), 
erhält man die gesuchten biorthonormierten Matrizen W* und B*; denn für sie gilt 
Wr B* — DW BG = DIN GT - 10661 = 6, 
Gl. (3) ist erfüllt. 


Die praktische Durchführung ist in £ 
Schema la, b gezeigt. Nach Berechnung von Br NEE I BZ 
Win Ia spaltet man in Ib nach (4) auf und v (0,n)- Ba 


berechnet im gleichen Arbeitsgang W* 

und %*, indem man W als rechte Seite 

eines Gleichungssystems der Matrix N be- 

handelt und dem Eliminationsprozeß unter- vo ||; 
wirft ; entsprechend ® als rechte Seite des 
transponierten Systems mit der Matrix W. 
Das Gleichungssystem wird dabei nur bis 
auf die gestaffelte Form umgerechnet mit 
den Dreiecksmatrizen Q, und ©, eine Aufrechnung der Unbekannten wie sonst beim Glei- 
chungssystem findet nicht statt. Die umgerechneten rechten Seiten sind dann die gesuchten 
biorthonormierten Matrizen W* und B*. Die Kehrmatrizen Q-! und &-! werden nicht explizit 
benötigt. Die Orthogonaleigenschaft (3) wird zweckmäßig nachkontrolliert, wie in Ib durch die 
Einheitsmatrix angedeutet. 

Die Systeme W*, 9* hängen noch von der besonderen Wahl der Dreiecksmatrizen Q und 
& ab, d.h. von deren Diagonalelementen, die für eine der beiden Matrizen noch frei wählbar 
sind. Indem wir hier g; = 1 wählen, liegt das System W*, 3* fest. Man kann jedoch die 
Biorthonormierung auch so vornehmen, daß eines der Ausgangssysteme, sagen wir ®, unver- 
ändert bleibt, worauf wir jedoch nicht näher eingehen wollen. 

Im Falle o = nrist unter der Voraussetzung linear unabhängiger vb, und iw; der Rang der 
Gramschen Matrix immer gleich n, diese Systeme können also immer biorthonormiert wer- 
den. Man spricht von einem vollständigen Biorthogonalsystem ([6] S. 74). 

Mit W —= ® erhält man die gewöhnliche Orthogonalisierung mit der symmetrischen posi- 
tiv definiten Matrix W%. Die Orthogonalisierung ist bei linear unabhängigen b; immer durch- 
führbar. Ihre Durchführung ist in der Biorthonormierung enthalten ([7]). 


(n,e)- 
Matrix Öchema Ia 


öchema Ib 


b)Gramsche Matrix mit Rangabfall 


Die Gramsche Matrix N besitze den Rang r< o, und die nicht verschwindende Deter- 
minante der Ordnung 7 möge bereits in der linken oberen Ecke von Xstehen. Wir teilen X und 


die Vektorsysteme B und W wie folgt auf: 


(a. IE 
No Nr’ u) To Rn ze fe, [210] 


Schema IIa. Nun läßt sich nur N. gemäß 

(4) aufspalten in zwei Dreiecksmatrizen, 
Road: (M 

Schema IIb. Ausdehnen des Eliminations- 


prozesses auf die Nachbarmatrizen NW, 
und Ro 2 ergibt die verwandelten 


Matrizen j 
RIND = Deere. rl), 


wobei wieder jede Spalte von N,, und Wo als rechte Seite eines Gleichungssystems der Koeffi- 
zientenmatrix N, aufgefaßt und dem Eliminationsprozeß unterworfen wird und entsprechend 
jede Zeile von N, und ®, als rechte Seite des transponierten Systems. Die Kehrmatrizen ©, ', 
85" werden wieder nicht explizit benötigt. 

Die so umgewandelten Teilsysteme 5 und WS sind biorthonormiert, es ist 


ee... (9) 


Schema Ta Schema ITb 


mit der r-reihigen Einheitsmatrix 6, 
Es fehlen noch die Teile ®f und Wf. Diese werden in IIb gebildet nach 


Breue a und H-—BNıa +... (0) 
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was eine Anwendung des Ga u ßschen Algorithmus auf Untermatrizen darstellt (vgl. [8] S. 3). 
Dann gilt, wie leicht nachprüfbar, 


re —-0, WM ea N SD 
(2? muß eine Nullmatrix sein, da andernfalls der Rang von X größer als 7 wäre. Die 
Vektoren von ®# sind orthogonal zu allen Vektoren aus W*, die von W] orthogonal zu allen 
aus B*. 
Insgesamt gilt damit 
Satz 2:| Hat die Gramsche Matrix (2) der Vektorsysteme (1) einen Rang 7<o, so führt 
eine Biorthonormierung auf 


Tr — Wr I, * -| 


al 4 el 


0.0 
mit der r-reihigen Einheitsmatrix &,. 
Auch hier kann die Normierung gegebenenfalls so vorgenommen werden, daß eines der 
Vektorsysteme, z. B. ®, unverändert bleibt. 
Die in den Schemata I aufgeführten Vektoren 3 und 3* werden in Abschnitt 5 erklärt. 


3. Bestimmung der Eigenvektoren x; und Y;, 
Für einen Eigenwert A, der Vielfachheit p; habe die charakteristische Matrix ® = YW—4; & 
einen Rangabfall 
G=n—r; mil. . ISIS 


Im Falle e; = p; besitzt die charakteristische Determinante nur lineare Elementarteiler zu 
) = 4, im Falle e;< p; auch nichtlineare. 

Es existieren dann bekanntlich e; linear unabhängige Eigenvektoren r;; zu (O1)(l=1,2, ..., e;) 
und ebensoviele Vektoren 4;; zu (02), die, worauf wir hier hinweisen möchten, beide in einem 
Arbeitsgang nach dem G au ßschen Algorithmus bestimmt werden können. Zu ihrer Berech- 
nung denkt man sich B = U —4; & aufgeteilt in vier Untermatrizen B,,, wobei die r;-reihige 
quadratische Teilmatrix 9}, gerade die nichtverschwindende Rangdeterminante besitzen 
möge, was nötigenfalls durch Zeilen- und Spaltenvertauschung erreicht werden kann. 

Das Gleichungssystem (O1) lautet dann: 
Ar 6)r (5 all: =) A NA re BE HEN EE 
Da Bay’ \Eın 
es spaltet sich auf in zwei Systeme, von denen nur das erste gebraucht wird; das zweite ist 
überzählig. 
Die Rangmatrix B;, läßt sich wieder aufspalten in untere und obere Dreiecksmatrix nach 


Bas ET ee In NIT Pe RE u 11, 


und bei diesem Eliminationsprozeß geht das erste System aus (Ola) über in das gestaffelte 
System 


GH SD Bin - ET REN Oh 


Indem man hier in bekannter Weise für tr beliebige, und zwar e; linear unabhängige Spalten 
einsetzt (z. B. die Spalten der e,-reihigen Einheitsmatrix), erhält man e; Teilvektoren r, und 


damit insgesamt e; linear unabhängige Lösungen %; 1, Li» - - - > Lie; - 
Zur Bestimmung der Lösungen d);; des transponierten Systems geht man aus von 
Bu ale) 
A ee Re 02a 
(ee Daprs Yıı ( 
und erhält aus der ersten Gleichung unter Berücksichtigung von (14) das gestaffelte System 
De ee BR eo 


dessen Matrix Q, aber schon vorliegt. Die Auflösung des transponierten Systems erfordert 
daher nur noch die zusätzliche Berechnung der e; rechten Seiten von (02b) aus beliebigen 
linear unabhängigen 4}; und das Aufrechnen der Unbekannten d,. Der Eliminationsprozeß 
(die Aufspaltung) der Gleichungsmatrix ®,, braucht für die beiden Systeme (Ola) und (02a) 
nur einmal durchgeführt zu werden. 


4. Eigenschaften von Hauptvektoren 
Um im Falle nichtlinearer Elementarteiler zu einer vollständigen Lösung mit 7; linear 
unabhängigen Vektoren zu gelangen, führt man neben den Eigenvektoren 


Gı=G Yamdyı mit, a2 mr 
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deren Ermittlung in 3. besprochen wurde, die sogenannten Hauptvektoren x7, und 1°, der 
Stufe o ein, definiert durch: 
AAN. lie), 
E en für ORDER: ti = Yir a) 

(W—A; Ey =dYir RER LOD). 
Eigenvektoren sind Hauptvektoren 1. Stufe (wenn möglich, werden diese auch im folgenden 
ohne oberen Index geschrieben.) 

Gehen wir von einem der e, Eigenvektoren r;, in (15) aus und berechnen nach (16a) 


Hauptvektoren aufsteigender Stufe o —=2,3,..., wobei jedesmal ein inhomogenes singuläres 
Gleichungssystem zu lösen ist, so sprechen wir bei 
ET al EI De EN EEE NN rt l17) 


voneiner Hauptvektorkette. In jeder dieser e;, Ketten kann man höchstens bis zur 
Endstufe q;, aufsteigen. Wir wollen annehmen, daß in (15) die Vektoren bereits so geordnet 
sind, daß man für ein 2 + 1 mindestens bis zur gleichen Stufe wie für l aufsteigen kann. Geht 
man also von %;, aus, so gelangt man zur niedrigsten, von r;,. aus zur höchsten Endstufe. 
Für die q;, legen wir demgemäß fest: £ 


rn A a Ey 
Weiter gilt für diese Endstufen: 
ei 
PH ee DE Sa EEE aa ee Fat en ce een ne (19). 
i=1 


Die Endstufe g;; der Kette Z kann jedoch nur erhalten werden, wenn sieh in dieser Kette keine 
Bestandteile von Ketten mit niederer Endstufe befinden. Im anderen Falle kann die tat- 
sächliche Endstufe erst nach einer noch zu besprechenden ‚‚Reinigung“ erhalten werden. 

Bei dem transponierten System (02) können wir in analoger Weise Vektorketten nach 
(16b) aufstellen, ausgehend von einem Vektor h;7, aus (15). 


Er (SO Pr EN en er l20): 
Für derartige Ketten (17) und (20) gilt nun die Beziehung: 
N nen. 


Diese erhält man, indem man 4,7, nach (16b) durch 9,7,+1 (U —A&) und dann (U —A6) 7, 
nach (16a) durch 177! ersetzt und das Vorgehen w-mal wiederholt. Die Werte der skalaren 
Produkte (21) bleiben also unverändert, wenn die Summe der oberen Indizes konstant bleibt. 

Beginnen wir damit, die erste Kette (17) mit 1 = 1 zu bestimmen, dann kann man sicher 
— g;, ist nach (18) die niedrigste Endstufe — bis zu g;, aufsteigen und in dem transponierten 
System sicher ein 4;j- angeben, das nach geeigneter Normierung der Beziehung genügt: 

DD BON RDR ar u. 220 (22), 
woraus nach (21) die gleiche Aussage für alle oberen Indizes besteht, deren Summe 1+ 4;, ist, 
also für 2+(,ı —1), 3+ (41 —2), usw. Da g;, auch für das transponierte System die 
niedrigste Endstufe darstellt, können wir uns die h,!,so geordnet denken, daß wir in (22)7 —1 
setzen können. Die Kette 9,7 steigt dann ebenfalls bis zu g;, auf. 

Vektorketten x?, und 9?,, die bis zur gleichen Endstufe g;, aufsteigen und durch eine 
Beziehung (22) miteinander verknüpft sind, bezeichnen wirals zugeordnete Vektor- 
kettenoderKettenpaare. Wir haben gezeigt, daß man das erste Paar mit 1 = l ohne 
weiteres ermitteln kann. Faßt man die Vektoren r?, und 47, eines solchen Kettenpaares zu je 
einer Matrix zusammen, wobei die r-Vektoren nach steigender, die y-Vektoren nach fallender 
Stufe geordnet sind und bildet mit ihnen die Gramsche Matrix 


eg rennen), 


so besitzt diese nach (21) und (22) den Rang g;,. Denn in W;, stehen unter der Hauptdiagonale 
nach (21) für o+r<gq;, nur Nullen, so daß nach (22) eine obere Dreiecksmatrix mit den 
Hauptdiagonalgliedern 1 vorliegt. Die endgültige Normierung zur Einheitsmatrix nach 2.a ist 
also besonders einfach. 
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Zur Ermittlung weiterer Ketten werden zunächst die Vektoren r;, und d;, für l=2, 3, ...,e, 
von Bestandteilen der ersten Kette ‚‚gereinigt“. Dehnt man diese Reinigung auch auf Vek- 
toren höherer Stufen aus, so kann man das Kettenpaar mit der Endstufe g;, bestimmen und 
normieren. Die Gramsche Matrix N;, analog (23) weist den Rang g;, auf. Fährt man so in 
der Bestimmung der Ketten fort, so erhält man e; Ketten, wobei man in der I-ten Kette, falls 
diese von allen Bestandteilen der früheren Ketten gereinigt ist, bis zur Endstufe g;, aufsteigen 
kann. 

Die Reinigung von bereits ermittelten Ketten niederer Endstufe, die in 6. ausführlich 
besprochen wird, geschieht dadurch, daß man durch geeignete Addition bereits normierter 
r-Vektoren die restlichen Ketten so abändert, daß diese orthogonal zu den Eigenvektoren 
des transponierten Systems werden. Damit kann man dann Hauptvektoren höherer Stufe 
ermitteln. Das Gleiche gilt für die )-Kette. 

Die Gramsche Matrix 


NEN Rn a Er 
gebildet aus den beiden zum Eigenwert A; gehörenden Hauptvektorsystemen 
a 2 4; i el 
Wald... [Hill ao) 


besitzt nach dem Gesagten den Rang 3 4; = p;. Dieser ergibt sich auch aus folgendem. 
i=1 


Hauptvektoren zu verschiedenen Eigenwerten A; #4; besitzen ähnlich wie Eigenvek- 
toren gewisse Orthogonalitätseigenschaften. Schreibt man nämlich (16b) füro =r,l1—=/ und 
) = 4, in transponierter Form auf, multipliziert von rechts mit za und subtrahiert die so ab- 


geänderte Gleichung von der mit ur vormultiplizierten Gl. (16a), so erhält man: 


Satz 3: Die Hauptvektoren 27, einer beliebigen Stufe o zum Eigenwert A; sind orthogonal 
zu allen Hauptvektoren Dr beliebiger Stufe 7 des transponierten Systems zum 
Eigenwert A, mit A, # 4: 

EA ig 

ah Barry | 

Il, 2,020, 

[% ms — 1,2 Pr) er 

Da die Gesamtheit aller zu (01) und zu (02) gehörenden Hauptvektoren eine Gramsche 

Matrix vom Range n besitzt, muß nach (26) die mit Y, und &, aus (25) gebildete Gr am sche 

Matrix den Rang 9; besitzen. Hauptvektorsysteme 9% und X bezeichnen wir dann als 

biorthonormiert, wenn die Gramsche Matrix zur Einheitsmatrix wird. 

Witte Bios ab Seisss hd Min Talale-0827) 


Die Durchführung der Normierung wird in 6. ausführlich besprochen. 


y.20, = 0 für #4 . (26). 


5. Biorthonormierung bei linearen Elementarteilern 


Hier erhält man nach der in 3. angegebenen Vorgehensweise zum p-fachen Eigenwert ?) 
4 gerade p linear unabhängige Eigenvektoren tr, zu (01) und ebensoviele Eigenvektoren 4; zu 
(02). Die Gramsche Matrix RN —= W X%, gebildet aus den beiden Vektorsystemen von je p 
Eigenvektoren, besitzt daher den vollen Rang 9, und die Normierung kann nach 2.a durch- 
geführt werden. Die Rechnung verläuft nach Schema Ia, b, indem man ® durch X und ® 
durch 9 ersetzt. X* = B* und Y* = W* sind die gesuchten biorthonormierten Systeme der 
Eigenvektoren. 

Das gleiche Schema läßt sich nun auch zur Bestimmung eines Vektors 3* heranziehen, 
der keinen zu A gehörigen Eigenvektor x; mehr enthält, der also zum Vektorsystem 9) ortho- 
gonal ist. Man geht aus von einem beliebigen Vektor 4, der nur linear unabhängig sein muß 
zum Vektorsystem &%. In Ia bildet man db = W’ 3 und berechnet in Ib 


"0-93 - 9% 
wo die letzte Beziehung als Kontrolle benutzt werden kann. Dann ist 
a Be A 


4) Wir lassen hier und im folgenden alle den Eigenwert A = A; kennzeichnenden Indizes i fort, da Ver- 
wechslungen nicht zu befürchten sind. 
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der gesuchte Vektor, Schema Ib. Denn hiermit wird 
Y z* > Y 3 — Yx*p* —d — YAOHTO-1H —dH— NN IH —(, 
wasin Ib durch die Kontrollspalte e geprüft wird. 


Wir fassen zusammen: Ein zum Vektorsystem & linear unabhängiger Vektor 3 wird von 
allen Komponenten aus & gereinigt, indem man ihn dem Gaußschen Algorithmus 
für Untermatrizen in der Form (28) unterzieht. 

Das gleiche Vorgehen läßt sich natürlich auch anwenden, wenn & die Gesamtheit der 
Eigenvektoren mehrerer verschiedener Eigenwerte A; mit durchweg linearen Elementarteilern 
darstellt, so daß 3* dann keinen dieser Eigenvektoren mehr enthält. 


6. Biorthonormierung bei nichtlinearen Elementarteilern 
In diesem Falle können nach 3. nur noch e= et! Eigenvektoren x, und 4, berechnet 


5 . . . „1 al . 5 = 
werden, die wir uns in den Matrizen &%ıo und 90 zusammengefaßt denken. Die mit ihnen 
gebildete Gramsche Matrix 


Re N 
weist einen Rangabfall von 
FREIEN. er 91% 2% 2& (30) 


auf°). Eine Biorthonormierung kann also nur noch im Sinne von 2.b vorgenommen werden, 
wobei die Gramsche Matrix die Form 
0,0 


annimmt. Gewisse Eigenvektoren bleiben orthogonal zu allen Vektoren des transponierten 
Systems, so daß zu ihnen nach (15a) bzw. (15b) Hauptvektoren zweiter Stufe ermittelt werden 
können. 

Zur Durchführung der Orthonormierung erweist es sich nun als zweckmäßig, die zu 
gleicher Endstufe gehörigen Vektorketten zu einer Matrixkette 

ET ne IE NIIT TERN FAN) 

von je o, Vektoren zusammenzufassen. Dabei gibt der obere Index die Stufe o der betreffenden 
Hauptvektoren an, während der untere Index s jetzt die Endstufe bezeichnet. Eine Kette 
heißt ‚‚leer‘‘, wenn keine Vektoren der Endstufe s existieren, wenn also 0, = ist. Die maxi- 
mal auftretende Stufe überhaupt, also die höchste Endstufe sei 8 und es ist 8 — gu (= 4; })- 
Es läuft dann in jeder Kette o von 1bis s und s läuft von 1 bis 8, und es gibt S Ketten, von 
denen freilich einige leer sein können. Dann gelten die Beziehungen 


S$ S$ 
eo [OR — el und a BEN); =» ae ne (32). 
s=l s=1 


Mit Hilfe der Matrixketten können alle Vektorketten gleicher Endstufe formelmäßig und 
auch numerisch gemeinsam behandelt werden, was das Vorgehen vereinfacht. 

Ebenso wie bei den Vektorketten führen wir auch für das transponierte System Matrix- 
ketten ein und sprechen wieder von zusammengehörigen Ketten oder Kettenpaaren, wenn 
beide gleiche Endstufe aufweisen. Die Gra m sche Matrix 


DR Sen ner ei einlansiee na » (33) 


hat dann den Rang o,. 
Für zwei beliebige Matrixketten 9} und &; gilt in Erweiterung von (21) die Beziehung 


a (la). 

Die Matrizen X und 9? müssen dabei der Kettenrelation genügen. 
Breeze)... 2.2. (16aa) 
Bene 2..08) 08.» (166bb). 


Die Untermatrizen (21a) der aus zwei vollständigen Ketten gebildeten Gr am schen Matrix 
sind für konstante obere Indexsumme o-+ r — konst untereinander gleich. Diese Gram- 


5) Die Größen ee stellen die sogenannten W eyrschen Charakteristiken des Eigenwertes A dar ([12] 
S. 212) und geben die Anzahl der Hauptvektoren der Stufe o an; esist 


Ss 
deo=p». 
o=1 


326 Unger, Zur Praxis der Biorthonormierung Bd. 88.Nr. ee 


sche Matrix ist quadratisch für s =t, wenn es sich also um ein Kettenpaar handelt, rechteckig 
mit mehr Spalten als Zeilen für s<t und mit mehr Zeilen als Spalten für s>t. Ordnet man 
bei der Bildung der Gramschen Matrix die &- und W-Glieder der Ketten analog wie in 
(23) an, 


(w 
N 94 “ 
| IC}, eis ED et a 


\yı ) 


so sind die auf parallel laufenden Schrägzeilen stehenden Untermatrizen untereinander gleich. 
Im Falle s=t stehen in der Hauptdiagonale s gleiche quadratische Untermatrizen vom 
Range o;. 

Die Durchführung der Normierung zeigen wir am einfachsten an einem speziellen Bei- 
spiel, das genügend allgemein gehalten ist, um den allgemeinen Sachverhalt hervortreten zu 
lassen. Wirnehmen an, daß zu dem p-fachen Eigenwert A drei nicht leere Matrixketten der End- 
stufen 1, 3 und 5 gehören, während die Ketten 2 und 4 leer seien. Die (von Null verschiedenen) 
Anzahlen o,, 03, 0, der in den Ketten 1, 3, 5 enthaltenen Eigenvektoren sind ohne Belang; es ist 
nur 0, +30,+4+50; = pund es ist 0, =0, = 0. Die Höchststufe ist 8 = 5. 

Es werden folgende Bezeichnungen gewählt: 

&%%o enthält die Hauptvektoren der Stufe o der s-ten und aller etwa noch folgenden Ketten 
sl, s2,2.., 
&%ı enthält die Hauptvektoren der Stufe o der s-ten Kette allein. 
Es gilt also 
Ks = = (%ı ’ Reit, 0)- 
% enthält schließlich die fertig normierten Hauptvektoren der Stufe o der s-ten Kette. 
Die gleichen Bezeichnungen gelten für die Y-Matrizen. Noch unbehandelte Ausgangsmatrizen 
(Rohmatrizen) werden mit einem Querstrich versehen. 


Zur Vereinfachung sei weiterhin angenommen, daß alle Vektoren bereits so geordnet 
seien, daß in den Gram schen Matrizen die Rangdeterminante jeweils in der linken oberen 
Ecke steht. 


Zur Durchführung der Biorthonormierung sind die beiden Schemata IIIa und IIIb vor- 
gesehen. Ausgehend von den Eigenvektoren werden in IIIa zunächst die Endstufenmatrizen 


%ı der &-Ketten bestimmt, so daß damit die Endstufen s selbst und die Anzahl der Ketten 


bekannt sind. Dabei erfährt die hierbei zunächst unverändert bleibende Matrix Yo eine den 
Ketten entsprechende Unterteilung. Nachträglich werden diese Eigenvektoren in IlIa noch 


so umgeformt zu Y,ı, daß letztere mit den X, bereits fertig orthonormiert sind. 
3 Die Matrizen %% und 9, bilden die Ausgangswerte für die Rechnungen in IIIb. Die zu 


X; gehörenden Matrizen niederer End-Stufe ze BL FH Ku werden durch Herunterrechnen 
nach (16aa) ermittelt, während die Y-Ketten nach (16 bb) von unten nach oben bestimmt 


werden müssen, ausgehend von 9, . Die endgültige Biorthonormierung erfolgt ebenfalls in IIIb. 
Im einzelnen verläuft die Rechnung folgendermaßen. 

Schema Illa. Ausgehend von &o und Yo mit je e! Eigenvektoren wirdin IlIladie Gram- 

sche Matrix (29) mit dem Rang g,< e! gebildet und die Rangmatrix %, die als bereits links 

oben stehend angenommen wird, aufgespalten in 


=D Oi Ay 


Die übrigen Teile von (29) werden nach dem Gau Asche Algorithmus in Untermatrizen 
(vgl. [8] S. 3) umgerechnet in: 


ZA T mi mil 
=; (1)- or a u Pe 
1 1 2 1 


2 und X sind Nullmatrizen, da im anderen Falle der Rang von (29) größer als g, wäre. 
Analog zu (35) werden die Eigenvektoren in &l, umgerechnet. 


=, ae U Ve 
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lc DI re e as e BT 2 i 
20 ist damit von der ersten Kette gereinigt und orthogonal zu Yo, so daß nach (16aa) eine 
Kettenmatrix zweiter Stufe bestimmt werden kann. 


7A “ (m, N,, R;, n,, 2, R,,n;) 


£4 


—r\ 


Fa 


Schema Ta 


Diese Kettenmatrix zweiter Stufe &3 wird in IIIa oben eingetragen unddie Gra m sche 
Matrix Yo Er gebildet. Die so entstehenden Matrizen 2%, ME und N2 werden analog g, mM; 
und Ns nach (35) behandelt, womit man &%, M2 und N? erhält. Die Matrix 2% bewirkt dabei 
die Reinigung von der ersten Kette. Auf Grund der gemachten Voraussetzung tritt keine 
Kette mit der Endstufe 2 auf, die zweite Kette ist leer. MZ und N werden damit zu Null. 
Berechnet man schließlich 


RE Ta RN U (37), 
so ist diese Matrix wieder orthogonal zu Yo , so daß X, — X}, gesetzt werden und eine Ketten- 
matrix 3. Stufe, nämlich x, bestimmt werden kann. 

Diese Matrix wird aufgeteilt in den Bestandteil a der Kette mit Endstufe 3 und den 
weiterführenden Bestandteil Ka 
Nach Bildung der Gra mschen Untermatrizen 
le Is und RE een. 0 2638) 
werden zunächst & und &/ nach (35) analog % umgerechnet zu %3 und £3. Weiter wird die 
Rangmatrix W& mit dem Rang o, gebildet und aufgespalten in zwei Dreiecksmatrizen. 


nr A EU ra De ı'\E 


Unterwirft man nun in entsprechender Weise die Matrizen N, Mi, Ni dem Gaußschen 
Algorithmus in Untermatrizen, so erhält man R3, M? und gemäß den Voraussetzungen Wi =. 


Die Umrechnung von 2 und x geschieht nach: 
Hi HE, Hu HMÜ-EM. (0). 

Die Ermittlung der dritten Kette ist damit zunächst abgeschlossen. Von %%, aus kann man, 
da diese Matrix wieder orthogonal zu Yıo ist, zu einer Kettenmatrix 4. Stufe aufsteigen, diese 
wieder von Bestandteilen früherer Ketten reinigen und schließlich zur 5. Stufe aufsteigen und 
die Endstufenmatrix der fünften Kette %, — die vierte Kette war wieder leer — ermitteln, 
womit in IIIa die Berechnung der &-Ketten beendet ist. Mr 

Zwischen den Endstufenmatrizen X}, X, und &53, und der Matrix Yo, durch die Ketten- 
bestimmung in SIR Ya, Yıı aufgeteilt, wurde zunächst nur eine „Halbnormierung‘“ er- 
zielt. Man erhält nämlich: 
Yıı a et 
dl, u | 0 |..:...... (a). 
Dsi N N 
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Um volle Normierung herbeizuführen, wird jetzt auch noch die re un si nach 
Midi, MM NY, MI MM NDN (42). 


Die (41) Dee Gramsche Matrix wird damit zur Einheitsmatrix. 


S 


"5 


SE 
“nal oln& 


are) 


I = =N,, 96, mitl<s 


2° hema IIIb. Aus Illa werden zunächst die Matrizen X],, %, und %, oben, die Matrizen 
11, Val und Y4 links eingetragen. Von den Endstufenmatrizen X, % und %, ausgehend 
werden nun die Matrizen X ( =s —1, s—2,...,1) durch ‚„Herunterrechnen‘ mittels 
der Kettenrelation (16aa) bestimmt (in IIIa sind die Matrizen nicht vorhanden!) und im 
Schemakopf eingetragen, der damit vollständig ist. Weiter liegen in der Gramschen Ge- 
samtmatrix alle Nullmatrizen und diejenigen mit den o,-reihigen Einheitsmatrizen €, bereits 
fest. 
Das erste Kettenpaar liegt mit X, = X! und Yl, = Qi, die also nach unten bzw. nach 
rechts übertragen werden können, bereits fertig normiert vor. 


Zur Ermittlung des dritten Kettenpaares wird zuerst Vi 3 nach rechts übertragen. 
Aus 9 wird nach (16bb) 97 bestimmt und mit der Gram len Teilmatrix 
. REIN RR 
von Bestandteilen der ersten Kette gereinigt zu: | 
DIEB RE ARE 


9), wird wieder rechts eingetragen und stellt eine Kettenmatrix zu W}dar. Es kann wiederum 
nach (16bb) eine Matrix 3. Stufe 95} ermittelt werden, die nach Reinigung mittels T#, in 9? 
übergeführt wird analog zu (43) und (44). Die Y-Kette liegt damit fertig vor. Um die end- 
gültige -X-Kette zu finden, sind die Matrizen 


Bea Mn, SV 900: (5) 
zu bilden, mit denen die Umrechnung vorgenommen wird zu: 
Mh, Bm —HC, MM — Re ie Zee 


Man erkennt, daß es sich um Kettenmatrizen handelt. Mit ‚l2B) ist die Biorthonormierung des 
dritten Kettenpaares beendet. 

Die Bestimmung des fünften Kettenpaares beginnt wieder mit der Berechnung der Y-Ket- 
ten. Ausgehend von Y,4 = 9,;', was sofort nach rechts übertragen werden kann, wird 9,? nach 
(16bb) bestimmt. Die Reinigungsmatrizen 


—h = Dit und TE N TREE 
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Dr eine Entfernung von Komponenten aus der ersten und dritten Kette, indem man 
ildet: 
ne... een: (AB), 
9; wird wieder rechts eingetragen als Kettenmatrix zu W! und 93 nach (16bb) gebildet. Da 
die hierzu gehörenden Kettenmatrizen Y,!und 9? sind, wird die Matrix 
De Re N en aa: (49) 


zu Null, was im Schema durch Einklammern gekennzeichnet ist. Die Matrix Tß) ist also nicht 
gleich der auf gleicher Schrägzeile stehenden Matrix T3,. Im allgemeinen von Null verschieden 
wird natürlich wieder 


Emo eT .. 2.2.2. (50), 
welche die Reinigung von Vektoren 1. Stufe aus der ersten bzw. dritten Kette bewerkstelligen. 
gel... .2.2.2..0.6l) 


wird wieder rechts eingetragen, dann 9; gebildet und so fortgefahren bis zur Matrix 9, womit 
die fünfte Y-Kette fertig vorliegt. 
Im gesamten Schema stehen nun unter der Hauptdiagonale nur Nullen. Nach Ermittlung 
der Matrizen ©,,und ©,,;ergibt sich die endgültige fünfte X-Kette durch folgende Umrechnung: 
= 8: = Ka ME, = & —- US, —- 55%: | 
x5 Bi &ı ins x S; guy x S arg x3 5 MT & ©, ’ 
B-M—- UN UA MA MC — U — MC | 
Die Gra m sche Matrix der unten stehenden Matrix &7 und der rechts stehenden Matrix Vy 
ist zur Einheitsmatrix geworden und die Biorthonormierung beendet. 


Bestimmung eines zu Vi orthogonalen Vektors 3* in IIIb 
Hierzu geht man von einem von den Vektoren in &7 linear unabhängigen Vektor 3 aus 
und bildet in IIIb 9 3, im Schema ist dafür eine Spalte vorgesehen. Dann ist der Vektor 
+ =: - Mid 
orthogonal in dem gewünschten Sinne, wovon man sich durch Vormultiplikation mit 97’ leicht 
überzeugen kann. 


7. Folgerungen 

In den Abschnitten 5 und 6 wurde die praktische Durchführung der Normierung und die 
Bestimmung eines Vektors 3* gezeigt, der keinen zu einem bestimmten Eigenwert gehörigen 
Eigenvektor oder Hauptvektor enthält. Mit diesen Ergebnissen können wir jetzt auch — wir 
gehen dabei wieder zu der in 4. gebrachten Darstellung der Vektorketten zurück — die in 
der Einleitung angeschnittenen Fragen allgemein beantworten. 

Es sei 

t?, ein Hauptvektor der Stufe o, der Kette 1zum Eigenwert A; und 


Den ein Hauptvektor der Stufe r, der Kette !’ zum Eigenwert 4; 
we a rer rise lies m en; 
2 ee A a 
Mit 27 und 177, bezeichnen wir die nach (27) normierten Vektorsysteme, die die jeweiligen 
reziproken Grundsysteme darstellen. Damit können wir die Entwicklungskoeffizienten bei der 


Entwicklung nach Eigen- und Hauptvektoren allgemein angeben in Erweiterung von (04) bzw. 
(05). In 


wb=ahritahtit tr Fa 4.0. (8) 
können die Entwicklungskoeffizienten aus 
Ben a an nenne . (54) 
T 1 


bestimmt werden. Im Falle einfacher Eigenwerte erhalten wir als Spezialfall (05). j 
Soll schließlich ein Vektor 3* bestimmt werden, der alle zu ,,43, . .., A, mit r< Mm gehören- 
den Eigen- und Hauptvektoren nicht enthält, dann müssen die Bedingungen erfüllt sein: 


Yrzr—=0 NDR ern. (55). 
Die Durchführung kann dabei ohne weiteres aus dem in 6. Gesagten für den Fall mehrerer 


Eigenwerte gefolgert werden. 
Wegen der Durchrechnung von Beispielen sei auf [10] verwiesen. 
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8. Zur Durchführung der Normierung bei komplexen Eigenwerten 


Im Falle eines komplexen Eigenwertes®) A, bei (01) bzw. (02) spalten wir in Real- und 
Imaginärteile auf: 


A=:+m Hetkttiv YmdYrıt tdi: +: (56). 
Mit 8 = WU — £&erhalten wir dann die schon in [11] angegebenen Beziehungen: 
1 
Bir+Nnir=d; Su Rimar, BLr of s20 rr$ gegen hen Re Keen e (57a), 
1 | ee 
—ntrt+ dur =; Erin Dei ee IINRNAFNERIEN BED 


Für das transponierte System haben wir nur Ö durch ®’ und r durch d zu ersetzen. Real- und 
Imaginärteil sind Lösungen des gleichen Problems: 


Brenn. Me 


Für die Matrix ®? ist —n? ein Eigenwert mit der doppelten Vielfachheit wie bei A geworden. 
DB? besitzt also für —n? den doppelten Rangabfall wie Yfür 4, =E&-+ in. 

Das reelle Eigenwertproblem (58) enthält daher 27, linear unabhängige Lösungen, näm- 
lich die Real- und Imaginärteile, die aber durch die Beziehung (57) miteinander verbunden sein 
müssen. 

Der Kürze halber gehen wir hier nur auf den Fall linearer Elementarteiler ein. Wir fassen 
die je ?, linear unabhängigen Lösungen tz, und d,, zusammen zu: 


&n = (Erıs Er Era)» lim Enseshm)- N). 

Entsprechend bilden wir W; und 9),. Auf Grund der Relation (57) gelten die Beziehungen: 
ES 1 ee 
Rd = Drinnen 


Gemäß dem in 2. Gesagten bezeichnen wir die komplexen Matrizen Xrk+:%5 und YR + :IF 
als biorthonormiert, wenn gilt: 


(dEHEIN) (Arksah)ebrn 2. a aloe 


Diese Bedingung können wir so abändern, daß wir es nur mit reellen Größen zu tun haben. 
Normieren wir nämlich, wie nebenstehende Anordnung 
zeigt, so daß 
YR&KR =— Y&r = 6/2 . . . . . (61a) 
DER] IF AR EU RFIRER  T O1 
wird, dann haben wir gerade die Forderung (61) erfüllt. 
Zur Durchführung der Normierung geht man von 9, 
linear unabhängigen Lösungen von (58) aus, bestimmt nach 
(57a) die dazugehörigen Imaginärteile so, daß insgesamt 27, linear unabhän- 
%* x gige Lösungen ermittelt sind und faßt diese zu (X, &,) zusammen. Dann 
eeıe O6, ermitteln wir aus (58), nachdem wir ®? durch 3’? ersetzt haben, ebenfalls 


2 p,linear unabhängige Lösungen, die wir aber, was jederzeit möglich ist, so 
zu p, Vektoren überlagern, daß 


VaE, = 0 un nu eR Jeh DIE OcEEEEE GE 


wird. Damit ist dann auch Y,&r =. 
Nach der in 2. mitgeteilten Methode wird %r mit Wr zu } & normiert. 
Ir = 6 VEN 
Dabei handelt es sich bei Y% um Überlagerungen, die die Relaiion (62) nicht stören. Die 
gleichen Operationen, die auf X} angewandt wurden, müssen nun auch auf X, angewandt wer- 
den, da sonst die Relation (57a) nicht in Ordnung ist. Auch diese Durchführung läßt aber die 
Beziehung (62) bestehen. 
E5 = Kr OTUME DIE SER ERSE NE ) Pa ee EEE 


Ersetzt man nun in (57) Ö durch ®’, dann kann Y, aus 9% ermittelt werden. Diese Matrix ist 
bereits die gesuchte 97. Denn sie genügt den Bedingungen (61a) und (61b), so daß die ge- 
wünschte Normierung beendet ist. 


°) Unter 2, ist ein beliebiger Eigenwert, nicht etwa der kleinste oder größte zu verstehen. 
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Schließlich sei noch auf die Bestimmung eines Vektors 3* eingegangen, der keinen zu 
ı = &+ inundA,— & —ingehörenden Eigenvektor mehr enthält. Dazu geht man von einem 
zu &r und %7 linear unabhängigen Vektor 3 aus und bildet Y% 3 und 9%3. Dann ist der Vektor 


a LI 
3* =; —2(&R, &) 6) 3 


so beschaffen, daß er keine Bestandteile eines zu A, oder A, gehörenden Eigenvektors besitzt. 
Der Vektor ist reell. 
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Hyperbolische Probleme der Gasdynamik mit mehr als 
zwei unabhängigen Veränderlichen 
Von R. Sauer in München 


Die theoretische Gasdynamik hat sich seit 1945 überwiegend mit Fragen befaßt, bei denen mehr als 2 unab- 
hängige Veränderliche auftreten, also mit nicht stationären 2-dimensionalen und mit stationären oder nicht- 
stationären 3-dimensionalen Strömungen. Die vorliegende Arbeit gibt einen systematischen Überblick über 
die hierbei verwendeten mathematischen Methoden. Die Untersuchung beschränkt sich auf reibungsfreie 
Medien und auf Probleme vom hyperbolischen Typus. 

Theory of compressible flow dealt, after 1945, in first line with questions involving more than 2 indepen- 
dent variables (non steady two-dimensional flow, steady and non steady three-dimensional flow). The present 
paper gives a systematic review of the mathematical methods applied for the solution of these questions. The 
review is restricted on nonviscous flow and on problems of hyperbolic type. 

Dans la theorie des &coulements des fluides compressibles on a &tudie, apres 1945, en premier lieu des 
questions oü le nombre des variables independantes esi superieur a 2, (Ecoulements bidimensionnels non per- 
manents, ecoulements tridimensionnels ou permanents ou non). Il est &tabli, dans la note prösente, un resum£ 
syst&matique des methodes mathematiques que l’on a employees pour la solution de ces questions. Le resume 
ne concerne que les fluides non visqueux et les problömes du type hyperbolique. 

Teoperuyeckas AUHAMUKA TA30B 3BAHNMAIach c 1945 ToMa IIPeHMyIIeCcTBeHHO BOIIPOCAMH, 
B KOTOPEIX 4ucHo PurypHMpyIomuXx He3aBHCHMBIX THEePeMEeHHBIX ÖOJIBIIe ABYX, — PAccMATpH- 
BAJIUCb, CHENOBATEJIBHO, HECTAIHOHAPHEIE ABYMepHble HU CTAIHOHAPHEIE HIH HeCTaIHOHAPHBIE 
TpexMepHsie reyenun. Hacrosımas pa6oTa NaeT cHCTeMaTuyeckHü 0030p YIOTPeÖNAeMEIX 
IPH 9TOM MATeMATUYeCKUX METONOB. Mccste1oBaHne OTPAHHYHBAETCH HAEAJIbHBIMH Cpelamu u 
ITPoÖ.IeMaMH THUIePOOAHYecKOrO TUTA. 


Die theoretische Gasdynamik der kompressiblen Medien hat sich seit 1945 großenteils auf 
Probleme mit mehr als zwei unabhängigen Veränderlichen konzentriert, also auf das Studium 
der zweidimensionalen nichtstationären und der dreidimensionalen stationären oder nicht- 
stationären Strömungen. Da in den Lehrbüchern der Gasdynamik (vgl. z. B. die Zitate [2], 
[3], [7b] am Schluß dieses Aufsatzes) diese Probleme noch nicht eingehend behandelt sind, 
wird in dem vorliegenden Bericht!) ein systematischer Überblick über dieses neuere Teilgebiet 


1) Der Bericht bildete den Gegenstand eines Vortrags bei der GaMM-Tagung in Braunschweig am 
5. Juni 1952. Eine ausführliche Behandlung in Buchform wird von meinem Mitarbeiter Herrn Dorfner vor- 


bereitet. 


; R Math. Mech. 
332 Sauer, Hyperbolische Probleme der Gasdynamik Bd. 33 Nr. ORT OkE Nov. 1058 


der Gasdynamik gegeben. Wir beschränken uns dabei auf solche Strömungsvorgänge, bei 
denen die Wärmeleitung und innere Reibung vernachlässigt werden können, und außerdem 
auf Probleme, bei denen die zugrundeliegenden Differentialgleichungen vom hyperbolischen 
Typus sind. Die weiten Gebiete der Grenzschichtfragen und der stationären Unterschall- 
strömungen bleiben also außer Betracht. Ein ausführliches Verzeichnis der einschlägigen Lite- 
ratur findet sich im Leitartikel des Aprilheftes der Applied Mechanics Reviews; 
vgl. [7a]. 

Der folgende Bericht ist nicht nach dem materiellen Inhalt der in Frage stehenden Pro- 
bleme, sondern nach den zu ihrer Lösung verwendeten mathematischen Methoden gegliedert. 
Die linearen Methoden (Abschnitt I) werden ausführlicher, die nichtlinearen Methoden (Ab- 
schnitt II) nur kurz besprochen. 


I. Lineare Methoden 


Mittels der linearen Methoden kann man die Strömung um ‚schlanke Körper“ (z. B. 
Drehkörper) und um ‚„‚dünne Körper‘ (Tragflügel) sowie um Kombinationen solcher Körper 
(Untersuchungen über die Wechselwirkung zwischen Rumpf und Flügel) behandeln. Hierbei 
ist der umströmte Körper entweder in Ruhe (stationäre Probleme) oder führt periodische 
Schwingungs- oder nichtperiodische Flatterbewegungen aus (nichtstationäre Probleme). 

Durch die Linearisierung werden die aus den Erhaltungssätzen der Masse, des Impulses 
und der Energie sich ergebenden nichtlinearen (quasilinearen) Different ialgleıc tungen im 
stationären Fall auf die lineare Wellengleichung 


L[vo] = Po, — do =0, ß=ctga, =JM2—1 =const SoRL ba 1 


reduziert. wist das Geschwindigkeitspotential p oder eine Ableitung von g. Die Veränderliche 
& ist die Ortskoordinate in der Richtung der Grundströmung (ungestörte Parallelströmung mit 


der konstanten Mach-Zahl M, = _ > 1 und dem konstanten Mach-Winkel a,); der 


Laplacesche Differentialausdruck Aw = w,, + @,, enthält die zweiten Ableitungen nach 
den übrigen Ortskoordinaten %, 2. 

Im nichtstationären Fall kommt man, wenn die Grundströmung durch eine Galilei- 
Transformation auf Ruhe transfor- 
miert wird, auf die lineare Wellenglei- 
chung 


1 
— 0 — do —=0 nzs (2); 
a2, 
t = Zeitkoordinate.. Aw = 0 + Oyy 
bzw. = Oo + Oyy + Bo: 


Die mathematische Aufgabe besteht 
darin, die Wellengleichung (1) bzw. 
(2) unter den durch die Gestalt und Be- 
wegung des umströmten Körpers be- 
stimmten Anfangs- und Randbedingun- 
gen zu integrieren. Es handelt sich hier- 
bei nicht um das Cauchysche An- 
fangswertproblem (vgl. [8]), sondern um 
verwickeltere Anfangs- und Randbe- 
dingungen. 

Wir erläutern die linearen Metho- 
den am Beispiel der stationären Über- 
schallströmung um den ruhenden Trag- 
flügel endlicher Breite (Bild 1); vgl. [4], 

Bild 1 [5]; [6]. 
1.Anfangs- und Randbedingungen 

Die Anfangs- und Randbedingungen werden folgendermaßen präzisiert: 

a) Der Flügel wird ersetzt durch seine Projektion (F) in der x, y-Ebene. Die Bedingung 
tangentialer Strömung an der Körperoberfläche führt dann zur Vorgabe der Geschwindigkeits- 
komponente 9, bzw. , auf der oberen und unteren Seite des Gebiets (F). 

b) Die Grundströmung bleibt ungestört bis zu der vom vorderen ‚‚Überschallteil‘“ AA’ 
des Randes von (F) ausgehenden stromabwärts laufenden charakteristischen Fläche (vordere 
Mac h-Fläche); vgl. das ‚‚Prinzip der verbotenen Signale“ in [1]. Die Zerlegung des Randes 


Z. angew. Math. Mech. P : 
Ba. 38 Nr. 10/11 Okt. /Nov. 1953 Sauer, Hyperbolische Probleme der Gasdynamik 333 


von (F) in die Unterschallteile ACB, A4’C’B’ und die Überschallteile 4A’, BB’ ist aus Bild 1 
unmittelbar verständlich; bei ‚‚nichtovalen‘‘ Flügelformen ergeben sich entsprechende andere 
Zerlegungen. 

c) Die Verdichtungsstöße und etwaige Verdünnungsfächer der strengen Theorie werden 
bei den linearen Methoden auf stromabwärts laufenden charakteristischen Flächen, z. B. der 
in b) eingeführten vorderen M a c h-Fläche, konzentriert ; auf diesen Flächen bleibt die tangen- 
tiale Geschwindigkeitskomponente stetig, die übrigen Ableitungen von aber weisen gewisse 
Unstetigkeiten (Sprünge oder Unendlichkeitsstellen) auf; vgl. das ‚‚Prinzip der konzentrierten 
Wirkung‘ in [1]. 

d) Wie in der Aerodynamik der inkompressiblen Medien wird dem Flügel eine Wirbel- 
schicht angeheftet und diese Wirbelschicht durch ihre Projektion ( W) in der x, y-Ebene ersetzt. 
(W) ist begrenzt von den beiden zur x-Achse parallelen Tangenten der Flügelprojektion (F). 
Das Geschwindigkeitspotential 9 und Ableitungen von p können in (F) sowie in (W) gewisse 
Unstetigkeiten besitzen. In (F) kann 9, nur auf Teilen des Randes unendlich werden. In (W) 
kann, da der Druck p stetig ist und die Strömung tangential zur Wirbelschicht verläuft, nur o, 
unstetig sein. 

e) Die Vernachlässigung der Reibungskräfte hat zur Folge, daß zur eindeutigen Fest- 
legung der Strömung ein ‚‚Auswahlprinzip‘‘ als zusätzliche Bedingung hinzugefügt werden muß. 
Es besteht in der Forderung, aus der Gesamtheit der mathematisch möglichen Lösungen eine 
Lösung mit ‚„‚möglichst schwachen Singularitäten‘‘ auszuwählen. Im Spezialfall der zwei- 
dimensionalen inkompressiblen Strömung ist eine Auswahlforderung dieser Art die Kutta- 
Joukowski-Bedingung des Endlichbleibens der Geschwindigkeit (glatter Abfluß) an der 
Hinterkante. Bei dem hier in Frage stehenden Tragflügelproblem stellen wir die Bedingung, 
daß p, auf den ‚‚Unterschallteilen‘‘ CB, C’B’ der Hinterkante stetig ist; dann nämlich ist dort 
auch der Druck 9 stetig. 

Zur Lösung der so präzisierten Anfangs-Randwert-Aufgabe stehen im wesentlichen drei 
lineare Methoden zur Verfügung, nämlich die Methode der Grundlösung (Ziff. 2), die Methode 
der homogenen Lösungen (Ziff. 3) und die Operatorenmethoden (Ziff. 4). 


2.Methode der Grundlösung 
Als Grundlösung (,,Überschall- Quelle‘‘) der Wellengleichung (1) wird die Funktion 


m Innern 
Grasyasense) Kae DR neh des Mac h-Kegels (3) 
der vom Aufpunkt (£, n, &) stromabwärts läuft (Bild 2), bezeichnet; dabei ist 
ea ee a le nd ee 
Die Grundlösung @ ist im Innern des Ma c h-Kegels reell und wird auf seiner Oberfläche un- 
endlich. 
v.Karman, Tsien und Andere 
(vgl. etwa [7b]) haben I1-parametrige An- 
ordnungen dieser Überschall- Quellen und 
der durch Differentiation nach { entstehen- 
den Überschall-Dipole zur Darstellung der 
Strömung um achsial oder schief angebla- 


Anström- 


richtung 


Bild 2 Bild 3 


sene Drehkörper benützt. Die in Ziff. 1 formulierte Tragflügelaufgabe läßt sich durch 2-para- 
metrige Anordnungen solcher Quellen und Dipole folgendermaßen behandeln: 
Die vom Überschallteil AA’ der Vorderkante der Tragflügelprojektion (F) (vgl. Bild 1) 
ausgehende vordere Mac h-Fläche und der von einem dahinter liegenden Aufpunkt (&,n,£) 
24 
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stromaufwärts laufende Mac h-Kegel schließen zu beiden Seiten der z, y-Ebene ein Raum- 
gebiet (V) bzw. (V’) ein (Bild 3). Durch Anwendung der Greenschen Formel auf (V) und 
(V’) wieinder Hadamardschen Theorie des Cauch yschen Anfangswertproblems (vgl. 
[5], [6], [8]) erhält man für den Wert der gesuchten Funktion im Aufpunkt 


on d=|[ [I ge) wo) 2) ds. dyıomdtsunls.td) 


Hierbei ist das Integral über das (in Bild 1 schraffierte) Gebiet (y) der z, y-Ebene zu er- 
strecken, an dem die Raumgebiete (V), (V’) aneinandergrenzen. Das Integral ist divergent und 
muß, wie das Hakenzeichen vor dem Integral andeuten soll, durch seinen Hadamardschen 
Hauptwert (‚‚partie finie‘‘) ersetzt werden. Der Hadamardsche Hauptwert ist eine mit 
dem bekannten Cauch y schen Hauptwert verwandte Begriffsbildung der Analysis. 

Es handelt sich nun darum, die in die Darstellungsformel (4) eingehenden Werte von © 
und ; auf der oberen Seite von (y) und von w’ und x auf der unteren Seite von (y) aus den 
Anfangs- und Randbedingungen des Problems zu ermitteln. Hierfür zerlegen wir die Strömung 
in üblicher Weise in einen symmetrischen und antisymmetrischen Bestandteil: 

Bei der symmetrischen Strömung 9 (8, %,2) = p(&, y, — 2) ist 9— 9’ = O in der ganzen 
&, y-Ebene und 9, — 9; =0in der z, y-Ebene außerhalb (F). Bei vorgegebenem 9, — ; in 
(F), also vorgegebener Geometrie des Tragflügels, liefert die Darstellungsformel (4) mit » = 
die gesuchte Strömung im ganzen Feld. 

Bei der antisymmetrischen Strömung g(z, %, 2) = — (2, Y — 2) ist 94, — fu = in 
der ganzen x, y-Ebene und 9, — 9 = 0 in der x, y-Ebene außerhalb (F). Bei vorgegebenem 
9x — Yr ın (F), also vorgegebener Auftriebsverteilung des Tragflügels, liefert die Darstellungs- 
formel (4) mit = 9, die gesuchte Strömung im ganzen Feld. Wesentlich schwieriger ist es 
hier, bei vorgegebenem 9, in (F), also vorgegebener Geometrie des Flügels, die Strömung zu 
ermitteln; denn in diesem Fall gehen in die Darstellungsformel (4) auch die von vornherein 
nicht bekannten Daten in den Teilgebieten von (y) ein, die nicht zu (F) gehören; vgl. in Bild 1 
das Gebiet APQ. Man muß dann diese Daten, unter Umständen mit Zuhilfenahme des Aus- 
wahlprinzips, erst berechnen. Im Falle einer ‚‚ovalen‘‘ Flügelform, d. h. einer wie in Bild 1 aus 
zwei Unter- und zwei Überschallteilen bestehenden Randes der Flügelprojektion (F), führt dies 
aufeine Abel sche Integralgleichung; vgl. [4], [5], [6]. 

Vorteile der Methode der Grundlösung; Anwendbarkeit auf sehr allgemeine Flügelfor- 
men; Aussagen über Existenz, Eindeutigkeit und Singularitäten der gesuchten Lösungen; Aus- 
sagen über das Problem der Umkehrung der Strömungsrichtung. 

Durch die neuerdings vonL. Sch war zentwickelte Theorie der ‚‚uneigentlichen Funk- 
tionen‘ (vgl. [9]) läßt sich die Methode der Grundlösung in den Rahmen eines sehr umfassenden 
Kalküls einordnen. 


3. Methode derhomogenenLösungen 


Ist das Geschwindigkeitspotential eine in z, %, z homogene Funktion vom Grade n, so 
führt der Ansatz 


« . p(%, Y; 2) = m"Öd(n,?), = (= 


auf eine Funktion ®(n,£) mit nur zwei unabhängi- 
gen Veränderlichen. Nach geeigneter Transforma- 
tion dieser Veränderlichen folgt aus der Wellen- 
gleichung (1) für die Ableitungen von 9(z, %, 2) 
von der Ordnung n in der Form Q(n,£) eine sehr 
einfache Differentialgleichung, nämlich 


A2(n*,&*) =0 im Innern \des Mach-Kreises 
PERe) r P+öÖ=tg’a 
—— —=0(0 im Auß 
DO nl Bilda; (5) 


Für n=1 ergibt sich als Spezialfall die Buse- 
mannsche kegelsymmetrische Überschallströ- 
mung; vgl. z.B. [7b]. Wie in diesem Sonderfall liegt 
auch bein > lim Innern des Mach - Kreises eine Randwertaufgabe vor, die mit den Hiilfs- 
mitteln der konformen Abbildung gelöst werden kann. 

Vorteile der Methode der homogenen Lösungen: Bequeme Anwendbarkeit auf praktisch 
wichtige Flügelformen (Dreieck, Rechteck, Pfeil, Trapez) und auf geeignete nichtstationäre 
Randbedingungen (,‚pseudostationäre Strömungen‘‘, rollender Dreieckflügel). 


Bild 4 


OÖ 
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- Statt wie im Falle des Rechtecks-, Pfeil- oder Trapezflügels nur endlich viele homogene 
Lösungen zu superponieren, kann man auch 1- oder 2-parametrige Lösungsscharen überlagern 
und dadurch sehr allgemeine Flügelformen erfassen. 


4.Operatorenmethoden 
Die Operatorenmethoden führen die Wellengleichung (1) bzw. (2) durch eineFourier- 
Transformation oder eine Laplac e-Transformation auf eine Differentialgleichung mit einer 
um 1 verminderten Anzahl der unabhängigen Veränderlichen zurück. Die Fourier-Trans- 
formation wurde insbesondere von v. Karman verwendet (‚Methode der akustischen Ana- 
logie“), Nähere Ausführungen zu den Operatorenmethoden finden sich in [1], [4] und [5]. 


5.Korrektionen dderlinearen Lösungen 
Korrektionen für die von der linearen Theorie gelieferten Näherungslösungen werden in 
dem Artikel von Goldstein eingehend erörtert; vgl. [4]. Einerseits werden Entwicklungen 
nach einem Dickenparameter aufgestellt, andrerseits wird das Verhalten der Lösungen in der 
Umgebung singulärer Ma c h-Linien und in großem Querabstand vom Körper untersucht. 


II. Nichtlineare Methoden 
Wir wenden uns jetzt zu den strengen, nicht linearisierenden Methoden der Gasdynamik. 
Sie beruhen auf der Charakteristikentheorie der Systeme quasilinearer Differentialgleichungen 
erster Ordnung vom hyperbolischen Typus. 


1. Halblineare Methoden 


Während die linearen Methoden von der ungestörten Parallelströmung als Grundströ- 
mung ausgehen und die Abweichungen von der Parallelströmung nur linear berücksichtigen, 
benützen die halblinearen Methoden geeignete allgemeinere Grundströmungen, z. B. die dreh- 
symmetrische stationäre Überschallströmung um einen Drehkörper von vorgegebener Form. 
Dieser Grundströmung, die sich mit einem der üblichen Charakteristikenverfahren für Probleme 
mit zwei unabhängigen Veränderlichen ermitteln läßt, wird eine linearisierte Zusatzströmung 
überlagert und die Anzahl der unabhängigen Veränderlichen für diese Zusatzströmung durch 
Separationsansätze ebenfalls auf zwei erniedrigt. Auf diese Weise wurde insbesondere die 
stationäre Überschallströmung um Drehkörper unter kleinem Anstellwinkel berechnet, die 
Methode ist aber auch für nichtstationäre Nachbarströmungen (z. B. pendelnde Drehkörper) 
verwendbar; vgl. z. B. [7b]. 

Sehr eingehend wurde von verschiedenen Autoren der Spezialfall des Drehkegels unter- 
sucht mit Berücksichtigung der Tatsache, daß im Falle eines nicht verschwindenden Anstell- 
winkels die Strömung hinter der Kopfwelle nicht mehr wirbelfrei ist. Ausführliche Tabellen für 
die Kegelströmung hat K o pa 1 [10] aufgestellt. 


2.Strenge Methoden 
Die strengen Charakteristikenmethoden erfordern bei mehr als zwei unabhängigen Ver- 
änderlichen einen Rechenaufwand, der ohne den Einsatz moderner programmgesteuerter 
Rechengeräte kaum zu bewältigen ist. 


3.GrenzfälleM>1und M>o 


In den Grenzfällen M— 1 (transonic flow) und M— © (hypersonic flow) versagen die 
Näherungsgleichungen der linearen Theorie; sie sind durch nichtlineare Näherungsgleichungen 
zu ersetzen, z.B. im Fall M—1 durch die Trico mische Gleichung. Diese Näherungs- 

leichungen wurden von v. Karman, Tsien und anderen Autoren zur Aufstellung von 
Ähnlichkeitsgesetzen verwendet, die für Modellversuche von großer praktischer Bedeutung 
sind. Neuerdings hat auch Oswatitsch Ähnlichkeitsgesetze entwickelt, wobei er insbe- 
sondere die Vorgänge im Verdichtungsstoß berücksichtigt. Literaturangaben siehe [7a]. 


4.Konvergenzfragen und Fehlerabschätzung 

Die Frage der Konvergenz der in der Gasdynamik gebräuchlichen Gitterkonstruktionen 
bzw. Differenzenverfahren der Charakteristikentheorie wird für Probleme mit zwei unab- 
hängigen Veränderlichen im Buche von Conrant-Friedrichs [2] angeschnitten. In 
einer neueren Arbeit vonR.CourantundP.Lax[11] wird für Systeme quasilinearer Diffe- 
rentialgleichungen vom hyperbolischen Typus mit zwei unabhängigen Veränderlichen ein ite- 
ratives Charakteristikenverfahren entwickelt. An praktisch brauchbaren Fehlerabschätzungen 
sowie an Verfeinerungen der Differenzenverfahren im Sinne der Methode von Runge und 
Kutta wird gegenwärtig von meinem Mitarbeiter Herrn Bauhuber gearbeitet. 
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Der vorliegende Bericht kann naturgemäß keinen Anspruch auf Vollständigkeit erheben. 
Insbesondere wurde auf eine Erörterung der Probleme mit Typenwechsel (Überschallströ- 
mungen um Körper mit nicht anliegender Kopfwelle, Unterschallströmungen um Körper mit 
anliegendem Überschallbereich) verzichtet. 
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Eingegangen am 20. Juli 1952. 


Die Energieversorgung der Wandturbulenz ”) 


Von H. Reichardt in Göttingen 
(Aus dem Max Planck-Institut für Strömungsforschung, Göttingen) 


Es werden die Gesetzmäßigkeiten der Energieversorgung der Wandturbulenz am Beispiel der ebenen Kanal- 
strömung behandelt. Ähnlich wie bei der Geschwindigkeitsverteilung kann man auch beim Energieumsatz 
zwischen einem Wandgesetz und einem Mittengesetz unterscheiden. 


As an example of ihe laws, governing the supply of energy of the turbulence near the wall, the plane 
channel flow is treated. Similarly to the distribution of velocity, one law for the region near the wall and one 
law for the middle region can be distinguished for the transformation of energy. 


Les lois gouvernant l’approvisionnement d’Energie de la turbulence murale sont traitees & l’ewemple du 
courant de canal plan. Analoguement & la distribution de la vitesse on peut distinguer & la transformation de 
l’energie une loi pour la region murale et une telle pour la region du milieu. 


JAKOHOMEPHOCTH 3HEPTOCHAÖSKEHHA TYPÖYJIEHTHOCTH Y CTEHKH TPAKTYITCH Ha IIPHMepe 
ILIOCKOCTHOTO KAHANIOBOTO TeyeHugn. Ilo10o0Ho cyyam pacıpenereHug CKOPocTeä, PasıımyarmT, 
Ipu NpeoÖpasoBaHuMm 3HePTUM, MeikAy 3AKOHOM IOTPAHHYHOH OÖJTACTH H 3AKOHOM CPeAHHHOÜ 
ODNACTH. 


1. Einleitung 


Bei den bisher durchgeführten Untersuchungen des Energiegleichgewichtes in turbulenten 
Reibungsschichten hat man sich vornehmlich mit den wandfernen Bezirken beschäftigt, in 
denen die Zähigkeitskräfte im Vergleich zu den Trägheitskräften vernachlässigbar klein sind. 
Die Energiebilanz in unmittelbarer Wandnähe, wo die Zähigkeitskräfte eine wesentliche Rolle 
spielen, ist anscheinend noch nicht behandelt worden. 

Im folgenden werden nun die Gesetzmäßigkeiten der Energieversorgung der Wandtur- 
bulenz behandelt und zwar unter besonderer Berücksichtigung des bemerkenswert hohen Ener- 
gieumsatzes in den wandnahen Schichten. Diese Betrachtungen werden für die ausgebildete 
Strömung durch einen flach rechteckigen Kanal durchgeführt. 


2. Energiegewinn aus der Hauptbewegung; Definition von M* 


In diesem relativ einfachen Beispiel einer Hauptbewegung besteht der Ausdruck für die 
Energieversorgung nur aus einem einzigen Glied. Die je Volumen- und Zeiteinheit von der 


*) Auszugsweise vorgetragen auf der GaMM-Tagung in Braunschweig am 6. Juni 1952. 
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Hauptbewegung geleistete Arbeit beträgt z, du/dy, wenn z, die turbulente Schubspannung und 
U die mittlere Strömungsgeschwindigkeit im Wandabstand % bedeutet (dieser Ausdruck werde 
mit M bezeichnet). Dafür läßt sich auch schreiben 


a A En reiht (1), 


wenn 7,, die molekulare Schubspannung und u den Zähigkeitskoeffizienten bedeutet. 

Formel (1) zeigt, daß die Turbulenzerregung durch das Zusammenwirken der 
turbulenten und der zähen Spannungen bedingt ist. 

Wir wollen nun eine dimensionslose Größe M* dadurch definieren, daß wir M mit u/12 
multiplizieren (7, sei die Wandschubspannung). Dann ist also M* = 1,1,/r2. Dafür läßt sich 


auch schreiben 
A/u = 
* ne _ _ ‘ 
Ska [aursp 7, be 2 


Hier ist A die turbulente Austauschgröße, u der Zähigkeitskoeffizient; A/wist identisch mit dem 
Verhältnis 7,/7„. M*ist also im wesentlichen abhängig vom Verhältnis A/u, das eine Funktion 
des dimensionslosen Wandabstandes n = yu*/v darstellt (u* — Yry/o — Schubspannungsge- 
schwindigkeit; v = u/o = kinematische Zähigkeit). 

Unter der Voraussetzung, daß näherungsweise 7 — r, gesetzt werden kann, läßt sich über 
die Größe M* auch ohne nähere Kenntnis der Funktion A/u(n) folgendes aussagen: 

1. M* hat das Maximum in dem Wandabstand n, wo A/u = list, wo also die turbulente 
Reibung ebenso groß ist wie die zähe Reibung. 

2. Das Maximum von M* beträgt 1/4. 

Da die n-Stelle, wo A = wist, bei hohen Re-Zahlen sehr nahe an der Wand liegt, so ist 
die vorausgesetzte Bedingung, daß 7 = r, sei, für nicht zu kleine Re-Zahlen praktisch erfüllt. 


3. Das Energiegleichgewicht. 


Bevor wir die Eigenschaften von M* im einzelnen erörtern, wollen wir die Energiebilanz 
der Nebenbewegung kurz diskutieren und zwar für jede Schwankungskomponente einzeln. 
Wenn man die Energiegleichungen dieser drei Komponenten für die stationäre ausgebildete 
Kanalströmung dimensionslos macht — und zwar ebenfalls durch Multiplikation mit u/d — 


und für fin 
Ru, ou, 2 jr u? 
re la) "D® 
usw. setzt, so erhält man die Gleichungen (3) bis (5): 
dv,W  u,0p' de u? o9u,\? ou,\? ou,\® E 
RR, | ee EI RER INU HN 
r 2 dn Tu0E  2dn? | 3 Bi) El (8) 
dv: v,0p' dw ov,\2 0v,\? av,\? 
£ - ee / : TITEL (A): 
2 dn on dns ) a Ze) (4) 


dv,w w,op' , d?w? ow, ig ow,\? | ow,\? 
2Adn Toß ä le eo) 0% 
Hier sind & = zu*/v; n = yu*/v und £ = zu*/v die dimensionslosen Raumkoordinaten und 
u, = u/u*; v, = v//u* und w, = w‘/u* die Komponenten der dimensionslosen Nebenbewegung 
in den Richtungen &, y und 2. 
Die Größe M*, die den Energiegewinn der Nebenbewegung darstellt, steht nur in der 
Gl. (3) für die Energie der Längsschwankungskomponente «. Die übrigen Glieder auf den 
linken Seiten der Gleichungen (3) bis (5) beschreiben die Energiediffusiont) der Nebenbewegung, 
wobei sich die sechs Glieder mit v, und p’ auf die durch die Schwankungen selbst hervorgerufene 
Diffusion beziehen, während die zweiten Ableitungen der drei Schwankungsquadrate einen 
durch die Zähigkeitsspannungen verursachten Transport von T urbulenzenergie darstellen. 
Auf den rechten Seiten der Gleichungen (3) bis (5) stehen die stets positiven Glieder der 
Dissipation der Schwankungsbewegungen. Allerdings handelt es sich hier nicht um die örtliche 
Dissipation, die in den gebräuchlichen, von O.Reynol ds entwickelten Energiegleichungen 
steht. Die Bezeichnung ‚‚Dissipationsglieder‘ soll lediglich zum Ausdruck bringen, daß diese 
Glieder die Gesamt dissipation der Nebenbewegung ergeben, wenn sie über den Raum 
bzw. über n integriert werden. 


1) Die Diffusionsglieder sind dadurch gekennzeichnet, daß ihre Raumintegrale verschwinden. 
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Um auf den rechten Seiten von (3) bis (5) die Ausdrücke für die Anteile der einzelnen Schwankungs- 
komponenten an der örtlichen Dissipation zu erhalten, müßte man — wie sich leicht zeigen läßt — 


s 5 E f F e d ou, 
das Gleichungssystem um gewisse Diffussionsglieder!) erweitern und zwar (3) beiderseits um dr v, DE 


(4) beiderseits un _4 (v, 22 und (5) beiderseits um 2 (o, =) (dus wäre gleichbedeutend mit einer Er- 
292 


dn on, 
weiterung der Summe der drei Gleichungen beiderseits um das Diffusionsglied 


v; 


dn? 


) . Diese formalen Ma- 


nipulationen hätten jedoch nur dann einen praktischen Wert, wenn weitere, über den Inhalt von (3) bis 
(5) hinausgehende Angaben über die lokale Dissipation möglich wären. Solange aber noch keine verbind- 
lichen Aussagen über die örtliche Dissipation der nichtisotropen Turbulenz existieren, besteht auch keine 
Veranlassung, die örtliche Dissipation in den Energiegleichungen zur Darstellung zu bringen, wenn die 
Formeln dadurch komplizierter werden. Daher werden wir die Energiegleichungen in der oben mit- 
geteilten Form (Gl. (3) bis (5)) verwenden. 


Der Energietransport durch Zähigkeitskräfte ist für die Wandnähe von Bedeutung, wo 
die Querschwankungen v’ praktisch verschwinden. Beispielsweise wird die Dissipation der 
«-Schwankungen in der zähen Wandschicht ausschließlich durch eine solche Energiezufuhr 
gedeckt. Für diese Grenzzone verschwinden in Gl. (3) die meisten Glieder und es bleibt übrig 


au _ (du) 
2dn? \dn 


links steht die Energiezufuhr durch Zähigkeitskräfte, rechts steht die Dissipation) ?). Diese 
Gleichung ist in der Tat erfüllt, da die w-Schwankungen in der zähen Wandschicht in jedem 
Augenblick dem Wandabstand n proportional sind. (Dasselbe gilt für die w-Schwankungen.) 

Wenn man die Vorgänge, die sich in größerer Entfernung von der Wand abspielen, be- 
trachten will, empfiehlt es sich, die Koordinaten &, n und Z durch die dimensionslosen Wand- 
abstände z/r, y/r und z/r zu ersetzen (wobei r die halbe Weite des Kanales bedeutet). Nach Ein- 
führung dieser neuen Koordinaten geht beispielsweise Gl. (3) über in 


vo, ww 1 Ru 1 (/ Qu, \® ou,\? al 
M*n. — DER = Fr FRAGEN) 
ne dy/r  Tuoa/r A n, 2 d(y/r)? lo) - (au) | o2/r} | Er 
Hier ist zu beachten, daß die Größe‘ M* mit dem Faktor n, = ru*/» behaftet ist, während vor 
den Dissipationsgliedern der rechten Seite 1/n, steht. 


4. Das Wandgesetz der Energiegewinnung 


Wir wollen nun die Diskussion der Eigenschaften von M* fortsetzen (s. Gl (2)). Diese 
Größe hängt von A/u ab und A/u ist eine Funktion des dimensionslosen Wandabstandes 
n = yu*/v. Der genaue Verlauf dieser Funktion ist noch unbekannt. A/u läßt sich aber nähe- 
rungsweise darstellen. Von den vom Verfasser erprobten Funktionen?) hat sich die folgende 
bisher am besten bewährt: 


Alu = nn — Tan (n/a) Tansta)| = 3 0) 


(hier istx = 0,4 und 7, = 7,15). Diese Formel ermöglicht es, sowohl die Geschwindigkeitsver- 
teilung in Wandnähe als auch den turbulenten Wärmeübergang zu berechnen und zwar in guter 
Übereinstimmung mit den Experimenten. 


Formel (6) ist in Bild 1 dargestellt; einzelne Zahlenwerte stehen in Tabelle 1. Demnach 
ist A/uin unmittelbarer Wandnähe (etwa bis 7 = 2) praktisch Null (hier ist das Gebiet der soge- 
nannten lJaminaren Wandschicht), steigt dann mit wachsendem Wandabstand allmählich an, 
um bei größeren n-Abständen (etwa bei 7 = 40) in eine Gerade mit der Neigung 0,4/1 überzu- 
gehen. (Der geradlinige Verlauf ist durch das logarithmische Geschwindigkeitsprofil bedingt). 

Bein= 10,7 ist A/u =1bzw. A= u. Die Stelle n = 10,7 stellt also die Grenze dar 
zwischen dem schmalen wandnahen Gebiet mit überwiegender Zähigkeitsreibung und dem aus- 
gedehnten wandfernen Bereich mit überwiegender Turbulenzreibung. Das ersieht man auch 
aus dem in Bild 1 dargestellten Verlauf des Verhältnisses 7,/7,.,, das sich für große Wandab- 
stände dem Grenzwert 1 asymptotisch nähert (dort wo A/u — list, ist T,/Tyes — 0,5.) 


?) Eine entsprechende Beziehung erhält man aus (5) für die ebenfalls parallel zur Wand verlaufende 
Schwankungskomponente w’. | 
.. »)H.Reichardt, Vollständige Darstellung der turbulenten Geschwindigkeitsverteilung in glatten 
Leitungen. ZaMM 31 (1951), S. 208/219. — H.Reichardt, Die Grundlagen des turbulenten Wärmeüber- 
ganges. Arch. ges. Wärmetechnik 1951, S. 129/142. 
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Bei Verwendung des dargestellten Tabelle l 
A/u ergibt Gl. (2), wenn T=r, gesetzt 
wird, den ebenfalls in Bild 1 aufgezeich- in 2% ur zus 
r ; & 

neten Verlauf von M*. Wie bereits gesagt, 20 0,0008 0,0008 0.0008 
liegt das Maximum von M* an der Stelle, 30 0.006 0.006 0.006 
wo die turbulente Reibung gleich der zähen 4,5 0,04 0,039 0,037 

6,0 

9,0 

7 


Reibung ist, und das Maximum von M* be- 0,13 0,115 0,102 
0,58 0,367 0,232 


ee : < 
Hi 0,25. en Zahlenwerte von M 10, 1.00 0.500 0'250 
stehen in Tabelle 1. 11 1,13 0,531 0,249 
Da 7 = r, vorausgesetzt wurde, und da 12 428 0,576 DR 

it die Funktion M* nur von n und nicht - vr we Ya 
somit die Funktion nn nic 18 3,46 0,776 0,174 
von y/r abhängt, kann man hier von einem 21 4,62 0,822 0,146 
„‚Wandgesetz‘‘ der Energiegewinnung aus 2% Bne0 are DS 
der Hauptbewegung sprechen. 48 15.4 0.939 0.057 


r 


025 


0,15 


005 


Bild 2 
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Bild 2 zeigt außer der Funktion M* die vom Verfasser gemessene Verteilung der Schwan- 


kungskomponente Ve 4) sowie das dazugehörige Diffusionsglied d? u7/2dn?, dessen Bedeutung 
im vorigen Kapitel dargelegt wurde. Diese Angaben über die Energiediffusion durch Zähig- 
keitskräfte sind allerdings nur sehr roh. Da man von einer fehlerhaften experimentellen Funk- 
tion nicht die zweite Ableitung bilden kann, ist hier umgekehrt die zweite Ableitung so gewählt 
worden, daß einerseits die theoretisch erforderliche Grenzbedingung an der Wand erfüllt wird 
(linearer Anstieg von Yaz) und daß andererseits die durch Integration gewonnene Vw-Kurve 
einigermaßen durch die Meßpunkte geht. Das Maximum von 277 liegt bei n = 22, ist also etwa 
doppelt so weit von der Wand entfernt wie das Maximum von M*5). 

M* + d?u2/2 dn? ist die Energie, die für die Deckung der Dissipation der w,-Schwan- 
kungen in Wandnähe zur Verfügung steht, wo die Glieder mit «, und 9° entfallen oder noch 
schwach sind. (Siehe Gl. (3)). Wie man auf dem Bilde sieht, erfolgt die Energieversorgung der 
Wandschicht (n< 8) auf Kosten der Energiebelieferung des weiter entfernt liegenden Ge- 
bietes. 


5. Das Mittengesetz der Energiegewinnung 


Während das Verhältnis A/u in kleinen Wandabständen y/r nur von n = yu*/v abhängt, 
ist der Ausdruck A/u n, in größeren Wandabständen nur eine Funktion von y/r (diese letztere 
Aussage ist die Grundlage für das sogenannte Mittengesetz der Geschwindigkeitsverteilung®)). 


400 


700 


0 0.05 01 0,75 02 0.25 03 
—— Yr 
Bild 3 


Wenn man nun in der Definitionsgleichung für M* (Gl. (2)) für größere Wandabstände, bei 
denen 1 gegenüber A/u vernachlässigbar ist, von A/u auf A/u n, übergeht, so erhält man 


(7%)? 
m. 
u Al 


Hier ist außer A/u n, auch 1/r, nur von y/r abhängig. Die Formel (2a) besagt also, daß M* „, in 
größeren Wandabständen eine von der Re-Zahl unabhängige Funktion des Wandabstandes y/r 


*) H. Reichardt, Messungen turbulenter Schwankungen, Z. Naturwissenschaften Bd. 26, 1938, 
S. 404/408. (In dieser Arbeit wurde bereits die Vermutung ausgesprochen, daß die Anfachung der Turbulenz 
im Übergangsgebiet zur Laminarschicht erfolge.) ; 


ER Sic se Messungen von Lauferliegt das Maximum von Yu? bei n=17(J.Laufer, NACA 


6) Siehe Fußnote 3 Seite 338. 


iterel23% 
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darstellt. Die Funktion M* n, beschreibt sozusagen das ‚‚Mittengesetz“ der Energiegewinnung 
aus der Hauptbewegung. 

Dieser Tatbestand kommt in Bild 3 zum Ausdruck. Hier ist M*n, über y/r aufgetragen 
und zwar für die Werte n, = 400; 800 und 1600 (das entspricht etwa den Re-Zahlen 13 000; 
30 000 und 63 000) (s. Tabelle 2). Wie man 


sieht, ist M*n, — von den kleinen Wand- Tabelle 2 
abständen y/r< 0,1 abgesehen — für die ee 
verschiedenen Re-Zahlen eine einheitliche ®  yir A 
Funktion von y/r. en ei De a Aa 
In Wandnähe dagegen fallen die Kur- 0,002 0,22 4,1 56,0 
ven auseinander. Das Maximum von 0,004 1,20 26,3 270 
M* n, steigt mit wachsender Re-Zahl und 0,006 6,00 76,2 Sr 
rückt zugleich immer näher an die Wand RE =. 1 nn 
heran. Da der in Gl. (2) auftretende Fak- 0,012 34.2 191 264 
tor (7/7), der bei der Erörterung des 0,014 47,8 194 219 
Wandgesetzes 1 gesetzt wurde, in Wirk- is RR 184 5 
lichkeit unter 1 liegt, so ist das Maximum 0.02 Er = 
von M* n,, insbesondere bei den. kleineren 0.03 92,6 100 89,3 
Re-Zahlen, auch kleiner als der aus dem 0,04 75,6 72,0 61,9 
Wandgesetz folgende Wert 0,25 n, (bei- ie 2 Be a. 
spielsweise beträgt das Maximum von M* 0.08 395 98 1 293 
n, bei n, = 400 nicht 100 sondern nur 96). 01 22.8 22,7 23.5 
Bild 3 läßt erkennen, daß die in die 0,2 10,7 11,2 11,4 
Nebenbewegung fließendeEnergie in Wand- n. en Mn Yen 
nähe ganz erheblich höher ist als in grö- 05 312 3.23 3.28 
ßBerer Entfernung von der Wand. In dem 0,6 2,03 2,10 Do 
mittleren Bereich des Strömungsquer- 0,7 1,17 1,21 1,24 
schnittes ist der Energiegewinn aus der us ne N I 
Hauptbewegung nur geringfügig. Wir 1.0 0 0% 0 


haben daher unsere Darstellung bei y/r —0,3 
abgebrochen (von wo ab M* n, allmählich auf Null abklingt). Das Mittengesetz der Energie- 
gewinnung dürfte daher im Vergleich zu dem zuvor behandelten Wandgesetz nur von 
geringer Bedeutung sein. 


6. Das Querschnittsmittel der Dissipation 

Wenn die Frage der örtlichen Dissipation nicht zur Debatte steht, sondern wenn lediglich 
die über dem Querschnitt gemittelte Dissipation interessiert, so kann man den durch Gl. (1) 
dargestellten Ausdruck der Energieproduktion M als ein Maß für die in der Volumen- und 
Zeiteinheit in Wärme umgewandelte Energie betrachten. Denn die durch die scheinbaren 
Spannungen aus der Hauptbewegung in die Nebenbewegung überführte Energie wird ja 
letzten Endes vollständig in Wärme verwandelt. Dabei ist es bei dieser Betrachtungsweise 
ohne Belang, wo und wie diese Umsetzung von Bewegungsenergie in Wärmeenergie erfolgt. 

Die durch die zähen Spannungen der Hauptbewegung entnommene Energie wird un- 
mittelbar in Wärme verwandelt. Nennt man diesen Dissipationsanteil je Volumen- und 
Zeiteinheit N, so gilt s 
I mc. (7). 

u 


Gesamtdissipation 
Im Rahmen dieser Mittelwertsbetrachtungen ist die Summe M + N ein Maß für die 
gesamte Dissipation der Haupt- und Nebenbewegung. Bezeichnet man M + N mit Z, so wird 
mit Rücksicht auf (1) und (7) L=r,„r/u. Macht man auch diesen Ausdruck durch Multipli- 
kation mit u/z3 — dimensionslos und setzt Z u/73 = L*, so wird, da r,,/7, = du/u*/dn ist, 
ee al naltg nein ba (8), 
To in % 


1 
Für das Querschnittsmittel | L* " das L}, heißen möge, folgt dann 


0 
Unax l Umasx 0 ] 
T du TU u AT 
ee en N E79), 


De 
ar Ar To U* To U* ur Ty 
Zr; 
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Hier verschwindet das Glied — es da an der Wand % = 0 und in der Kanalmitte 7 = 0 ist. 
0 


Wird noch die für die ausgebildete Strömung typische Beziehung 7/z, = 1 — y/r eingeführt, 
so folgt schließlich für die Gesamtdissipation der Kanalströmung 


1 
udYy Um 


* nn 
I Nr Fe au* Y U u* 
f 


(Gesamtverluste) . . . (10). 


(Hier bedeutet u, = [ “I die über dem Strömungsquerschnitt gemittelte Geschwindigkeit) 4 


Dieser ohne Einschränkung geltende Ausdruck für die mittlere Gesamtdissipation läßt 
sich auch aus dem Leistungsbedarf der Kanalströmung Ap F w„ entwickeln, der von einer 
Pumpe aufgebracht werden muß (hier ist Ap der Druckunterschied über der Kanallänge |, 
F der Strömungsquerschnitt und «, das Querschnittsmittel der Geschwindigkeit). Bedeutet 
L,„ den auf das Flüssigkeitsvolumen Fl bezogenen mittleren Leistungsbedarf , so folgt Z,, = 
AP Um/l. Da im flachrechteckigen Kanal Ap/l zry/r ist, wird Z, = 7,u„/r. Führt man 
noch I# = L„u/7 ein, so folgt I, = W. u/r ty Hieraus erhält man schließlich, wenn 
7%, = o u*? gesetzt wird, LA, n, = Un/u*. (Demgegenüber wird für die Rohrströmung LA, n, = 
2 u„/w* erhalten, da im Rohr Ap/l = 2 1y/r ist.) Of 


Dissipation der Hauptbewegung 


Nach Gl. (7) ist die unmittelbare Dissipation der Hauptbewegung je Volumen- und Zeit- 
einheit N = r„,/u. Wird diese Beziehung durch Multiplikation mit w/73 dimensionslos gemacht 
und N u/t? mit N* bezeichnet, so folgt, da 7/z, =1 + A/uist, 


le /% i 
(Sl (1a ee 


w* I'n’an 


RD 
(rn)? 


N VE 


Bild4 


Wir bilden das Querschnittsmittel N#, = j N* z und integrieren abschnittsweise und 


zwar zunächst über die Wandzone (O<n< Ar wo die Wandgesetze auf der Basis = 
ae und dann über den mittleren Bereich (ya <Y <r), wo die Mittengesetze gelten, und 
erhalten: 


* 7/t, \?d 
an ar any® ae re 


I 


?) Eine entsprechende Ableitung für die Rohrströmung ergibt: ee Nr = 2 um|u*. 


EEE 
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Der Integrand des ersten Integrales klingt stark mit n ab (siehe Bild 4). Dieses Integral 
wächst somit oberhalb eines gewissen n-Wertes nur noch schwach mit n®). Beispielsweise 
liefert die graphische Integration für I den Wert 8,830, wenn sie bis n, — 36 durchgeführt 
wird. Integriert man weiter bis 7. = 50, so wird I — 8,866. Das ist also nur eine geringfügige 
Erhöhung. 


Bei der Berechnung des Integrales II ist zu beachten, daß 7/ =1 — y/r und A/un, 


näherungsweise n (1 & gesetzt werden kann (die letztere Aussage beinhaltet, daß die 
logarithmische Geschwindigkeitsverteilung bis zur Kanalmitte reichend angenommen wird). 
Aus dem zweiten Integral von (12) folgt dann: 


1 1 
Be Falke SE (13) 


Da x = 0,4 ist, erhält man für n, = 36 II = 0,173 und für n, = 50 II = 0,125. Bilden 
wir nun die Summe I + II, so wird für , =36 I + II = 9,003 und für, = I+II = 
8,991. I + Il ist also praktisch unabhängig von dem gewählten n, ziemlich genau 9. Somit 
folgt aus (12) für die Dissipation der Hauptbewegung in der ebenen Kanalströmung bei hin- 
reichend hohen Re-Zahlen die Näherungsformel 


Nmn, =9 (Verluste der Hauptbewegung) . . . . (14) 
(Demgegenüber wird für die Rohrströmung N?, n, = 18 erhalten.) 


Das Verhältnis II/I liegt je nach Festlegung der etwas willkürlichen Grenze n, zwischen 
0,01 und 0,02. Bei der Hauptbewegung ist also die Dissipation im ausgedehnten mittleren 
Bereich des Strömungsquerschnittes nur ein kleiner Bruchteil der Dissipation in den dünnen 
Wandschichten. Mit anderen Worten, der Energieverbrauch der Hauptbewegung findet 
praktisch in unmittelbarer Wandnähe statt. 


Dissipationder Nebenbewegung 


Das Querschnittsmittel der Dissipation der Nebenbewegung M7, ist gleich dem Quer- 
schnittsmittel der Gesamtdissipation Z7, vermindert um das Querschnittsmittel der Dissipation 
der Hauptbewegung N}. Wir erhalten daher für die Dissipation der Nebenbewegung der 
ebenen Kanalströmung bei hinreichend hohen Re-Zahlen aus (10) und (14) die Näherungs- 
formel 
Um 
u* 


Mann =— —)I (Turbulenzenergie) ..... » (15). 


Die entsprechende Formel für die Rohrströmung lautet: M7,n,=2 m 9) \ 


Formel (15) ist zugleich ein Ausdruck für das Querschnittsmittel der je Volumen- und 
Zeiteinheit aus der Hauptbewegungge wonnenen Energie der Nebenbewegung. Da My, n,; 
dimensionslos ist (durch Bezugnahme auf r?/u), so wird der Betrag an Turbulenzenergie im 
Flüssigkeitsvolumen Fl der Kanalströmung dadurch erhalten, daß man die rechte Seite von 
(15) mit FI 75/u multipliziert. 

Das durch Gl. (15) beschriebene Verhalten von My, n, ist leicht zu überblicken. Da u/u* 
mit n ansteigt, wachsen %,.,/u* und u„/u* mit n, bzw. mit der Reynoldsschen Zahl. 
Demnach wächst auch M#,n, mit der Re-Zahl. Verwendet man für w„/u* = Re/2n, die von 
J. Nikuradse?°) gemessenen Werte, so erhält man die in Bild 5 über log Re aufgetragenen 
Mn Werte. L 

Diese Kurve beginnt an der unteren Grenze der Reynoldsschen Zahl Re, = 3000 
mit dem Werte M#,n, =5. Allerdings sind die My, n,-Werte in dem unteren Re-Zahlbereich 
ziemlich ungenau, da die benutzte Gl. (15) hohe Re-Zahlen voraussetzt 9). 


8) Es ist daher innerhalb gewisser Grenzen ziemlich unwesentlich, welcher n,-Wert als „Grenze“ der 
Wandschicht festgelegt wird. Wesentlich ist jedoch, daß die Re-Zahl bzw. das Verhältnis n,/na hinreichend 
großist, damit die für den Bereich 0 < n< n, vorausgesetzte Bedingung 7 = 7, auch erfüllt wird. 3 

®) J. Nikuradse, Gesetzmäßigkeiten der turbulenten Strömung in glatten Rohren. VDI-For- 
schungsheft 356, 1932. \ , \ j 

10) Das in Gl. (15) zum Abzug gebrachte N*,», ist bei den kleineren Re-Zahlen kleiner als 9, da das in 
die N*-Formel (11) eingehende (r/z,)? (das wir für hohe Re-Zahlen 1 setzen durften) merklich kleiner ist als 1. 
Demnach sind die wirklichen M* n,-Werte in dem unteren Re-Zahlbereich größer als die hier angegebenen 
Werte. 
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%. Erlöschen der Turbulenz 


Wir wollen zum Abschluß die schon oft behandelte Frage erörtern, weshalb der turbu- 
lente Strömungszustand unterhalb einer bestimmten Reynoldsschen Zahl nicht mehr 
existenzfähig ist. Zu dem Zweck betrachten wir die Energiegleichungen in der durch Gl. (3a) 
dargestellten Form. Wir denken uns diese Gleichungen addiert und über y/r integriert. Da die 
Integrale der Diffusionsglieder fortfallen, erscheint dann auf der linken Seite der Ausdruck 


1 
| M* = 7, = M5,n,, den wir im vorigen Kapitel durch die Formel (15) dargestellt haben. Die 
0 


Gleichung für das Querschnittsmittel der Turbulenzenergie der ausgebildeten Kanalströmung 
lautet daher: 


2 \2 \ 2 2 \dı 
I CAESC EIER EEnEezt) 
0 


Daß die Turbulenz unterhalb einer gewissen Reynoldsschen Zahl plötzlich ver- 
schwindet, läßt sich aus dem Verhalten von M7, n, nicht erklären. Diese Größe ist bei der kriti- 
schen Re-Zahl keineswegs Null, sondern hat min- 
destens den Wert 5. Das Erlöschen der Turbulenz 
kann daher nur auf das Verhalten der Glieder der 
rechten Seite der Gl. (16) zurück zuführen sein. 

Der auf der rechten Seite von (16) stehende 
Faktor 1/n, wächst mit abnehmendem n, bzw. mit 
abnehmendem Re. Hingegen sind die dimensions- 
losen Schwankungsgeschwindigkeiten uw = w/u* usw. 
nach unseren experimentellen Feststellungen von der 
Re-Zahl ziemlich unabhängig. Eine Verkleinerung 
von Mn, mit abnehmender Re-Zahl, wie sie durch 
Formel (15) vorgeschrieben wird (siehe Bild 5), kann 
daher nur von einer Vergrößerung der Abmessungen 
Yn, /r, Yn, y/r und Yn, 2/r der Turbulenzelemente mit 
abnehmendem n, herrühren. 

Genauere Aussagen über die komplizierte nicht- 
6 isotrope Wandturbulenz sind zur Zeit nicht möglich. 

Setzt man, um überhaupt zu einer Formel zu ge- 

langen, isotrope Turbulenz voraus, so läßt sich der 
bekannte, von G. I. Taylor entwickelte Dissipationsausdruck verwenden und Gl. (16) 
geht über in 


4 
log Re —— 
Bild 5 


5 uw 
n. (A/r)? 
Hier ist A ein ungefähres Maß für den Durchmesser der kleinsten Wirbel, deren Bewegungs- 


energie in Wärmeenergie umgesetzt wird. 
Aus (15) und (17) folgt dann weiter 


11/2 
u] ja; u SM, one Wise Io runrEs 


r/A kann gedeutet werden als die Anzahl der kleinsten Wirbel, die auf einem Halbmesser Platz 
haben. 

Nach unseren Feststellungen ist das Querschnittsmittel von u? etwa 1,25. Mit der Ver- 
wendung dieses Wertes ergeben sich dann aus (18) folgende »/A-Werte: 


Re | 3.103 10: 105 106 
ia | 5,3 11 3947 431 


Wenn auch diese Angaben für die nichtisotrope Wandturbulenz keineswegs verbindlich 
sind, so geht doch aus der vorstehenden Tabelle deutlich die Tendenz der Gl. (16) hervor, daß 
mit abnehmender Re-Zahl die Ausdehnung der Wirbel wächst, bzw. die Anzahl der Wirbel 
geringer wird. Selbstverständlich können die kleinsten Wirbel nicht größer sein als die Stärke 


M%n,= a a7 
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der Reibungsschicht (es ist also allein aus Platzgründen r/A> 1). Auch das Auftreten einzelner 
Wirbel in der Reibungsschicht bedeutet noch keine Turbulenz, sondern der turbulente Strö- 
mungszustand ist gerade durch das Auftreten zahlreicher größerer und kleinerer Wirbelgebilde 
gekennzeichnet. 

Dieser Strömungszustand verändert sich mit abnehmender Re-Zahl in der Weise, daß 
vornehmlich die kleinen energieverbrauchenden Wirbel verschwinden. Schließlich bleiben nur 
noch wenige grobe Wirbelballen übrig, die nicht mehr in der Lage sind, einen turbulenten Aus- 
tausch zu bewerkstelligen und der Hauptbewegung soviel Energie zu entziehen, wie zur 
Aufrechterhaltung des turbulenten Zustandes erforderlich wäre. 


Eingegangen am 26. Juli 1952. 


The physical components of vectors and tensors 
By C. Truesdellin Bloomington, Indiana, USA }) 


Es werden physikalische Komponenten von Vektoren und Tensoren, bezogen auf allgemeine krumm- 
linige Koordinaten, definiert, und es wird gezeigt, daß diese Komponenten Größen darstellen, die die natür- 
lichen physikalischen Dimensionen des Feldes besitzen und einer direkten physikalischen Interpretation 
fähig sind. Weiter wird gezeigt, daß skalare Funktionen von Tensoren zweiter Ordnung sich direkt durch 
die Invarianten der physikalischen Komponenten ausdrücken lassen. Für den Fall orthogonaler Koor- 
dinaten wird eine Formel für die physikalischen Komponenten der kovarianten Ableitung eines Tensors 
beliebiger Ordnung abgeleitet. 

Physical components of vectors and tensors referred to general curvilinear co-ordinates are defined and 
shown to represent quantities possessed of the natural physical dimensions of the field and capable of imme- 
diate physical interpretation. It is shown that scalar functions of second order tensors may be expressed 
directly in terms of the invariants of Ihe physical components. For the case of orthogonal co-ordinates, & 
general formula for the physial components of the covariant derivative of @ tensor of any order is derived. 


Des composantes physicales de vecteurs et tenseurs refertes ü des coordonnees generales curvilignes 
sont definies, et il est montr& que ces composantes reprösentent des quantit&s qui possedent les dimensions 
naturelles physicales du champs et peuvent ätre soumises & une interpr &tation physicale directe. De plus, il 
est montr& que des fonctions scalaires de tenseurs de second ordre peuvent Etre exprimees directement par les 
invariantes des composantes physicales. Pour le cas de coordonnees orthogonales, une formule pour les com- 
posantes physicales de la derivation co-variante d’un tenseur de n’importe quel ordre est derivee. 


B eraTbe ompepensirtca $HCcHYecKHe KOMIOHEHTEI BEKTOPOB H TEH30POB, OTHECEHHEIE K 
OÖIIHM KPHBOJIUHHeÜHBIM KOOPAHHATAM; HOKA3EIBAeTCH, YTO 3TH KOMHOHEHTEI IIPeACTABIAOT 
co6oU BEIIHYHHEI, OONAAAIOIMME ECTECTBEHHEIMH PHSMYeCKUMU PA3MEPHOCTAMH MONA H NOILYCK&- 
IOINMHe HeNOCpeACcTBeHHyIW dPH3nNYecKyW® HHTEePpNpeTaumm. Jlaree ImOoRasaH0, YTO CKANAPHEIE 
ÖYHEIHH TeH30POB BTOPOTO HOPAAKA MOTYT ÖBITb HeIMOCcPe]CTBeHHO BEIPA3KeHBI Tepe HHBAPH- 
AHTbBI DUSMHYeCKHX KOMIIOHeHT. Ha cayyali OPTOTOHANABHEIX KOOPAHHAT BEIBONHTCA POopMyıla 
AA BH3UYecKHX KOMIIOHEHT KOBAPHAHTHOH IIPOM3SBOAHOH TEeH30P& IIPOH3BOJIBHOTO HOPAMKA,. 


1. Nature and history of the subjeet 


In mathematical physics and in elementary Euclidean geometry a vector Zisregarded as a 
directed line segment. Its components in a given co-ordinate system (z!,..., =”) are also line 
segments; the &! component, for example, which we may call A (?), is a line segment tangent to 
the z! co-ordinate curve at the initial point of the vector and extending thence to the plane 
which passes through the terminal point of the vector and is paralleltothe plane of the remai- 
ning co-ordinate directions #°,..., x” at the initial point (Fig. 1). 

Such components add by the parallellogram rule to form the vector itself. It was these 
quantities which Ricci & Levi-Civita in their great memoir on tensor analysis (1901, see Ch. I, $4) 
called ‘‘composantes du vecteur selon les lignes coordonn&es.”’ In books designed for mathema- 
ticians these components are no longer mentioned, however; what Ricci and Levi-Civita called 
“syst&mes covariants et contrevariants’’ are now called the covariant and contravariant com- 
ponents A, and A of the field, distinguished respectively by the transformation laws 


- 0: RECLZEN 
ey, ne ER E70 4} 
under change of co-ordinates Ü=f!(#,..., ="). 


The contravariant and covariant components differ in general from the components 
used in mathematical physics; indeed they are sometimes difficult to interpret physically and 
they often fail to bear the physical dimension of the vector field to which they belong. Consider, 


1) Professor of Mathematics, Graduate Institute for Applied Mathematics, Indiana University. This 
paper was written for the Applied Mathematics Branch, U. S. Naval Research Laboratory, Washington, D. C. 
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for example, the gradient of a scalar f(x!, 22) in plane polar co-ordinates !=r, ?—=©. 
For the covariant components f;, we have 


OPEN 
m RL EN. Se 


If f be a velocity-potential, say, then Af/Or is of the proper dimension (length)/(time) for a 

velocity, but we obtain an azimuthal “velocity” 0%/00 whose dimensions are (length)?/ (time), 

% which is certainly not a velocity 

2 Aa in the physical sense. In books 

on physics or elementary vec- 

tor analysis one finds the for- 
mula 


r e ef 
srad f = tg (113): 


instead of 0/80 we have now 
(Y/r) 8/80, which bears the 
same physical dimensions as 
x!curve Of/Or, and the quantities 'Of/Or 
and (1/r) 07/80 are indeed the 
components in the geometrical 
sense illustrated in Fig. 1. The 
contravariant components 0f/or, 
(1/r)? 0f/00, are different again 
from either of thetwo previous. 


In their treatment of ten- 
sors of order higher than 1, de- 
fined bythe now familiar laws of 
transformation analogous to (1.1), Ricci and Levi-Civita did not attempt to determine quanti- 
ties susceptible of direct physical interpretration. In their applications to physics (1901, Ch. VI) 
they were content to derive forms of various classical equations valid in any co-ordinate system, 
not mentioning that those forms are not the same as those usually employed in mathematical 
physies. While in each case their end results depart from the usual forms only by an insignifi- 
cant factor, the intermediate steps are often of a wholly different type. E. g., their formula for 
the stress-strain relations in a non-isotropic elastic solid if written out in a curvilinear co- 
ordinate system would look quite other than that usually employed in elasticity books, since it 
gives contravariant stress components in terms of covariant strain components and material 
constants forming a contravariant tensor of fourth order. It is not immediately obvious how 
any of the components of these three tensors could be expressed in terms of measurable 
quantities. 

While the problem of applying tensor formulae to classical physics was taken up by 
Murnaghan (1928), he did not try to relate tensor component to the components used in 
elementary vector analysis, being content simply to find quantities “physically useful”. For 
vectors he discussed the two types of orthogonal projections which had been calculated long 
before by Ricei & Levi-Civita, but he did not mention what they had called “composantes du 
vecteur sur leslignes coordonne&es,’ which are really the “components” in the elementary sense. 
His treatment of tensors of higher order is difficult to follow; his development of elasticity, for 
example, besides being elaborate, confuses matters by employing a tensor of third order to 
represent stress instead of the usual and quite adequate second order tensor. 

A satisfactory explanation of the relation between tensor analysis and Lame’s theory of 
orthogonal curvilinear co-ordinates was first given by McConnell (1931, Appendix, $ 2). 
His method might be explained as follows. Given a tensor field referred to an orthogonal 
curvilinear co-ordinate system, let its physical components at a given point be defined as its 
components in a rectangular Cartesian system whose axes are parallel to the curvilinear co- 
ordinate curves at the point. These physical components are exactly the quantities which are 
called simply components in works on vector analysis or physiss. McConnell’s definition 
is perfectly clear and perfectly adequate, so long as the co-ordinate system be orthogonal. In 
some cases, however, it is desirable to have quantities susceptible of direct physical measure- 
ment even in an oblique co-ordinate system. 

Synge&Schild(1949, 85.1), Ollendorff(1950, III5),andGreen&Zerna 
(1950, equation (5.16)) define physical components of a tensor of the second order. Their 


Fig 1 
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definitions differ one from another except in the case of an orthogonal system, when they re- 
duce tt McConnell’; I do not follow the remarks motivating these definitions and I do 
not think that any of them represent what physicists mean by tensor components in an 
oblique curvilinear system. 

In 1944 I first had the puzzling experience of getting an apparently wrong form of the 
equations of elasticity in polar coordinates by specializing the perfectly correct general tensor 
equations to the case at hand. A few hours of thought were required to find and remove the 
cause of the trouble. Since that time I have had to explain the matter many times, both to 
colleagues and to students. Since the definition and method I employ differs from those cited 
above, except in the orthogonal case, and since it seems perfectly direct and natural, it may be 
worthwhile to publish it. 

The reader of this paper is assumed to be familiar with the elements of tensor analysis, as 
presented e.g. in Eisenhart (1940). 


2. Definition of Physical Components of Absolute Vector Fields ?) 

In any co-ordinate system «in a Riemannian space of any number of dimensions 
denotethe covariant components of the metric tensor by g;,, so that the squared element of arc 
length is given by 

ds? = 9, dal dad BER IE rn. 2. ns (21), 


where, as henceforth, we employ the summation convention for diagonally repeated indices. 
The connection with notations used in the theory of orthogonal curvilinear co-ordinates is as 
follows: European writers usually put A,”?= g;,, while American writers put A? =g;;,. In 
orthogonal co-ordinate systems we have also the simplifying relations 


 luirnanp, gi SACHE ARTEN 
The treatment in the present paper is not confined to the orthogonal case, but only real co- 
cordinate systems in aspace where (2.1) is positive definite are considered. 

Let A? be the contravariant components of an absolute vector field in the co-ordinate 
system x. Then the squared length of the field is the scalar g,,4°4. A well known theorem of 
tensor analysis (e. g. (1940), [16.6]) tells us that this length equals the length of the diagonal 
of a parallelepiped whose sides are parallel to the co-ordinate curves and whose edges are of 
length yg; A (unsummed). InaEuclidean space the result is perfectly general, but in a 
non-Euclidean Riemannian space one must restricet attention to infinitesimal 
fields A, since only infinitesimal parallelepipeds exist. Comparing the foregoing theorem with 
the opening paragraph of $ 1 shows that analytical definition of the physical components A (i) 
must be 

Moe Von A (unaummed) Pen nd rn. . (28). 


The physical components, of course, do not obey the simple laws of transformation (1.1) 
which are enjoyed by the contravariant or covariant components. From (2.3) and (1.1), we 
have the more complicated transformation law 


in - VE (unse ae RE UNE LEN, 


whose form is not independent of the co-ordinate system. It is possible to interpret the coeffi- 
cients in (2.4) in terms of the angles between the co-ordinate curves inthe zand # systems. 

Suppose now we wish an expression for the physical components A(i) in terms of the 
covariant components A, . Since Ai = gÜj A,, we have from (2.3) 


A)—=Yg:g7A (unsummed) ......2.... (25). 
Thus the physical components are linear combinations of the covariant components. In the 


special case of an orthogonal co-ordinate system, however, by (2.2) we may put (2.5) into the 
simpler form 
PIE en er ire2.6). 
Ye; 

2) The formula (2.3) and its interpretation were given by Ricci & Levi-Civita (1901, Ch. I, $ 4). They 
caleulated in fact four sets of lengths connected with a vector field: “Nous allons determiner pour des coor- 
donnees gön6rales x,, 25, 2; les expressions des projections R,, et En, et des composantes R;, et R„, selon les 
tangentes aux lignes et selon les normales aux surfaces coordonnges”’. These lengths are given by AlVgii, 


ZilVgii, Vai: At, Vgii A; (unsummed). 
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The formulae (2.3) and (2.6) for the orthogonal case are given by McConnelland 
more recent writers. They yield, for example, an explanation of the difference between (1.2) 
and (1.3), since in plane polar co-ordinates we have 9; = 1, 912 = 9a =, 92 = r?, so that 
from (1.2) and either (2.3) or (2.6) we obtain (1.3). 


3. Dot and cross produets 
Let ® and »i denote, respectively, the dot (scalar) and cross (vector) products of the two 
fields A‘, ui; the latter can be formed only in a three-dimensional space. For the tensor com- 
ponents we have 
Behr aeiyagkimg’hu, 
ijk en 
ee v, = Yg Ein it uf 
where eÜ/k — g,,,.isthe usual permutation symbol and g = det g;,. These formulae show that 
the simple rules for caleulating dot and cross products in Cartesian co-ordinates do not carry 
over to the curvilinear case, so long as covariant or contravariant components are used. 
If we use the physical components (2.3) the situation at first sight appears to be worse: 


® DE v(i) = Han FH aa) au u re 

9:95; Irr Irı 
When the co-ordinate system is orthogonal, however, by (2.2) these awkward formulae reduce to 
D=Ali)uti), vi) as All 2 en a 


These equations embody 


Theorem 1: In an orthogonal curvilinear co-ordinate system the dot and cross products may be 
calculated from the physical components by the usual formulae of elementary vector analysis. 


4. Physical dimensions of the physical components of vectors 

From (2.3) and (2.6) it is plain that in a rectangular Cartesian system the distinction bet- 
ween covariant, contravariant, and physical components disappears, and we may speak of the 
components. Now why are physical components used at all? The only valid reason is that they 
carry the proper physical dimensions, i.e. the same dimensions as the rectangular Cartesian 
components, while the contravariant or covariant components generally do not. We now go 
into this point in detail. 

Rectangular Cartesian co-ordinates are a special case of the Riemannian normal 
co-ordinates (see e.g. Veblen (1927), Ch. VI, $ 16) available locallyinany Riemannian 
space. Let such co-ordinates be denoted by a bar: &,. Since the element of arc isthen 


dei (FE un hf Br Tea a Ar 
i=1 


„all the normal co-ordinates & have the dimension of length, and, consonant with usage in the 
Euclidean case, we shall suppose that for any field A? having physical significance all the 
normal components A? be assigned a common physical dimension, which we shall call the natural 
dimension of the field. From (1.1), we have then 


ox8_ dim «° -. 
dim 2 = dim (02,7) — im. uf era Aue 
ri 92) 7 Tengeh: and er 
Thus except in the special case when the co-ordinates @° possess the dimension of length the 


contravariant components A! will possess physical dimensions different from dim Ai and diffe- 
rent one from another (e. g., in plane polars dim «! = dim r =length, dim =? = dim 9 = 0). 
But from (2.1) we have 


dim Yg;, dim a length 2.27 308% v a BE ER Et 5 
in any co-ordinate system «‘. By taking dimensions in (2.3) and by using both (4.2) and (4.3) 
we obtain then 

length dim x 


dim A(i) = dim yg;; dim Ai — ) 
Feen im Y9s: dimA dim «* length 


dimi—=dim#..... (4.4). 
Thus we have. 


Theorem 2: Ar the physical components A(i) have a common physical dimension, and that 
dimension is the natural physical dimension dim At of the field. 
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5. Definition of physical components of tensors of the second order 

To find a natural and proper definition of physical components of a tensor of second order 
one must understand how such tensors arise in physics. Now indeed, how do tensors occur ? A 
little reflection reveals that in classical physics, only scalars and vectors have primary physical 
significance. That is, whenever a tensor occurs, it is defined in terms of scalars and vectors. 
Now a vector can be obtained from a tensor of the second order in one of three ways: by con- 
tracting the outer product of the tensor with a vector, or by taking the divergence of the tensor, 
or by a more elaborate operation derived from repetition and combination of these two basic 
operations. A similar remark applies to scalars. 

Both operations can be illustrated by the stress tensor in continuum mechanics. Accor- 
ding to one method, it may be introduced as a tensor t“, such that if n’ be the unit normal to a 
surface, then 


ME... ld) 
is the force per unit area acting across this surface. We could write instead 
re ER agents ae ld), 


but then 2; would naturally be a different tensor from t',;, except in the special case when either 
is symmetric. In both cases, fi and n’ are absolute contravariant vector fields. The second 
method of introducing stress would start from Cauch y’s statical equation. 


and state that 1“ is that tensor whose divergence equals the resultant per unit volume of all 
surface forcesinthe medium. Alternatively we could write f? ‚ for this effective force, and then 
again the stress tensor would have different meanings in the two cases, unless it happened to be 
symmetric. 

Since in all cases the interpretation of tensors can be reduced to the interpretation of 
vectors and scalars, and since physical components for the former have already been given, we 
could dispense altogether with physical components for tensors. To do so, however, would be 
rather awkward sometimes. In books on dyadics the relations (5.1) and (5.2), respectively, are 
written 


ee KR en 1 I A . (9.4). 


What is meant by T in the two cases ? Of course it is an ordered set of n? quantities, or, if we 
prefer, like 13, ti,,t/, and t“/anditisann X n matrix, but what are the entries in that matrix? 
The answer to this question is plain from the foregoing example. Let a‘; be a given abso- 
lute mixed field, let 4 be a given absolute contravariant field; then absolute contravariant 
fields „a“ and vi are defined by 
wW=a,N, vw=Naf. A a N 
Physical components for the fields A, u, v are defined by (2.3). By right physical components 
a(ij) of the field ai, we shall mean a set of quantities such that 


le PATE ne an nee. 86); 


that is, right physical components a(ij) of a’, are a set of quantities such that scalar multiplication 
on the right by physical components A(j) of a field A} yields the physical components u(i) of the field u 
defined by uw‘ = a‘, )). Left physical components (ij)@are analogously defined. 

Thus two sets of physical components correspond to each mixed tensor a‘;. Which is to 
be employed depends on the way the tensor is introduced. In the case of the tress tensor, for 
example, if (5.1) be the defining equation then right physical components £(ij) will be what the 
physicist means by stresses; while if (5.2) be employed instead, then left physical components 
(ij)t are required. The difference is a matter only of notation, but consistencey must be main- 
tained. If there were general agreement about the preferred form of tensor equations in phy- 
sics, only one type of physical component would be needed. 

To obtain a formula for a(ij), we begin by expressing both u‘ and A in (5.5) in terms of the 
physical components w(i) and A(j) by (2.3). Then 


u) _ ,; AU) 


— =al, — (funsummed) . ...... 60.7) 
Vi V9;; 


or 
un = 420) (i unsummed).. . . - . - (5.8). 
ii 
25 
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Hence the only possible set of quantities @(ij) which can satisfy (5.6) is given by 


alij)= a, (unsummed) «i kudleben un@9): 


Similar reasoning can be applied to find the left physical components. The foregoing analysis 
is summarized in ’ 
Theorem 3: Let a‘; be a given mixed absolute tensor field, and let A be an arbitrary contra- 
variant vector field. Physical components of vectors being defined by (2.3), in order that the physical 
components of al, 2 (A ajt) be atij) A(j) (Ali) (i)a) in any co-ordinate system, it is necessary and 
suffievent that 


a (ina= a) (unsummed) . . . . ..(5.10). 
95; I; 


39 
The result (5.10) shows that 
ai) =.at (i)amia)- Kunsummed) Zn. 
That is, diagonal mixed tensor components are right or left physical components, as the case 


may be, but off-diagonal mixed tensor components in general are not physical components. 
By summing (5.11) on ? we obtain 


Ki) Bahia ee Deere 


aresult we may express as 
Theorem 4: The four matrices a',, a, (ij) a, a(ij) have the same trace. 
The definitions of physical components (5.10) are easily expressed in matrix form. Let 
A, stand for the matrix ||a,||, „4 for the matrix ||a,||, @ for ||g;;||, where in each case the left 
hand index is the row index and the right hand index the column index. Then ‚A =@4,@". 
Write 
RB =] 70, (unsummedy 4. 22. ie SUMME en 
where ö,;isthe Kronecker symbol. Now putting A, for the matrix ||«(ij)|| of right 
physical components, „A for the matrix ||(ij) a|| of left physical components, we have 
A, = RAR LAN RAW Lg Re 


Theorem 4 could also be proved by using (5.13) and an easy theorem on matrices. 
The line of further development suggested by Theorem 4 is taken up again in $ 10. 


6. Relations between right and left eomponents. Symmetrie tensors 


From (5.10), and (5.10), we have 


na Hg, g’!a*, (i,7 unsummed), | 
AR EEE EN 


For Giga 
I; Ir 
This formula connects the left and right physical components of a given field «‘,. 
By (2.2) we conclude from inspection of (6.1) 
Theorem 5: In an orthogonal co-ordinate system the left physical components (ij)a of any tensor 
field a’; equal the right components a(ij), both in turn being given by MeConnell's formulae 


9 ’akkı) (%, 7 unsummed) | 


(ia alij) = 5. == V& a’ Yg ga! = BR (unsummed) . . (6.2). 
I; ii 19 

Returning to a general co-ordinate system, consider the special case when the tensor «a;,is 
symmetric: a,; = a,;. It is-plain that in general the matrix of right (left) physical components 
is not asymmetrie matrix: a(ij)#a(ji) ((ij) a (ji)a). We have, in fact, for a symmetric field 


Fri /g:: 
na Ha Bor-aun . ie 


This result expresses 


Theorem 6: If the tensor a,, be symmetric, then the matrix of its left physical components is the 
transpose of the matrix of its right physical components. 
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In the case of a symmetric stress tensor, for example, as we mentioned earlier it would 
make no difference whether (5.1) or (5.2) were used to introduce it. This conclusion follows also 
from Theorem 6, since we have then, in physical components, 


IN HMMM. 2:8 22. rn (ÖA). 
Combining Theorem 5 with Theorem 6 yields 
Theorem 7: In an orthogonal co-ordinate system the matrix of left physical components of a 
Symmetrie tensor is equal to the matrix of right physical components and is a symmetric matrix. 
The strain tensors, for example, are symmetric by definition, and the stress tensor is 
symmetric if there be no extraneous couples. These tensors are symmetric in the usual tensor 
sense: t;;, —1;;, etc. Hence their right and left physical components generally do not form 
symmetric matrices in an oblique co-ordinate system. 


7. Tensor divergences 

We return to the second way in which tensors are introduced into classical physics, viz., 
through their divergences. Since covariant differentiation may be regarded formally as multi- 
plication by a symbolic vector p,, the rules of $5 may be applied in this case also. If it is 
@‘,;—=y;@‘ which has physical meaning, then left physical components (ji)@ should be 
introduced; if it is «7, which occurs, then right physical components a(ij) should be used. 

For the special case of an orthogonal system, when the difference between right and left 
physical components disappears, covariant differentiation will be discussed from a more general 
viewpoint in $ 12. 


8. Physical dimensions of the physical components of second order tensors 


Asin $4, let bars denote aRiemannian normal co-ordinate system, so that dim a; is 
the proper physical dimension of the field «‘,. Then since the tensor transformation law gives 


oaloT _. xt 00 a 
= Er Hyatt ee em tar (8.1), 
we have 
KLLIT GB a U Ge De un a re (0,2), 
but 
: dim our m 
dim al, — dm „; im a‘; NER SDE EN URNENT EN 0 U, WB, (8.3). 


Thus the diagonal mixed components a!,, ..., @*„ have the proper physical dimensions, but 
the off-diagonal components do not. But for the right physical components «(ij) we have from 
(5.10), (8.2) and (4.3) 

dim g;; 
dim g;; 


dim a(ij) = dim «@‘; 


dim # dim& . ’ u 
ne leer dnkay =dımay. „2. nua v (Bd), 


so that all the right physical components enjoy the natural physical dimensions of the field. 
Since a similar analysis can be carried through for the left components, we have 

Theorem 8: For any absolute field a‘, in any co-ordinate system, all the right and the left physical 
components have a common physical dimension, and that dimension is the natural physical dimen- 
sion of the field. 

From the foregoing theorems it follows, for example, that to obtain from the general 
tensor equations the formulae concerning stress and straincomponents in orthogonal curvilinear 
co-ordinates given in books on elasticity, for the diagonal components one may simply take the 
mixed tensor components, but for the off-diagonal components one must multiply the mixed 


tensor components by suitable factors Yg;,/9;;- 


9. Proper values and prineipal axes 
The proper values and principal axes of second order tensors often have physical signifi- 
cance. We now establish the connection between these quantities as calculated from the four 
matrices a‘, a, (ij)a, a(ij). 
Now wis a proper value and Alis a left proper vector of a if 
ee... (91); 


Im 


g43) 
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u;is a right proper vector if 
a,’ Mi; EIN ki Ad . . . . . . . . . . . . ® . . (9.2). 
In either case x must satisfy 
det (a? —& INNEN uichanke: Tee Ferse (9.3); 
The direction of Alis that of a left principal axis corresponding to the proper value a; that of u; 
is that of aright principal axis. In the case of a symmetriec tensor a/, the distinetion between 
right and left proper vectors and principal axes disappears. 
By raising and lowering indices in (9.1) and (9.2) we get 
Kor ai, u, WON“. ee a ee De Pte EiL: (9.4). 
Since Aland A; (u; and u‘) represent the same vector, by comparing (9.1), (9.2) and (9.4) we get 
Lemma 1. The matrices a} and ai; have common proper values and common left (right) prin- 


cipal axes. 
If we write (9.1) in terms of the physical components (2.3) and (5.10), we have 


ae Veina=a a a eu 
Y: 9 V9; 
whence 
AGD eÜaeee 5 Ir In ae Te: 


Comparing (9.1) and (9.6), and carrying out a similar manipulation starting with (9.2), we obtain 
Lemma 2. The matrices a} and (ij)a (a? and a(ij)) have common proper values and common left 
(right) principal azes. 

Combining lemmas 1 and 2 yields 
Theorem 9: The four matries aj, a';, (ij)a, alij) have common proper values. The left (right) 
principal awes of a} coincide with those of the matrix of its left (right) physical components (ij)@ (a(ij)). 
If a} be symmetric, then the four matrices a}, a’; , (ij)@, a(ij) have common principal axes. 


10. Invariants 


Braun (1951) has noted that in many cases physical components may be used asifthey 
were the mixed tensorial components. Precise statements expressing this idea will now be given. 
First we have the rather trivial 


Theorem 10: The matrix of right (left) physical components of the tensor whose matrix of mixed 
components is the inner product of the matrices of the mized components a‘; and bi, (a; and b/) 
(a', and b,) vs itself the product of the corresponding matrices of right (left) (right and left. respec- 
tively) physical components a(ij) ) and b(ij) ((ij)a and (ij)b) (a(ij) and (ij)b). 

The proof is immediate from the definitions (5.10): 


y& ai; bi, - as 2077 Zend > N ON 
ya Au Weka y@ or; all) IHR. ee 102 


Theorem 10 is also immediately evident from the matrix formulation (5.13). 

From the theorem we have 

Corollary 1. The inverse of the matriz of right (left) physical ee of ai;(af), if it 
ewist, is the matrix of right (left) physical components of the inverse of ai, (a); and 

Corollary 2. The matrix of the right (left) physical components of any power, positive or 
negative, of the matrix of mixed components of a tensor a‘; (a;}) is the same power of the matrix of 
right (left) physical components a(ij) ((ij) «). 

We have also the following 
Theorem 11: Any scalar function of the mixed components a‘; and a;' of certain temsor nz may 
be regarded alternatively as the same function of the physical components a(ij) and (ji)a 


For this theorem and for the proof which follows now I am indebted to my colleague 
G.Whaples;Ihad originally derived a weaker result in a more elaborate way. 
Lemma. Let subseript naughts indicate evaluations at a particular point. Then it is pos- 
sible to find a co-ordinate system % such that at any single arbitrarily selected point we have 


(2), = (atij))., W(afseliha) - a 


Also 
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For if we take 


= (195) 5 unsummed) ah. und ler autom \\ (10.4), 
then, since (Yg;;)o is a constant, we have 
om es — dx %= Ö,; r 
ee allge; Er (a (unsunmmed)‘ .*.. . ....(10.5), 


so that by the tensor law of transformation 


ea Oki dam 9; \ 
al, = na = | Ku) 0; Kunsunmed)r 20. 2er 2.2 (10.6) 
Evaluating the result at the point in question, by (5.10), we derive (10.3),. A similar analysis 
yields (10.3),. 

The truth of theorem 11 now becomes evident, since a scalar function f(a',, aj') has the 
same value in all coordinate systems, and in particular we shall have 


Fa)» (a’)) =Fl@), (W) =FlatidI, (GdA)ı): 
but since this result holds at any arbitrarily selected point, it holds at all points, and the 
naughts may be dropped. 


11. Tensors of order greater than 2 
While it is not difficult to write out in full the definition of the various types of physical 
components for tensors of order greater than 2, it is more instructive to illustrate it by a special 
case. Consider the tress-strain relations in elastieity theory, 


Bar ne ah ll), 


where t', and ei, are the stress and strain tensors, both symmetric. We seek physical compo- 
nents for the tetradic C",*, of elastic constants of the material. By (5.10) we may write (11.1) 
in the form 


= Er or N ec De el er urn 8 E95 
9;; 911 

Since ti; and ei, are symmetric, use of left instead of right physical components would have 
led to a relation of the same form. Without loss of generality we may suppose (ÜRrl— (ik!, 
Cijkl— Oüilk, hut in general there are no symmetry relations connecting the first with the 
second pair of indices. From (11.2) it follows that suitable physical components for C,*, are 
given by 


Olijkı) Vene 0, ER nn ILL), 
9; 9ı 

These physical components will bear the dimension of stress in any co-ordinate system and 

allow us to write (11.1) in the form 

BE LTE TR ea, 5 (IA). 

It is the quantities occuring in (11.4), not those in (11.1), which the experimenter may measure 

directly in observing the elastic properties of cerystals. 

Another set of physical components is 


Veen 222222220002 019) 
Ir 9ıı 
These, too, are of the correct dimension, but they are not the components desired, because they 
do not satisfy (11.4) except in the special case of an orthogonal co-ordinate system. 
This example is particularly instructive, since as noted in $1 Ricciand Levi- 
Civita wrote 
an... nn(ik6), 


legitimate tensorially and equivalent to (11.1), but difficult to interpret physically. er 
To obtain the physical components appropriate to a given problem, one takes the defining 
equations and, starting with the tensors of lowest order occurring, replaces all tensor components 
for which suitable physical components are already defined by expressions in those components. 
The final result then indicates the proper definition of physical components for the quantity 
defined by the equations. In this way one always gets physical components which have the 
right physical interpretation and yet can be manipulated according to the rules of tensor algebra. 
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12. Covariant derivative (gradient) in the case of orthogonal co-ordinates 


All authors who have discussed this subject agree that when the co-ordinate system is 
orthogonal, physical components are to be defined by 


j r Iaü I Jipip ger ei 
ERBN An en Pers, 


Inner * Tigig (12.1) 
= Yo % Bir Gjgig a“ TIER 


a 
net 
Vi -- e 77 ’q : 


For simplieity restricting our attention to this case, we may find a rule for calculating physica 
components of the gradient of a tensor directly in terms of the physical components of the 
tensor itsellf. 


For the covariant derivative we have 
ee rasen en ee 
1 q lei, 2 ; ET . N ee x 
r d ER? Tü,, + [40 Piperj, Te (12.2), 


Oak 


the Tü,,beingChristoffel symbols. 
First we define 


a‘ DES HERE 


reelle (Anstmmed) aaıbing al Baal 
V9359«r 


These quantities are dimensionless in all co-ordinate systems. and in fact it is easy to show that 


1 es. — 
ae Yi- ee re Wi, 


or a; 00 
23 ‘ V9;; . (12.4). 
T(ij)= yo, Bail08 Y9; . T(ijk) =0 (unsummed, £, 7, k#) 
ii 
By (12.1) and (12.3) we may put (12.22) into the form 
ee 
Gıuü: LTR 
0) je was 9 Inn: Gas 1/Irr Iıı 
= — aliyz: jo) Ve ati... Ver ram+- 12:5: 
Ink | iin Iigin 1 zu. ER (iz. ..Ip) oh (i1kl) ( 
1 lie I; al Ela etalie dei a Inn Tütkj) —--- (in: j, unsummed) 
IGuür- 911 
whence 
re 1 data... RE 
@(i1.» .jp>k) Be 2 4 lo 08 Shih "Sgig 
V9rx 07 on 02 Inu Fi, 
+ allig.. 23) Tlink) + — air. iplier nd) Tikjy) —  (unsummed) (12.6). 


The method of forming the gradient indicated by (12.6) is expressed in 
Ehegten I The co-ordinate system being supposed orthogonal, to find the physical components 


alh.. k) of the gradient of a tensor whose physical components are all ...1,) , the following 
steps oil be carried out: 
1. Divide the indices I,...1, arbitrarily into two sets iy...Üps J1r++9g: (The result is 
independent of the choice of elements Jor the two sets, one of which may be Be if convenient.) 
2. In the tensor formula for the covariant derivative a" "'v,. ji, Make the: following 
replacements: 
I a" ... 'r, ig by aliy in) 
1 
11. i | 


Fr Be Perl? 
Il, 0m Er N imn pls 
where the latter symbols are defined by (12.3). 
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3. To the result of step 2 add the quantity 


ih ing: Ind Figig 
Ygz, 08" TERN PILZ 
FERNER ar 
The result of step 3 is ally* + -1,,k). 
As will appear in the succeeding examples, the term added in step three usually cancels some 
of the terms arising from the Christoffel symbols, with the result that expressions for 


the physical components are usually simpler than the corresponding fully explieit formulae for 
the tensorial components. 


As a first example we calculate the gradient of a vector. From theorem 12 we obtain 
either of the two results 


u 1 Al) Al) 9 ar 
Ali, j) = — = + —- — log / — +4) Tiüjn, 
Yg; 9 Yg;; 9 gi | 


m1oıy 
1 Ai 4 0 n_ ( 
= im er Kr 57108 Yan an) Dtji) (unsummed) | 
’ ji 
By (12.4), both forms reduce to 
1 Mi j — 
Aare GN annldt.uE ‚logYg; (I); 
Yg; 97 Ygi; 9 
Ali,j) = ms e er HRFL2:B). 
Pi a + BEN) TErF „108 Ygi; (© =) (unsummed) 
ii li hen 


These expressions yield at once Lame’s formulae for strain components in orthogonal curvi- 
linear co-ordinates, etc. 


As a second example we take the gradient of a second order tensor. Beginning with one 
of the three possibilities given by theorem 12, we obtain 


1 r) .. ade; 
— —— Beh E— a . log Vi 9; + ed) Tiikt) + acil) I'(jkl) (unsummed) (12.9). 
Yen 9% V9xr 9® 
and hence by (12.4) it follows that 
1 datij) (kj) 0 de od a: 
= AD eher rl Yrr (I,k#), 
nr * irre 


1 odalij) , alii) —al(jj) 9 — a(il) 0 — wa. 
a aml0eVe;+ 2, — log Yan; (k=7j,17%)), 


@(ij,k) 


a(ik 


Y9;; 07 Vi 127 Y9ıı 
1 odatlij) , ajj)—eatii) 9 — a(lj) 0 — BAAGUN 
Alt ik) = er ee log Yg: + — -logYgi: (k=1,i7j), | (12.10). 

Yi 98 V9;; 0x P2 Ynı 9% ( 
1 odatii) tik) -+akki) 0 — Kr: i 

—— > -logYr (@=j,kri), 

9er 02 Vi 0% 

1 o9atii) ala train 9 


z — I —_ JogYgi; (i=j=k) (unsummed) 
Ya; 97 FRi Yı 9 


By summing and combination one may deduce from (12.8) and (12.10) formulae for the 
divergence of a second order tensor, curl of a vector, second gradient ofa vector, Laplacian 
‘ofa vector, curl of the curl of a vector, etc. 


13. Coneluding remarks 


The problems set and solved in this paper were suggested by my own work in general 
continuum mechanics. I believe that with the apparatus and results given here one can quickly 
and easily put any tensor equation occurring in classical physics into a form in which all 
quantities therein have immediate and natural physical significance. Except for the remarks 
on eross-products in $ 3, all results and methods in this paper are validina Riemannian 
space of n dimensions. 
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Eingegangen am 30. Juli 1952, 


Die Spannungsfunktionen des dreidimensionalen Kontinuums 
und des elastischen Körpers *) 
Von Hermann Schaefer in Braunschweig 


Das Tensorfeld der Spannungen, dessen Divergenz durch die Gleichgewichtsbedingungen und dessen 
Rotation durch die Verträglichkeitsbedingungen vorgeschrieben sind, wird mit Hilfe eines Tensorpotentials, 
dem Tensor der Spannungsfunktionen, bestimmt. Für die kovarianten Rechnungen sind enge Analogien 
zur Einsteinschen Gravitationstheorie richtungsweisend. Im Sonderfall des isotropen elastischen, 
Körpers stößt man auf die Boussinesq-Neuberschen Darstellungen der Spannungen und Ver- 
schiebungen durch Potentialfunktionen. 


The tensor field of tensions, the divergence of which is prescribed by the conditions of equilibrium, and the 
curl of which is prescribed by the conditions of compatability, is determined by aid of a tensor potential, the 
tensor of the stress function. For Ihe co-variant computations, close analogies to Einsteins theory of 
grawitation are representative. In the special case of the isotropic elastic body, one obtains the Boussi- 
nesq-Neuber representation of tensions and displacements by potential functions. 


Le champ de tenseur des tensions, la divergence duquel est prescrite par les conditions d’equilibre, ei la 
rotation duquel est prescrite par les conditions de compatibilite, est determine a l’aide d’un potentiel de tenseur, 
c’est-a-dire du tenseur des fonctions de tension. Pour les calculs co-variants il faut se servir d’analogies 
semblables a la theorie de gravitation dEinstein. Dans le cas special du corps isotropique &lastique on 
rencontre les reprösentations dd Boussinesq-N euber des tensions et deplacements par des fonctions 
potentielles. 


TeH3opHoe Hose HANPSUKEHMÜ, PACXOKAEHHME KOTOPOTO 3ANAETCH PABHOBECHEIMH YCJIOBUAMH, 
& BUXPb — YCJIOBHAMH COBMECTUMOCTH, ONPeNenAeTch IIpH IIOMOLLMH TEH3OPHOTO HOTEHNHAIA, — 
TeH30pa PyHRIHH HanpsuKeHHga. Asa KOBAPHAHTHEIX PACYeTOB PYKOBOAATCH AHAJIOTAMH TEOPHH 
TAroTeHHa JüHnmreüma. B yacrHoMm cAIy4ae H3OTPONHOTO YIPYrOTO Tea HOAIYYAMTCA IPeN- 
eragıeuna Byenunecka-Heüdepa Ana HampspKeHui U CABHKOB IIOCPeACTBOM HOTEHIHAJIB- 
HBEIX PYHRIMM. 


1. Bezeichnungen 


In den nachstehenden Ausführungen werden durchweg rechtwinklige Parallelkoordinaten 


x% (* =1,2,3) benutzt. Zur Abkürzung schreiben wir = — 0;; Summenzeichen werden 
nr 


fortgelassen. Der Operator A bedeute wie üblich 0,0, = FF 
T ig 


1T= 


2. Aufgabenstellung 


Die Grundgleichungen der klassischen Elastizitätstheorie lassen sich bekanntlich in den 
Spannungen alleine formulieren. Von den Randbedingungen abgesehen und bei Abwesenheit 
von Volumkräften hat der symmetrische Spannungstensor o,;, den 3 Gleichgewichtsbedingungen 


gr = ln. Se 


*) Vorgetragen auf der Jahrestagung der GaMM in Braunschweig 1952. 

Nach meinem Vortrage machte mich Herr H.Richter, Haltingen, auf die Arbeit vonBloch, V.I: 
Die Spannungsfunktionen in der Elastizitätstheorie, Priklad. Mat. Mech., Moskva 14, 8. 415—422 (1950) [Rus- 
sisch] aufmerksam. Die Besprechung dieser Arbeit im Zentralblatt für Mathematik 39 (1951), durch Herrn 
Richter war mir entgangen. Inzwischen besaß Herr Richter die Liebenswürdigkeit, mir Einblick in die 
Blochsche Abhandlung zu geben. Bei derselben Zielsetzung ist die Bl o ch sche Formulierung und Darstel- 
lung (Nabla-Symbolik) des Problems so grundverschieden von der meinigen, daß ich die Veröffentlichung 
meines Vortrages nicht für überflüssig halte. Die Ne ubersche Darstellung der Spannungen, zu der man die 
Blochschen Ergebnisse vereinfachen kann, scheint im russischen Schrifttum nicht bekannt zu sein. 
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und den 6 Verträglichkeitsbedingungen (Beltramische Gleichungen) 


mM 
(m =Poissonsche Querkontraktionszahl; o = o,,) zu genügen [1]. 

Diese insgesamt 9 Gleichungen für die sechs Komponenten des Spannungstensors sind 
nicht unabhängig voneinander, vielmehr bestehen zwischen ihnen 3 Differentialidentitäten. 
Auf diesen Umstand kommen wir an späterer Stelle noch zurück. 

Wenn der Spannungszustand von einer Koordinate, etwa x,, unabhängig ist, gelangt man 
bekanntlich durch Einführung von Spannungsfunktionen — die Air ysche und die der Tor- 
sion — zu einer übersichtlichen Formulierung des Integrationsproblems [2]. 

Es erhebt sich die Frage, ob man im Falle des allgemeinen dreidimensionalen Spannungs- 
zustandes entsprechend verfahren kann [3]. Bei einem solchen Integrationsverfahren würde es 
sich also zunächst darum handeln, die Gleichgewichtsbedingungen (2.1) eines Kontinuums mit 
Hilfe von Spannungsfunktionen identisch zu befriedigen. Es ist plausibel, daß man hierzu 
mindestens 3 Spannungsfunktionen benötigt. (Nach einer treffenden Ausdrucksweise von G. 
Prangeist das Gleichgewichtsproblem des Kontinuums ‚‚dreifach funktional unbestimmt‘‘). 
Im Schrifttum der klassischen Elastizitätstheorie findet man die Ansätze vonMaxwellund 
Morera][4], die in einem noch später zu erläuternden Sinne komplementär sind. In beiden 
Ansätzen sind die Spannungen durch Summen zweiter Ableitungen von je drei Spannungs- 
funktionen dargestellt, und man kann zeigen, daß jeder Spannungszustand im Gleichgewicht 
mit diesen Ansätzen erfaßbar ist. 

Nachdem nun auf solche Weise die drei Gleichgewichtsbedingungen (2.1) abgegolten sind, 
hätte man mit einem dieser 3-Spannungsfunktionen-Ansätze in die 6 Verträglichkeitsbe- 
dingungen (2.2) einzugehen. Für die drei Spannungsfunktionen erhält man dann 6 partielle 
Differentialgleichungen 4. Ordnung, die wieder nicht unabhängig voneinander sind. Sowohl 
für dn Maxwellschen als auch den Moreraschen Fall sind diese Gleichungen aufge- 
stellt worden [5]. Sie sehen jedoch so verwickelt aus, daß ein Versuch, für sie eine allgemeine 
Integrationsmethode zu finden, von vornherein als hoffnungslos erscheint. 

Die Frage der identischen Erfüllung der Gleichgewichtsbedingungen durch Spannungs- 
funktionen ist von C. W e be rerneut aufgegriffen worden [6]. Er konnte zeigen, daß die Über- 
lagerung der beiden Ansätze vonMaxwellundMorera--ein Ansatz also, der 6 Span- 
nungsfunktionen enthält — bei einer orthogonalen Transformation seine Gestalt bewahrt, so- 
fern man verlangt, daß die 6 Spannungsfunktionen sich wie die Komponenten eines symmetri- 
schen Tensors 2. Stufe transformieren. Hieraus schließt man, daß jeder der beiden Ansätze 
vonMaxwellundMorera für sich nicht kovariant ist, womit ihre Unbrauchbarkeit für 
ein allgemeines Integrationsverfahren verständlich sein dürfte. 

Aber auch der kovariante Ansatz von C. Weber hat noch nicht die Form, die seine 
Weiterverwendung erfolgversprechend machen könnte. Wir werden uns deshalb nach rationel- 
leren kovarianten Darstellungen umzusehen haben. 


3. Die Verträglichkeitsbedingungen des infinitesimalen Verzerrungszustandes 
Aus bald ersichtlichen Gründen nehmen wir die Betrachtung der Verträglichkeits- 


bedingungen vorweg. 
Mit den drei Komponenten u; des infinitesimalen Verschiebungsvektors u wird der sym- 


metrische Verzerrungstensor &;; durch 

&5:=3 (94 %%) 2 ee Kr (3.1) 
definiert. Wir nennen ihn den symmetrischen Gradiententensor des Vektors u und schreiben 
symbolisch 

(&;£) = Gradi. 

Offenbar kann nicht jeder symmetrische Tensor 2. Stufe ein Verzerrungstensor sein. Dieser 
hat vielmehr den Verträglichkeitsbedingungen, einem System von Feldgleichungen 2. Ordnung 
zu genügen, die man am einfachsten durch Elimination der Komponenten des Verschiebungs- 


vektors in (3.1) gewinnt. 
Abweichend von der in der Literatur üblichen Weise führen wir die Elimination folgen- 
dermaßen aus. Zunächst bilden wir durch Linearkombination der Gleichungen (3.1) 


er treu) — Ir + 2) 
(65. =Kronekersches Symbol; e = e,,), woraus wir durch Divergenzbildung 
0% 247 nn 4 Au, BIER FERN UN A A A A (3.3) 
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erhalten. Die Anwendung des Operators A auf (3.2) ergibt 
— de, ++3(9 Ay +, Aw) — 409,9 Au =0 . -.... 0. . (8.4), 
und nach Einführung von (3.3) in (3.4) ist die Elimination bereits vollzogen. Die Verträglich- 
keitsbedingungen erscheinen somit in der Gestalt 
Ar: = — Act N Hıt RI Hmm = I ::: = - B5). 
Ihre linken Seiten bilden einen symmetrischen Tensor 2. Stufe, den wir mit A;; bezeichnen. 


Diese 6 Feldgleichungen für die e;, (oder auch für die &;,) sind nicht unabhängig voneinander. 
Eine einfache Rechnung bestätigt, daß 4 
[op 4A;r = 0 . . . © . . . . ° . . . . . . . (3.6), 
der Tensor A;, demnach divergenzfrei ist. 
Im Schrifttum sind die Verträglichkeitsbedingungen in der Gestalt 
R; k,im = 0; 0, E&mt 0% Om El wä 0% I, Eim Mm 0; Om Ekı —=0 PHiYRIer SEIT (3.7) 
angegeben. In dieser Darstellung bilden die linken Seiten einen Tensor 4. Stufe R;; 1, der je- 
doch wegen seiner Symmetrien 
Ram kim: Time Hk: Rirım > Bumik de (3.8) 
in drei Dimensionen nur 6 algebraisch voneinander unabhängige Komponenten besitzt. (Im 
zweidimensionalen Kontinuum bleibt nur die eine Gleichung R;3, 12 = 0). 
Im verzerrten Kontinuum hat das Bogenelement zwar die Gestalt 
ds? = (0,, 1.2 82) U dia ee Eee 


aber die Maßbestimmung ist euklidisch geblieben. Es muß also der Riemann-Chri- 
stoffelsche Tensor R;; m von (3.9) verschwinden, und dies führt bei infinitesimalen Ver- 
zerrungen auf (3.7). 

Man bestätigt leicht durch Ausrechnen der Komponenten (3.5) den Zusammenhang 


Aus — Ra+1, «+2; ß+1,ß+2 ae bear Fa 0, N set Ze Heuer (3.10). 


(Jeder Index größer als 3 soll durch seinen Rest nach dem Modul 3 ersetzt werden.) 
Aus R;;, /m entsteht durch Verjüngung der symmetrische Krümmungstensor 


BER ER EEE BRENNEN) 
und durch weitere, Verjüngung das Krümmungsmaß 
Ri Bryan alle. r0E 07 =» 10 M-Dau TER ZER LEN 
Wiederum am einfachsten durch Aufstellung der einzelnen Komponenten ergibt sich, daß 
A: 0 Beurer aeg een 


also der Einsteinsche Tensor des Bogenelementes (3.9) ist. Der Riemann-Chri- 
stoffelsche Tensor kann somit im dreidimensionalen Raume durch den Einsteinschen 
Tensor ersetzt werden. Die Einführung der e;;, an Stelle der e;; in (3.2) und (3.5) findet man 
bereits beider Einsteinschen Theorie schwacher Gravitationsfelder [7]. 

Aufschlußreich ist die Darstellung des Tensors A;; durch ein Matrizenprodukt: 


0 9; —05\ (Eıı Ei £ıs 0 —8; 0, 
(4;:) = - 0 0 Oıll&sı &o2 © TERN PET 203127, 
9 — ud & Ef) \-h 9 0 
Die Matrix des Verzerrungstensors ist rechts mit der Operatormatrix der Rotorbildung, links 
mit der dazu inversen multipliziert. Wir wollen deshalb den Tensor A;; den symmetrischen 
Rotortensor des symmetrischen Tensors &;, nennen und symbolisch schreiben 
(A;s)= Rot (eu) htm se valmı- Ida an eat. 


Daß nach (3.6) die vektorielle Divergenz von A;, verschwindet, bringen wir durch die Schreib- 


weise 
Div Rot (&;) = 0 je a SE NO ER 5 TUE (3.16) 


zum Ausdruck. 
Die Verträglichkeitsbedingungen (3.5), die wir jetzt 


Rot (u). 0 NEL Te er en 


schreiben können, sind die notwendigen und bekanntlich auch hinreichenden Integrabilitäts- 
bedingungen des Gleichungssystems (3.1). Dann und nur dann ist ein symmetrischer Tensor &;; 
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symmetrischer Gradiententensor eines Verschiebungsvektors u, wenn sein symmetrischer Ro- 
tortensor verschwindet. Es gilt also die Identität 


a ante (BB) 


und der Satz, daß jeder symmetrische Tensor, der (3.17) befriedigt, sich als symmetrischer 
Gradiententensor eines Vektors u darstellen läßt. 


4. Kovariante Ansätze zur Erfüllung der Gleiehgewichtsbedingungen dureh Spannungs- 
funktionen 
Die Gleichgewichtsbedingungen (2.1) des Kontinuums verlangen das Verschwinden der 
Divergenz des Spannungstensors. In Analogie zu (3.16) lassen sie sich identisch erfüllen in den 
Komponenten F,; des symmetrischen Tensors der Spannungsfunktionen durch den Ansatz 


(0; 5) == Rot (F; ) r ° . . . . . A . . . . . . (4.1), 
oder nach (3.5), wenn man die Spannungsfunktionen 


(F = F,,) einführt, durch 
9, = — Id; + Pat ID — IA MmDim 99: (AB). 


Aus (3.18) folgt, daß für einen vorgegebenen Gleichgewichtszustand die Spannungsfunktionen 
F;, bis auf einen beliebigen symmetrischen Gradiententensor Grad v bestimmt sind. Ent- 
sprechendes gilt für die ®,,.. Drei der sechs F};, (bzw. ®;,.), allerdings nicht beliebige, kann man 
durch geeignete Wahl des Vektors v zu Null machen. So wird man z. B. auf die nicht kova- 
rianten Ansätze vnMaxwell (fs; =Fıs = Fa =0)undMorera(Fı = Fa = Fa = 0) 
geführt. 

Auf eine kovariante Normierung der ®,; kommen wir weiter unten zu sprechen. 

Die enge Kopplung der Spannungsfunktionen mit den Verträglichkeitsbedingungen wird 
durch das Prinzip der virtuellen Verrückungen deutlich. Wir begnügen uns damit, diese Zu- 
sammenhänge zu skizzieren. Dabei sollen der Einfachheit halber die Verschiebungen an der 
Oberfläche des Kontinuums Null gesetzt werden. Im Prinzip der virtuellen Verrückungen kann 
man die Verzerrungen &;, unabhängig voneinander variieren, wenn man die Verträglichkeits- 
bedingungen, mit Lagrangeschen Faktoren versehen, dem Integranden des Voluminte- 
grals beigefügt: 


[Kein dei + Tin, ım öRir, ım) dV = 0 0 ee N Er, Se ee (4.4). 


Definitionsgemäß haben die I';;,ım denselben Symmetrien wie R;;,,m (vgl. 3.8) zu genügen. 
Nach zweimaliger partieller Integration des zweiten Summanden in (4.4) — die Oberflächen- 
integrale heben sich gegenseitig fort — erhält man, da die öe,, unabhängig voneinander sind, 

Okm — 0,0, Later . . . . . . . . . . . . . (4.5), 
und man überzeugt sich, daß wegen der Symmetrien von J';z,ım die Gleichgewichtsbedingungen 
erfüllt sind. Die 6 Spannungsfunktionen T';;,;m bilden einen Tensor 4. Stufe, aus dem durch 
zweimalige Divergenzbildung der Spannungstensor entsteht. 

Mit den Spannungsfunktionen ®;, bzw. F};, bestehen die Zusammenhänge 

Tin ım — Om Drı + Ökı D;m — 657 Dim dm D;ı ee A a (4.6) 
bzw. 
Bee UNSEREN. zur Sezuauyd, 144,7) 
T',;. ist die einmalige, 7’ die zweimalige Verjüngung von 7';,1m)- (4.6) und (4.7) gelten nur in 
drei Dimensionen. Im zweidimensionalen Fall bleibt nur die eine Komponente 72,13, die 
Airysche Spannungsfunktion. 

Jeder Verträglichkeitsbedingungisteine Spannungsfunk- 
tion zugeordnet. Die Spannungsfunktionen sind die Reaktionen auf die geometrische 
Zwangsbedingung, daß bei einer Verformung des Kontinuums die Maßbestimmung euklidisch 
bleibt. 


5. Die Spannungsfunktionen eines durch Volumkräfte belasteten Körpers 
Sind Volumkräfte X, vorhanden, so lauten die Gleichgewichtsbedingungen 
uallsijanii ss hiru.d ıpıyun (Bil). 
Eine kovariante Partikularlösung gewinnt man durch den Ansatz 
=D —5ND .::: rn: (2) 
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der in (5.1) eingeführt, die Bestimmungsgleichung 
ADEIFEIN OPML SE 


liefert. Mit einer Partikularlösung ®, dieser Poissonschen Gleichung kann man nach (4.3) 
die allgemeinste Lösung von (5.1) bilden: 


0 = — 49,+ 9 (dr + D)+ (a B;ı+ Di) — dr 9ılmPımt Di) - - - 4). 


Bei Abwesenheit von Volumkräften sind die ®; harmonische Funktionen, die in (5.4) fort- 
gelassen werden können, weil (5.2) dann einfach der Sonderfall A®;; = 0 des Ansatzes (4.3) 
ist. Mit den D,, sind nämlich auch die 9, 2,, harmonisch. 


6. Analogien zur Einstein schen allgemeinen Relativitätstheorie 
Im Fall fehlender Volumkräfte besteht eine enge Analogie der vorstehenden Entwick- 
lungen zur Einsteinschen Theorie der Gravitation im Sonderfalle schwacher Materie- 
felder. Die Maßbestimmung 


d2= (0,42 F.)an; de, ae em a 


weiche infinitesimal von der euklidischen ab. Der Spannungstensor ist dann gleich dem Ein - 
steinschen Tensor des Bogenelementes (6.1): 

ö;5 + 2 F;, entspricht den Schwerepotentialen g;, und der Spannungstensor o;; dem Impuls- 
Energie-Tensor der Einsteinschen Theorie, während unsere Gleichgewichtsbedingungen 
das Analogon der Erhaltungssätze von Impuls und Energie sind. Im Rahmen dieser Analogie 
entspricht (6.2) den Feldgleichungen der Gravitation. Der Integration der Feldgleichung wie- 
derum entspricht die Aufgabe, zu einem gegebenen Gleichgewichtszustande die Spannungs- 
funktionen ®;, oder F;, zu ermitteln. Die Einsteinsche Integrationsmethode für die 
linearisierten Feldgleichungen ist, auf diese Aufgabenstellung übertragen, kurz die folgende [8]: 
Da die F,, nur bis auf einen Gradiententensor bestimmbar sind — entsprechendes gilt für die 
®;, — kann man immer die kovariante Normierung 


Dal. ee MEN. 2 in 
erzwingen. Dann zerfällt aber (4.3) in die beiden Gleichungssysteme (6.3) und 
AO Be AUEHER . (6.4). 


Aus (6.4) folgt wegen 0,0;;, — 0, daß für jedes Integral von (6. N 2% ©. harmonisch ist. Man 
kann also durch geeignete Überlagerung harmonischer Funktionen nach Integration von (6.4) 
immer (6.3) erfüllen. 

Auch die kugelsymmetrische, im Mittelpunkt singuläre Lösung der Feldgleichungen der 
Gravitation für den leeren Raum, die als Schwerefeld eines Massenpunktes gedeutet wird, hat 
in unserem Problemkreis ihre Analogie. Ihr entspricht nämlich die Lösung 


PD, =1g Ya+ Ei DB, = Da = Di: = Dis = Da, = 0 
der Gleichungen 
Ad. =0; Pr —=0, 


die als Spannungsfunktionentensor eines dünnen, in der z;-Achse gespannten Fadens anzusehen 
ist. Im Gegensatz zur vierdimensionalen Relativitätstheorie ist jedoch die Maßbestimmung 
(6.1) im dreidimensionalen Raume außerhalb des Kontinuums euklidisch. In diesem Beispiele 
lassen sich deshalb die Spannungsfunktionen als Gradiententensor eines unendlich vieldeutigen 
auf der &;-Achsesingulären Vektors darstellen. 

Im Falle eines endlich ausgedehnten Körpers sind die wirklich auftretenden Spannungen 
Verträglichkeitsbedingungen unterworfen. Erst solche ‚‚Verträglichkeitsbedingungen‘‘ für den 
Impuls-Energie-Tensor, die den Aufbau der Materie durch die Elementarteilchen beschreiben, 
würden die Relativitätstheorie zu einer in sich geschlossenen Theorie machen. 


7. Der elastische Körper 


Wir kommen auf unsere eigentliche Aufgabenstellung zurück. Das Tensorfeld der Span- 
nungen des isotropen elastischen Körpers ist aus seiner vektoriellen Divergenz (Gleichgewichts- 
bedingungen) und tensoriellen Rotation (Verträglichkeitsbedingungen) zu berechnen: 


Div (0) = — Kinos walls IE Zn 
Rot (85) = 0m 10.1 Sn RER N BRITEN 
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wobei o;;, und e,, durch das H o o ke sche Gesetz 


1 


2 GE, = Or — eig 


ÖSBORIER TEN EN Ab FUN IBAN 1753) 
zusammenhängen. (@ = Gleitmodul.) 

Es sei an die bekannte Aufgabe der mathematischen Physik erinnert, ein Vektorfeld aus 
seiner skalaren Divergenz und vektoriellen Rotation zu bestimmen. Dort wird man auf den 
Begriff des Vektorpotentials geführt. Hier spielt der Tensor der Spannungsfunktionen die Rolle 
eines Tensorpotentials. 

Wir behandeln zunächst den Fall fehlender Volumkräfte (X, = 0) und schreiben (4.3) 
nach Einführung des Operators A* 


DE Fee (7A), 
Mit demselben Operator lauten die Verträglichkeitsbedingungen (3.5) 
a Bea ed u al7), 


und o;, und e;;, hängen nach (7.3) durch 

1 5 m 
DEFEKT 1 ad 
zusammen. Nun ist, wie man sofort bestätigt, für einen beliebigen symmetrischen Tensor H;; 


As He A AH ae Keraee e 7). 


Für ®,;, machen wir den Ansatz 
De NO ne en. ie, wi rlZ8), 


wobei wir uns das Verfügungsrecht über die skalare Funktion 2 vorbehalten. Aus (7.4) bis 
(7.8) folgt dann 


1 1 Mm 
— Ar (AP + 842 — ker o)-0 70, 
Jetzt wird über Qso verfügt, daß 
mM 
ee a a ee et A (7.10), 
und es bleibt 
AI nenn (711). 


Die AY;, genügen demnach den Verträglichkeitsbedingungen, lassen sich also mit Hilfe eines 
Vektors Am durch 


AP; = 0, Aw. + 0, Aw; — 6,750, Aw SW re ee ie (7.12) 
darstellen. Die Integration ergibt 
Vr = (, w.+ 0, Ww; — 0559 OT, (7.13), 


wobei 40;;, =0. Die ©; = 0,0;, sind dann ebenfalls harmonische Funktionen. Wir setzen 
(7.13) und (7.8) in (7.4) ein, wobei sich die Ableitungen der w, herausheben, und erhalten 


= — 2 +IEA2 + + RI — IHN. +. - (7.14). 


Sind Volumkräfte X, vorhanden, so lassen sie sich, wie ein Blick auf (5.2) und (5.3) lehrt, sofort 
berücksichtigen, indem man die ©; in (7.14) den Gleichungen 


unterwirft. Aus (7.14) berechnet man 
EA  S.anunwunnar nt. (7.16), 
und mit (7.10) wird 
Mm 
AR= 99 N ai 087.17). 


Für den Verschiebungsvektor u des elastischen Körpers gewinnt man aus (7.14) durch einfache 


Rechnung die Darstellung 


8. Schlußbetrachtung 
Die Resultate des letzten Abschnittes findet man bereits beiBoussinesq[9], der sie 
von den Differentialgleichungen der Verschiebungen ausgehend gewonnen hat. Auf demselben 
Wege sind sie von Neu ber [10] wiederentdeckt worden, der sie dann durch zahlreiche Unter- 
suchungen räumlicher Spannungszustände bekannt gemacht hat. 
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Daß unsere obigen Betrachtungen über ein Kapitel der klassischen Elastizitätstheorie in 
bekannte Gleise einmünden, kann nicht verwunderlich sein und war zu erwarten. Daß sie uns 
aber gerade auf die Boussinesq-Neuberschen Darstellungen von Spannungen und 
Verschiebungen führen würden, konnte nicht von vornherein vermutet werden. 

Man weiß, daß die klassische Mechanik, insbesondere die Hamilton- Jacobische 
Theorie, bereits das Gerüst der Atommechanik in sich trägt. Daß sich auch in der Elastizitäts- 
theorie, diesem klassischen Beispiele einer Feldtheorie, deutlich die KonturenderEinstein- 
schen Relativitätstheorie abzeichnen, scheint eine nicht unwesentliche Nebenerkenntnis unserer 
vorstehenden Untersuchungen zu sein. 
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Eingegangen am 9. August 1952. 


Der Wärmeaustausch zwischen einem geheizten Band 
und einer Konvektionsströmung 
Von M. Herbeck, Max-Planck-Institut für Strömungsforschung, Göttingen 


Es werden Lösungen bei verschiedenen Randbedingungen für die partielle elliptische Differentialglei- 
chung Tes+ Tnn =nTe aufgestellt, die unter gewissen vereinfachenden Annahmen das Temperaturfeld 
beim Wärmeaustausch zwischen einem geheizten Band und einer ihm parallelen Konvektionsströmung 
bestimmt. Der Vergleich zwischen Theorie und Experiment ergibt gute Übereinstimmung, er ist allerdings 
wegen der nicht genau bekannten Randbedingungen nur bei nicht zu kleinen Windgeschwindigkeiten mög- 
lich. Jedoch lassen sich umgekehrt mit Hilfe der Messung und der Lösungen der Differentialgleichung die 
Randbedingungen abschätzen, da die der Rechnung zugrunde liegenden Randbedingungen alle einen frei 
wählbaren Parameter enthalten. 

Solutions of the partial elliptical differential equation Tee + Tyan =nTe are stated for different boundary 
conditions. On certain simplifying assumptions, the equation determines the field of temperature in the case 
of heat exchange between a heated tape and a parallel convection flow. T’he comparison of theory and ewperi- 
ment shows good agreement. Because of the boundary conditions that are not exactly known, the comparison 
is possible only for not too small velocities of Ihe wind. But, as the boundary conditions, on which the compu- 
tations are based, contain a free parameter, it is, in Ihe reverse way, possible to estimante the boundary 
conditions by aid of the experiment and the solutions of the differential equation. 

Des solutions de I’ &quation differentielle partielle elliptique Tee + Tnn = n Te sont proposees pour de 
differentes conditions marginales. A de certaines suppositions simplifiantes cette &quation determine le 
champs de temperature a1’ &change de chaleur entre un ruban chauffe et un courant de convection parallele. 
La comparaison entre la theorie et l’ experiment montre un bon accord, mais, les conditions marginales 
n’&tant pas connues exactement, on ne peut faire cette comparaison qu’ a des vitesses de vent pas trop petites. 
Dans le cas oppos£, les conditions marginales peuvent ätre Evaluees, cependant, & l’ aide du mesurage et des 
solutions de I’ £quation differentielle, parce que toutes les conditions marginales sur lesquelles le calcul est 
base contiennent un parametre libre. 


YCTaHaBIHBAWTCA PeMeHHA A.1IA SILIHUTHUYeCKOTO AußdepeHlmalbHorO YPABHECHHA B YACTHBIX 
HPOHSBOAHBIX Tec+ Tyn=nT: Ipu pasımYHBIX Kp&AeBbIX YCJIOBHAX; YPABHeHHe 3TO OIpenel- 
AeT, B CIIY4AE HEKOTOPEIX YIPOINAHOIMHX IIPEAIOCBLIOK, TeMIEPATypHOe IONe IIPH TEIIOO0- 
MeHe MeiKıy HarpeBaeMof ‚IeHTOÜ MH KOHBEKIIHOHHBIM TEYeHHEM, IAPAJLIENIBHEIM ei. 
ÜpasneHne TeOpHH C OIBITOM YKA3bIBAeT HA XOpomee COBHANEeHHEe; B BUY ?Ke TOTO, YTO kpae- 
BbIe YCIOBHA B TOYHOCTH He H3BECTHEI, CPABHeHHe BO3MO:KHO IHM IIPH He CJIHIIKOM MAJIBIX 
CKOPOCTAX BeTpa. KpaeBsie YcıIOBuA, ONH&AKO, MOTYT ÖBITb OMeHeHEI O6PATHO IIPH TIOMONIM 
N3MepeHus, M pemeHuii AuddepeHnuasıbHorO YPABHeHHA, TAK KAK KpäeBble YCJIOBHA, HA KO- 
TOPbIX OCHOBBIBAETCH BEIYHCIEHHE, BCE CONEP:KAT IIPOH3BOJIBHO BEIOHPAeMEIH ONEePATOp. 


Problemstellung 
Der Zweck der vorliegenden Arbeit ist es, unter gewissen vereinfachenden Voraussetzun- 
gen für den stationären Fall den Wärmeaustausch zwischen einer geheizten Halbebene bzw. 


*) Ein Teil dieser Arbeit ist bereits in stark gekürzter Form als Vortragsauszug in ZaMM 30 (1950), 
Heft 8—9, S. 259, erschienen. 


| 
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einem geheizten Band der Breite b und einer inkompressiblen Konvektionsströmung parallel 
zur Bandoberfläche zu bestimmen. Das Geschwindigkeitsprofil habe die Form 


T 
ER Sorksragannan + Ih 


wie esim Wandabstand y einer idealisierten Reibungsschicht mit der Wandschubspannung 7, 
der Dichte o und der kinematischen Zähigkeit » des strömenden Mediums entspricht. Die 
geheizte Halbebene ist zwar physikalisch unmittelbar weniger von Bedeutung, sie wird im 
wesentlichen aus mathematischen Gründen behandelt, da, wie später gezeigt wird, bei Kennt- 
nis des Temperaturfeldes der Halbebene durch einfache Superposition zweier Halbebenen- 
Lösungen für jede Bandbreite das Temperaturfeld gefunden werden kann. 

Bild 1 zeigt die prinzipielle An- 
ordnung. Das Band bzw. die Halb- Strömungsvektor 
ebene sind in z-Richtung unendlich 
ausgedehnt angenommen, um ein 
zweidimensionales Problem zu er- 
halten, und schließen mit der ebenen 
Wand eines Kanales ab. Die Halb- 
ebene beginne bei 2 =0, das Heiz- f IE 
band liege symmetrisch zur y-Achse. ES ENT LBINNEN 

Die vorliegende Anordnung ge- Bild. 1 Prinzip der Meßanordnung 
winnt in neuerer Zeit meßtechnisch 
zur Bestimmung der Wandschubspannung eine ständig wachsende Bedeutung [1], [2]. Die 
notwendige Klärung damit verknüpfter Fragen gab Anlaß zu den folgenden Rechnungen. 
Näherungslösungen zur Bestimmung des Temperaturfeldes bzw. der Wärmeabgabe existieren 
bisher unter folgenden vereinfachenden Voraussetzungen [1], [2], [3], [4]: 

1. Die auftretenden Temperaturdifferenzen sind klein, sodaß die Stoffwerte als konstant 
angesehen werden können. Damit entfällt auch der Anteil der freien Konvektion. 

2. Die entstehende Reibungswärme wird vernachlässigt. 

3. Die Wärmeleitung in Strömungsrichtung wird gegenüber dem Konvektionsanteil 

vernachlässigt. 
Hier soll dagegen die letzte Voraussetzung fallen gelassen werden, also die Wärmeleitung in 
Strömungsrichtung, deren Vernachlässigung bei kleinen Wandschubspannungen und kleinen 
Bandbreiten nicht statthaft ist, berücksichtigt werden. Bei der Berechnung des Wärmeaus- 
tausches wird als weitere Vereinfachung angenommen, daß nur die Bandoberfläche Wärme an 
die Strömung abgibt, nicht aber auch die Seitenflächen und die Unterseite über das Zwischen- 
medium der Kanalwand. Damit ist für die Bestimmung der Wärmeabgabe nur die Kenntnis 
des Temperaturgradienten an der Oberfläche der geheizten Fläche resp. seines 
Integralwertes notwendig. Als Kennwert für die Wärmeabgabe soll dabei im folgenden die 
Nusseltsche Zahl Nu genommen werden, d.h. die Wärmeabgabe eines Heizbandes bzw. 
eines Streifens der Halbebene pro Längeneinheit in z-Richtung, dividiert durch die Temperatur- 
differenz zwischen Heizfläche und ungestörter Strömung und durch die Wärmeleitfähigkeit 
des strömenden Mediums. Da aber bei den praktisch benutzten Strömungsfeldern das Ge- 
schwindigkeitsgesetz (1) nur in Wandnähe gilt und eine Kontrolle vorhanden sein muß, ob der 
Temperaturabfall in dieser wandnahen Reibungsschicht erfolgt, werden auch Formeln für das 
Temperaturfeld angegeben, deren numerische Auswertung allerdings im allgemeinen 
zeimlich zeitraubend ist. 


Strömungsvektor 


1. Aufstellung der Differentialgleichung 
In den folgenden Gleichungen für das Temperaturfeld wird mit dimensionslosen Tempe- 
raturgrößen 


a nn an 8 one. (2) 
O5 = 9, 
gerechnet, wo 
[077 die konstant angenommene Band- bzw. Heizebenentemperatur in °C oder °K, 
u, die Temperatur des ungestörten strömenden Mediums, 


© (x, y) die von den Ortskoordinaten z, y abhängige Temperatur im Strömungsfeld bei 
Störung durch das Heizband bzw. die geheizte Halbebene bedeutet. 

Die allgemeine Differentialgleichung für ein stationäres Temperaturfeld ohne Wärmequellen, 

das durch Wärmeleitung und Konvektion bestimmt wird, lautet 


0:c, (wgrad T)— div (Agrad T) — 0 
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Cy spezifische Wärme des Strömungsmediums bei konstantem Druck, 

). seine Wärmeleitfähigkeit, 


m Strömungsvektor. 
Hieraus erhält man für den zweidimensionalen Fall mit einer nur in «-Richtung vorhandenen 
Geschwindigkeitskomponente 


Yeczy! 
17 () 
und der Forderung der Temperaturunabhängigkeit der Stoffgrößen die partielle ellip- 
tische Differentialgleichung 


nn 
mit 
( 


B=3Z, vepttn,). 4 Bin 


&,n sind wie 7 dimensionslose Größen, da ß!!? die Dimension einer reziproken Länge hat. Der 
elliptische Typus der Differentialgl. (3) bringt gegenüber den bisher vorliegenden angenäherten 
Behandlungen des Problems durch eineparabolische Differentialgl. [3] eine wesentliche 
Erschwerung. 

Für die partielle Differentialgl. (3) müssen Lösungen bei vorgegebenen Randbedingungen 
bestimmt werden. Aus Bild 1 sowie (2) ergibt sich für die Randbedingungen der Halbebene 
unmittelbar: 


DD Trturre Un 08 | i 
T= (,für.&. beliebig, 1.7 = O0.) an ou porn 9)» 
und für die Randbedingungen des Heizbandes mit der dimensionslosen Breite 1 = Pl!" b: 
T=1für ®<P/4, n=0,| R 
T:=0 für £sbeliebie, ınenoenle mi une a! 
Diese Randbedingungen sind jedoch zur Bestimmung des Temperaturfeldes 7 nicht ausreichend, 
es muß ferner T für n =0 und &< 0 (Halbebene) bzw. & > 12/4 (Heizband) bekannt sein. 
Da die Bestimmung dieser Randbedingungen mit zusätzlichen mathematischen Schwierig- 
keiten verknüpft ist, werden im Folgenden für den Temperaturverlauf bein =0,&< 0 bzw. 
&>12/4 nur plausible Annahmen gemacht. Es sollen hier vier verschiedenartige Randbedingungen 
untersucht werden, die die Eigenschaft haben, die Lösungen mathematisch verhältnismäßig 
einfach zu gestalten. Wegen dieser gewissen Unbestimmtheit in den Randbedingungen und der 
eingangs erwähnten physikalischen Vereinfachungen können die folgenden Ergebnisse noch 
nicht als abschließende Behandlung des vorliegenden physikalischen Problems angesehen 
werden. 


2. Lösungen für Temperatursprung am Rande 


21. Die Randbedingungen 
Als erste Näherung für den in Wirklichkeit kontinuierlichen Temperaturabfall am Rande 
der geheizten Halbebene oder des Heizbandes seieinTemperatursprung angenommen. 
Damit lauten die Randbedingungen bei der Halbebene: 
= 1 tur. | 
TEEN HEN (7) 
T=0 für & beliebig, n oo] 
und bei dem Heizband: 
T=.1.für&< 1/4, = be 
T=0 für 2>P/, n=0 Wr Dr SEE 
T=0 für & beliebig, n—® | 
Das Temperaturfeld des Heizbandes ist dabei leicht als Differenz zweier um + 1/2 und — 1/2 
verschobener Lösungen der Halbebene zu finden: 
Ty(E,n) = TE +Y2,) —- Te—Y2,n) 2. un u), 
wo der Index b auf das Heizband deutet, während die Lösung für die Halbebene ohne Index 
geschrieben ist. 
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2.2. Verbesserung derL&evöqueschenLösung 
Für die Randbedingungen (7) erhielt Lev&que[3]durch Vernachlässigung der Wärme- 
leitung in Strömungsrichtung über dieparabolische Differentialgleichung 
N ar.. (10) 
die Lösung 


1! 
3% 


m-1ı | 7 eo . (11) 
0 


im 1. Quadranten und 7° —0Oim 2. Quadranten. Mit (9) und (11) ist auch das Temperaturfeld 
des Heizbandes bestimmt. Aus (11) folgt für den dimensionslosen Temperaturgradienten der 
Halbebene an der Wand 


Bee nn: un (12). 
Hieraus berechnen sich die Nusseltsche Zahl eines Streifens der Halbebene der Breite 1 
1 
Nu = | (Taı=o d£ ER tee (13), 
wie auch eines Heizbandes der Breite / zu 
« x\1/3 
Nu= —0,807 ls — — 0,807 (® ”) Fur alas Hi 


d.h. die Wärmeabgabe ist proportional der Potenz 1/3 der Wandschubspannung und der 
Potenz 2/3 der Bandbreite. Dies stimmt für große I-Werte bis auf eine additive Konstante mit 
den bisherigen Beobachtungen [2] überein. 

Zur ‚‚Verbesserung‘‘ der Lösung (11) gemäß der Differentialgleichung (3) wird gesetzt: 


nee 7 7 es), 
wobei 7W der Differentialgleichung 
DPI TIE ST8,2 N para Ara al) 
und den Randbedingungen 
TW—=0 für n=0 und n—x genügen soll. 
Man erhält mit dem: Ansatz 
T® = & f(z), wo wieder 2= n/(9 E)® ist, im 1. Quadranten 


_ . p— + 2 
De ne | 154246232805. CRN m 13,438 >| 
3 


Ba te 
(17), 
ua N 


' 2 —- 


und im 2. Quadranten: TV —0, 


(1F, (a, ß, y) ist die Funktion von Kum- 
mer[5], S.111). Wegen &-#® in (17) 
ist diese Näherungslösung bei sehr klei- 
nen &-Werten nicht brauchbar! Der Ver- 
gleich mit der in 2.34 berechneten exak- 
ten Lösung, siehe auch Bild 2, wird 
zeigen, daß sie ungefährab > 1,n<£ 
anwendbar ist. Da 


(TP).-0= — 0,132 85% . . (18) 


ist, würde die Nusseltsche Zahl un- 
endlich werden, ihre Berechnung ist aber Bild2. Das Temperaturfeld der geheizten Halbebene bei Temperatur- 


Bi iılti 1 sprung am Rande 
auch wegender beschränkten Gültigkeit rra nF ung —_— Lösung von Levöque 
von (17) nicht sinnvoll ! — + — ‚verbesserte‘ L& v&quesche Lösung 


2.3. Lösungen mit Hilfe des Fourier-Integrals 
2.31. Die Integraldarstellung für das Temperaturfeld 
Da die Differentialgl. (3) die partikuläre Lösung 
T,= ei#+8) . (a +ian)!2- Zus + ion] Bares rl) 
26 
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hat, (Z1/3 = Zylinderfunktion der Ordnung 1/3; «,&, beliebige Konstanten), lassen sich, wenn 
längs der &-Achse gegeben ist und ferner die Bedingung T—0 für n—oo und beliebiges & 
existiert, mit Hilfe der Fourier- Integraldarstellung Lösungen formulieren. 

Für die Halbebene mit den Randbedingungen (7) ergibt sich 


1 1 eiat 
Te —hlan)da » 2.0... nn. (20) 


mit 


(a? +iom)!2- Ay ee (a? + va Di 


2 
el (Sala) 


Für Zı,s muß die 1. Hankel-Funktion wegen T=0 für n — oo genommen werden! 


h(a,n)= 421) 


2.32. Der Temperaturgradientander Wand(1.,= 

2.3212 (7 no tür Klermer Werte 

Wegen des großen Rechenaufwandes wird jedoch aus Formel (20) 7 selbst nicht berechnet, 
sondern nur der Temperaturgradient an der Wandoberfläche (7,)„-o , dessen Kenntnis sowohl 
für die Berechnung der Wärmeabgabe wie für die Bestimmung des Temperaturfeldes nach 
anderen Methoden notwendig ist. Bei der Halbebene ist 


a 2 
AH (Zalal 
1 (5 ) 
(T,)n=0 = lim >— eisen I. 5 dass: naurımcheid(22) 
ef a (Zalal) 


Die Schreibweise als Limes, die zur Konvergenz notwendig ist, wird im Folgenden fallen 
gelassen, bei der Auswertung von (22) und verwandten Formeln ist aber immer der entspre- 
chende Grenzwert genommen! 

Auf Grund der Umlaufrelationen der Zylinderfunktionen gilt 


(Tyn-—— | (Acosa&— Bsin ae) da. I nn tes ER 
7, 


u (303) 


en Mn 


Zur Durchführung der Integration werden 3 Integrationsbereiche unterschieden: 
a) Va 
BI 1 
Char ds 


mit 


In b) wird die Simpsonsche Regel angewandt unter Reihenentwicklung von sin«a& und 
cos &&. Die numerischen Werte für A und B werden mit Hilfe des Tabellenwerkes [6], das die 
Funktionen Zi) und 4), enthält, berechnet. 

Zur Integration in den Bereichen a) und ce) wird gesetzt 


2 
h(a) = MW (> fa 2). N N A (25). 
‚Für h(a) gilt die Differentialgleichung 


4 
On a 


8 
Man erhält mit dem Ansatz h(a) = = 4, a” eine asymptotische Lösung für h(«), wobei 
a, =i durch die asymptotischen Knemarkiuneer der Hankel-Funktionen bestimmt ist. 
Damit wird 

5% Eh, 15 11054 
3201 64a 1048 08 


hko)=A+iBri— at ... (27). 
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Darstellung (27) für h (a) unter Berücksichtigung der vier ersten Glieder wird im :Integra- 
tionsgebiet c) benutzt; diese Näherung bringt einen Fehler von weniger als 1% in die Berech- 
nung von (T,)„-o. Dabei ist 


1 Kerze 1 (cos2& ‚ 
=) Pe da = I 5 Hesi2e) EU I GE 


usw. 
Für den Integrationsbereich a) wird (26) durch die Transformation 


ander ehayeule). >... ea (29) 
in die spezielle Ricattische Differentialgleichung 
DE en N 2 u Ser EN) TE 2 . (30) 
übergeführt. 
Mit dem Ansatz u = 3 a, 2”, wo a, = — 0,36453 + 0,63138 i durch die Reihenentwick- 


n=0 


lung der Zylinderfunktionen bestimmt ist, erhält man für « < 1,8908 (die dem Nullpunkt 
nächstgelegene Nullstelle von Zn (3) ist & = 1,8908 e-'/2, siehe [6], S. XXXV]): 
h(a)= A+iB mit 


A = — 0,36453 a72/3+ 0,26576 a? — 0,11249 a2 
+ 0,019749 «1% + 0,007748 «1/3 — 0,006404 u$ 
+ 0,001335 #34 0,000569 a3 —.0,000481 a! . .-. =. . (31). 
+ 0,000103 &?*/? + 0,000043 &?®/? — 0,000038 «!! 


= 


+ 0,000008 at 3... 
B = + 0,63138 a-2/3-+ 0,46031 «2/3 

— 0,034206 «10/3 + 0,013421 a1*3 \ 

—0,002311 «22/3 0,000985 «2613 | De nun ul). 


— 0,000178 &?*/? + 0,000076 a®®®?+ ... 


Zur Integration werden die Reihenentwicklungen von sin x & und cos «& genommen und dann 
gliedweise integriert. Die Zusammenfassung der 3 Integrationsbereiche ergibt: 


als Is EP Ay &* 
(0 Et 8 100m alelong > lan 3 2 
mit 
%= — 0,381 G = 0103 
4, = — 0,0625 a, = — 0,056 
4, = — 0,059 Ag = — 0,093 
as= 0,062 2, 0,2 
4= 0,066 a 0,18 
d, = 0,041 dı = — 0,32 


Das 1. Glied in (33), welches das Verhalten bei sehr kleinen & bestimmt, entspricht dem 
Temperaturgradienten bei reiner Wärmeleitung ohne Konvektion. 

Mit einem Fehler von wenigen Prozent ist (33) mit der gegebenen Glieder- und Stellen- 
zahl anwendbar für — 16 <E<s2,. 


2.322. (T,)-o für'große $-Werte 
Da die Reihenentwicklung (33) von (T,)„-o für große £-Werte nicht brauchbar ist, wird 
noch eine asymptotische Lösung bestimmt. Durch schrittweise Zerlegung von A(o)und B(a): 


A = 4-—0,3645 a2 = A — 0,3645 028 + 0,2658 02° . 2.2.2... (834) 
B = B + 0,6314 021% = B + 0,6314 a 28 + 0,4603 a? . . 2» +... (35) 
und durch partielle Integration erfolgt: 
oo 1 x bs 
ı[ A cos a&da = — 0,269 |&|=4/3 — 0,066 |&| 37? — 78 | 4’ sinaE&da.. . (36). 
Are ö 
i 26* 
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Eine Abschätzung des letzten Integrals ergibt: 


1 Ar — 10.085 
2 4'"sina&da 5 eye: 
Entsprechende Rechnungen liefern für 
_. ı Bsin a & da = — 0,269 EU? — 0,066 &3/3 +0 >] : . (38). 
0 
Damit ist 
, 0,15 
(T,)n=0 = — 0,538 E43 — 0,132 8-58 +0|—; . (39) 
E>0 2r 
und 
0,15 
bei Temperatursprung bei exponentiellem (Die 0o= (0) Pe . . . . (40), 
Temperaturabfall c= 5 €E<0 


d.h. für £&<0 wird nur eine obere Schranke er- 
halten. Im Bereich 1S£< 2,5 stimmen Reihen- 
entwicklung (33) und asymptotische Darstellung (39) 
bis auf einen Fehler von höchstens 2,5%, überein. 
Man beachteferner, daß der Temperaturgradient der 
„verbesserten“ Lev&queschen Lösung (15) und 
(39) bis auf das Restglied übereinstimmen. Zum 
Verlauf von (7,)n-o siehe Bild 3. 

Im Folgenden sei der Temperaturgradient der 
Halbebene bei einem Temperatursprung mit 


(Da JE) ee ‚Ber en äh 


bezeichnet. Damit ist der Gradient des zur n-Achse 
symmetrischen Heizbandes 


Strecke 20 = I 


L 
| Temperaturabfall auf | 


& 


oT, 
()fe+VD EV) @) 
Te E -. [025 
wobei bei breiten Bändern auch für Er „ Mur 
= 
s h "< 1/4 
Bild 3. Der Temperaturgradient (Tn)n=o der Halbebene eine obere Schranke existiert, sobald für E<\/2 —1,5 
0,15 
ft& + 1/2) in die Größenordnung von 0,02 bzw. 172 —E)> oder darunter kommt. 
2.33. Die Nusseltsche Zahl Nu 
Die Nusseltschen Zahlen für einen Streifen der Breite Z der Halbebene 
l © 
L B 
Nu= [soa=- | Keosaı = + sina1z dass ae (AS) 
0 dr 0 $ z 
und für ein Heizband der Breite / 
y2 A 
vo |). -2 [osuin2as RE RE . (44) 
2 la" 


werden wegen des Wärmeleitungsgliedes a in (T,)n=o (33) unendlich. Setzt man jedoch 


Nu=,Nut,Nu, 22 Ge ee 


wo ‚Nu die Wärmeabgabe bei reiner Wärmeleitung, d.h. der Windgeschwindigkeit Null ist, 
»Nu diezusätzliche Wärmeabgabe durch Konvektion ist, so ist „N« endlich. Für kleine 
Werte von l erhält man für die Halbebene 


2 
„Nun — 0,381 1425 — 0,010 B + 0,003 8 — 0,002 Kin 1 +... Ara (40), 
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und für das Heizband 
I2 
NW — 2 + 0,007 4 — 0,004 Ylnl+... RBINIATHET]. (47)! 
Für I-Werte > 1 wird unter Benutzung von (39) ab&=1 bzw. £ BL = 1 für’den 
Streifen der Halbebene 2 

Nu = 0,253 — 0,807 12/3 + 0,198 ma + mi selen:- Kunney 6(48) 

und für das Heizband 
Nu, = 0,510 — 0,807 12/3 + 0,198 I— 2/3 + NP EIER URH NN AQ), 


Damit geht für große I-Werte in beiden Fällen „Nu asymptotisch in — 0,807 123, dieNusselt- 
sche Zahl der Lev&queschen Näherungslösung (14), über. 


Bild 4. Berechnete Nu-Werte als Funktion von 0,807 12/3 


a) Lösung von Lev&que 

b) Geheizte Halbebene mit exponentiellem Temperaturabfall am Rande für c = 6,12 
c) Wie b), jedoch ce = 0,612 

d) Heizband mit exponentiellem Temperaturabfallfür ce = 0,612 

e) Geheizte Halbebene mit Temperaturabfall auf endlicher Strecke 2D = 1,63 

f) Geheizte Halbebene mit Temperatursprung (ohne Wärmeleitungsanteil!) 

g) Heizband mit Temperatursprung (ohne Wärmeleitungsanteil) 


Wegen der unendlichen Größe von ‚N ist die vorangehende Berechnung von „Nu nicht 
recht befriedigend. Es zeigt sich jedoch bei etwas varriierten Randbedingungen, wo 7 keine 
Sprungstelle im Nullpunkt hat, und, wie erwähnt, durch das Experiment, daß sie zu einem 
brauchbaren Ergebnis geführt hat. Verlauf von ‚N« siehe Bild 4. 


2.34. Die Berechnung des Temperaturfeldes 


Für die Berechnung des Temperaturfeldes der Halbebene werden zwei verschiedene 
Methoden angegeben, als Erstes eine Reihenentwicklung in & und n. Mit dem Ansatz 


= n 
Dar Zn MEN 55 te Wopassort nn „=“; 30) 
wo 4 = T (8,0), = (T,)n= o ist, geht man in (3) hinein und erhält die Bestimmungsgleichung 
nenn _ı ne een ie Un Walt Fe We (51). 


Die Anwendung soll hier auf große positive &-Werte beschränkt werden, wo sie wegen der 
einfachen asymptotischen Formel für t, (vgl. (39)) besonders günstig ist. Hier ist 


Te 1 = 1—n (0,538 E-18 + 0,132 &-5/8) + 7? (0,0398 €-7/3 + 0,098 €-131%) 
nz IT sag ıT 2 CIUTEIRIET 
+ 7% (0,0149 &—213 -+ 0,0183 &-8/3) — 5(0,0155 &1818 + 0,084 8-17)... . . (52), 


wobei die Glieder (n&=1/3)* in der L&ev&queschen Lösung (11) enthalten sind, und die 
unterstrichenen Glieder bei der ‚‚verbesserten‘‘ Lösung (17) dazukommen. 

Bezeichnet man mit der Temperaturgrenzschichtdicke ö den Wandabstand, bei dem die 
Temperatur auf den e-ten Teil gefallen ist, so folgt aus (52) für & >1: 


. 1/3 
2124?) Bumsp In, BloydasnnEW Sijılandy(Bg), 
p 
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Zweitens gewinnt man eine Lösung für das Temperaturfeld bei kleinen Werten von &°+ n? 


. 
nach Transformation von e auf Polarkoordinaten r = (&+ n?)"? und 9=arc 187 : 


T Too —-sin®cosdr T,— sind T5.. ...% 2. (54) 
durch den Ansatz: 
Kae ICH udn nr] er FRE Haar 
n=0 
Mit Hilfe der aus (54) resultierenden, leicht lösbaren Differentialgleichungen 
Kt sin er d (na es sin 2 ah 2 ne ee 
tn m= IE 2 — 2) m—2+ SE SET OB = Se DT Re ea 


erhält man 


Tr, 9) =1 — +r:a,:sin®#+ r2|a,sin29-+ ZB an sin 2 2 
+73, sin3d+ ri 2, sin4d+ sin a — N = = | 
er nn 1%, sin5d-+ cos ® sin® 0 (58) 
+ r% a, sin 6 9 nn A | nn H cos? # sin! an +4) 
+ sin® der = = 4 a [sine a — 5 sin? d cos? o| +.. 


(58) ist dabei als Störung des reinen Wärmeleitungsfeldes T = 1 — ö/r aufzufassen. 
Die hier noch unbestimmten Koeffizienten @, müssen zur Erfüllung der Randbedingungen bei 
rsin® =n— m dienen, was jedoch in (58) "praktisch nicht möglich ist. Da aber 


= Ik: Et 
(= Fer: = (In)n=o Ne ein ke: a VE (59) 


ist und (7,),-o aus (33) bekannt ist, können durch Koeffizientenvergleich die @, bestimmt 
werden: 


Deals a, = 1/32 
a, = — 0,0095 @, = 0,0014 . (60). 
0.5=.0,0003 a, = — 0,0001 usw. | 


Es sei noch darauf hingewiesen, daß es und (7,)n-o nur dann bis zur Potenz r?"—1 
u £>0 

bzw. &2r—1, den Potenzen mit den In-Gliedern bzw. mit |&?”-1|, formal übereinstimmen, wenn 

von der asymptotischen Lösung (27) für h (&) das Glied mit «—?” berücksichtigt wird. 

Mit den Formeln (52) und (58) ist das Temperaturfeld für kleine 7 und positive & und 
kleine negative &-Werte bestimmt. Für große negative £ kann T näherungsweise gleich Null 
gesetzt werden, wie Bild 2 zeigt. 

Hieraus kann mit Hilfe der Superpositionsformel (9) das Temperaturfeld des Heizbandes 
in der Umgebung des Bandes bestimmt werden. Der Fehler, der bei breiten Bändern durch 
Unkenntnis des Temperaturfeldes der Halbebene für & < — 1,6 auftritt, fällt dabei praktisch 
nicht ins Gewicht, da dort immer 7(& + 1/2, n) > T(£ — 1/2, n) ist. 

Wie aus der Einleitung hervorgeht, genügt für das vorliegende physikalische Problem die 
Kenntnis von Tin dem beschränkten Gebiet. 


” 


3. Lösungen für exponentiellen Temperaturabfall am Rande 
3.1.Die Integraldarstellung für das Temperaturfeld 


Da die gesamte Wärmeabgabe bei Temperatursprung am Rande unendlich wird, während 
die beobachtete Wärmeabgabe bei den tatsächlichen, wenn auch nicht bekannten Rand- 
bedingungen endlich bleibt, sollen jetzt Wärmeabgabe und Temperaturfeld für exponentiellen 
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Temperaturabfall am Rande berechnet werden, worauf sich eine endliche Wärmeabgabe er- 
geben wird. Dabei gelten für das Heizband folgende Randbedingungen: 


et für 2° = 7/4, n=U, } 

T = erelle| —U2) led, a er. (le 

T=0 » £ beliebig, n— © | 
c ist eine positive Konstante. Weil die folgenden Formeln für beliebige c gelten, können durch 
E- 1/2 
Überlagerung eine Anzahl weiterer Lösungen konstruiert werden. Da & = .[=*) ist, fällt die 


Temperatur am Rande um so steiler ab, je größer die Wandschubspannung ist, wenn c inner- 
halb einer gegebenen Versuchsanordnung als konstant angenommen wird. 

Um wieder das Temperaturfeld für verschiedene Heizbandbreiten aus Lösungen der Halb- 
ebene zu berechnen, ist hier von zwei verschiedenen Ausgangslösungen der Halbebene 
auszugehen, die erste mit den Randbedingungen: 


T=i1 für E>0, n==(, | 
T=eet er: nl, 4 (62), 
T=0 „ £ beliebig, no “| 

und die zweite mit 
T=1—- ee IUsz& =D, n=U, | 
N) et, a Selen); 
T=V0 „ £ beliebig, n> | 


Die Randbedingungen (63) stellen also eine Ausnahme gegenüber der sonst als konstant 
angenommenen Temperatur der geheizten Halbebene dar. Die zu (62) gehörenden Lösungen 
erhalten im folgenden den Index d, die zu (63) den Index a. Damit gilt für das Temperaturfeld 
des Heizbandes 

Ders SH 2) let UaN)er) insert (6): 


Die Fourier- Integraldarstellung für das Temperaturfeld zu (62) lautet: 


oo 


1 Bf 1 In: se 
L=z+ K >) einsn(a, m) a ee er) 


DT ce— ia 


= 09 


und zu (63) 


1 1 1 —1 bA..9 4 | 
= bsrulkumantaz et he (21). . (66 
18 +35 I& ARE )e h(a, n) da h(a, n) siehe (21) (66) 


— 0 


(65) geht in (66) über, wenn ce durch — c ersetzt wird! Für c — oo gehen (65) und (66) in (20), 
der Lösung für Temperatursprung am Rande, über! 


3.2.Der Temperaturgradient (7,»-o 
Wie in 2. soll zuerst (7,,)„- 0, anschließend die Nusseltsche Zahl und das Temperatur- 


feld berechnet werden. 
Durch Differenzieren von (65) und (66) wird erhalten: 


oo 


oT, 1 da 37 
e =| FFzlesast@4— ac) —sinaseB+ uch] . (7) 


0 


oo 


Be) = RE) u | 2 [cos a&(®A+acB)—sina&d(eB—ocA4A)] F (68). 


on /n=0 m 

Differenzieren von (67) und (68) nach & ergibt die Differentialgleichungen 
ge haan u...» (69), 
a Erlen ans wor ie ce # - (70), 


wobei f(£) nach (41) die Lösung mit einem Temperatursprung am Rande ist, 
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Die Lösungen von g(£) und h(£) müssen für c— und &— ooin f(£) übergehen, deshalb 


08) = oe [ersf@)de . 2 WE 
und 
h&=ce-: [ercf(e)de. ee N 


Die Auswertung von (71) und (72) muß sich den Werten von & und c anpassen, es werden des- 
halb verschiedene Formeln bzw. Berechnungsmethoden angegeben. 

Da in (71) für —1,6 SE <s 1,5 und c >2 im wesentlichen das Verhalten von f(&) im 
Geltungsbereich von (33) entscheidet, gilt dort für & > 0: 


C 5 5 In |&| (c She 55 
a ee pe ec Ben. 
elz0r (1 73% E a an sap > 3: 192nd 
1 25 ( 5 1 ) 1 (@ ) Äh 2 e 
Ai, = tee 
+85 rer | j92 0 10246) 1024 nn pt er (73) 
mit 
3 
h„(c) De a, siehe (33) ur cn a 
Für &<0Oistin 1. Näherung: 
1 A r 
le) ee 8 su 
hinzuzufügen. 
Analoge Betrachtungen ergeben für h(&) im Bereich 


(die Integrationsgrenzen in (72) beschränken den Bereich nach unten!) eine (73) entsprechende 
Formel, in der ce durch —c ersetzt ist und le —1-+.c£) hinzugefügt ist. Für &<0 ist 
dieses Glied wieder zu streichen. (76). 

Das 1. Glied in (73) und (76): + = e=°® Eu( c &), welches das Verhalten für sehr kleine & 


bestimmt, entspricht dem Temperaturgradienten bei reiner Wärmeleitung für die gleichen 
Randbedingungen! 

Für kleine &-Werte und beliebige c-Werte, jedoch nur für kleine c-Werte von Interesse, 
erhält mandurchdirekte Auswertung von (67), wennim Integrationsbereich 0<a<2für AundB 
(31) und (32) bzw. Tabellenwerk [6] benutzt wird, und füra > 2 zur Herabsetzung des Rechen- 


aufwandes näherungsweise A» — und B*1 gesetzt wird: 


402 
oT, Fee 
en [eimtkl yY4+c)+ regt ha a  - (77) 
>0 
(In y ist die Eulersche Konstante) mit 
1 
H=tdn—z- 
ec? 
= — 1 pP 
2c 3° € 
nz Zn u  F (78) 
2.02. Io 
a a 
2 


& 
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Es ist 
2 
c [cAo®* — Ba’"+l 
de,2n ee 2 2 da 
T a + c 
0 ( 
4 nn a a a) 
ce [eBa?r1+ A 0?" 
d., ZN © 2 2 da 
TC a -+c 


Es braucht aber jeweils nur 


berechnet zu werden, da 


2 
Aa" " A a?n-2 j 
[28 B du = Aue | Garde EN) 


x 


usw. ist und AHA a?”r—2da usw. aus den Berechnungen für den Temperatursprung 


vorliegt. 
4A, = — 1,0749 Die 19251778 
A, = — 0,3237 Bez Ahlen 
A, = — 0,5780 BERIOBSÜN 0... 2: 
As = — 1,4576 B, = 32,453 
As = — 4,2381 B, = 103,69 


Für negative & ist in (77) wieder (75) hinzuzufügen. 
Die (77) entsprechende Formel für Temperaturanstieg wird erhalten, wenn ce durch —c 
ersetzt wird. Das Zusatzglied von (76) für £> 0 fällt hier weg! (82). 


Für & > 1 und beliebige c ergibt die asymptotische Entwicklung (39) von f(£) für (71): 
g(e)=— ce* fe (0,538 E43 + 0,132 8-5) JE+O[ET] - - - - - (83). 

Istauche£ >1, dann 
EEE 


. (84). 
— 0,132 &- Als: ve). it nt ce)? — a | 
Für (72) gilt wegen der Integrationsgrenzen nur für& > 1: 
& 
kK(&)=— ce! Zar es: (0,538 E13 + 0,132 EB) dE +0 [E?])- - - - - (85). 


Ist auch e &> 1, so gilt eine (84) entsprechende Formel, in der c durch —c ersetzt ist. Es gehen 
also beide asymptotische Lösungen (83) und (85), wie die Lösung bei einem Temperatursprung, 
für große £Eundc £indieL&v&quesche Lösung (12) über! 

Für den Temperaturgradienten des Heizbandes gilt: 


N N. HE 
on )„=o 
A } ? B . (07, : hrank « 
wobei, analog zu (42), bei breiten Bändern auch für re nur eine obere Schranke er 
AZ 
halten wird, sobald für & < a o(&+ =) in die Größenordnung von k(—1) oder darunter 


kommt. 
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3.3.Die Nusseltsche Zahl Nu 
Bei den Randbedingungen (62) berechnet sich Nu zu: 


1 ? da [cosal—1 snal 
Nu=-— il 5 (® B+oc Zn Gr] TER 
Ö 
Für die Randbedingungen (63) ist in (87) e durch — c zu ersetzen. (88) 
Für das Heizband mit (61) erhält man 
2 [cBda/[cosal—1 
v2 (Se 3 oma), en N 


Eine Entwicklung von (87) und (88) für kleine Werte von I gewinnt man sofort durch Inte- 
gration von (73), (76), (77) und (82). 
Für das Heizband (89) gilt für kleine !- und große c-Werte 


Nsälite ee 
on Hza)le ey Hop Er RT; 


5ln] 4 In @m—4 v2 ju+1 x (90), 
a ‘ AR ) mat! SRG 
16 7 > m! = en +1)! % 

wobei die h, aus h, (74) durch Streichen der Glieder mit @, „ entstehen. 

Für kleine /- und c-Werte erhält man: 
Alle ) 20212 3 ) 
Nu = — U nperla-+ o) IE ler 32) 
i (91) 


2c® N 25 0? 2ch 137 c 
MN - ei =. — N) oa —  ——— ee — 
el ss B)t Aln | a ++! a 60 2) 


2 


"onB 

Dt, Diem aa, 
2 2 
joaöte 


Für große l- und beliebige c-Werte ergibt sich für (87): 


© 12 
Nus= Nur(1) + ee £2/3 4 0,198 &?%) a agent 
3 1 
für (88): 
. € I 
Nus= Nur(1,5) + Bi [ &* (— 0,807 &2/®-+ 0,198 €) ae as 
—o 1.5 


und für (89) die (92) formal entsprechende Formel, nur ist für Nur, dem durch Reihen- 


entwicklung gewonnenen Wert von Nu, für 1 =1, bei (92) der der Halbebene und hier der 


des Heizbandes zu nehmen. (94).- 


Bei genügender Größe von I und« c I gehen (92), (93) und (94) in 
Nu =.Const ++0,80742/Aı» (14) vs dig 2.2 SH EI 


über, wobei die Konstante einen großen negativen Wert bei sehr steilem Temperaturabfall und 
einen großen positiven Wert < 0,807 12/3 bei sehr flachem Temperaturabfall hat, ungefähr bei 
ce = 1 ist die Konstante praktisch Null. Asymptotisch geht (95) wieder in den Lev&que- 
schen Wert (14) über. Zum Verlauf von Nu siehe Bild 4. 


3.4.Die Berechnung des Temperaturfeldes 


Wie in Abschnitt 2.34 kann auch hier das Temperaturfeld der Halbebene nur in einem 


Streifen längs der positiven &-Achse und für ein kleines Stück längs der negativen &- Achse be- 
stimmt werden. 


Für große positive &-Werte ist die durch Formel (50) und (51) bestimmte Methode zu be- 


nutzen, wobei für t, (62) bzw. (63) und für £, (83) bzw. (85) zu nehmen ist. 


EEE ON Okr INor; 10958 Herbeck, Wärmeaustausch zwischen einem geheizten Band usw. 375 


Bei kleinen Werten von &? + n? ergibt die Transformation in Polarkoordinaten mit Hilfe 
von (54) bis AL für die Halbebene mit exponentiellem Temperaturabfall 


n=1 


aD, cos nd-+Inrsin em 


ur 3 gm °0 Osin®d + sin ) + r? e sin? d — —— bin 49% — 4 sin 29) 
SE REN (cos 4# — cos Ra ind Inr 
144 77 Han 
2 (96). 
+r aa (sin 59 — 3sin 30 +2 sin®) +7 ” cos sint d 
: 396 
Eresr (cos 59 — cos 3 "FSoißr SoIsz sin3d — 1608 sin® 
BEE 4 
sin 59+,, „cos 9 sin ’) 
Hierbei ist T’ als Störung des reinen Wärmeleitungsfeldes 
 (er\®f 35h) 
D= Een. Inoren sinnd-+ cos nd — sin Da | lan&7: (97) 


aufzufassen. 
(98) Bei exponentiellem Temperaturanstieg erhält man eine (96) analoge Formel, wenn c 
durch — cersetzt wird und folgende Glieder hinzugefügt werden: 


N 3 n 2 2 
— 2, (— 9)" le 9 — Zsin 204°, sin3o) 
n=0 


sind+ 7 gsi) + ..:., 


Halo 


Das zu (98) gehörende reine Wärmeleitungsfeld ist: 


1 (-er)" :0E buxek pr 
T= l -sinnd+(d—n)cosnd—sinnd > —|. . . (99). 
=> = Inen sinnd+( z) cos N > ( 


Zur Erfüllung der bisher nicht berücksichtigten Randbedingungen bei r sin = Nn—> © 
müssen wieder die noch unbekannten Koeffizienten @,, die eine Funktion des jeweiligen e und 


N 
in (96) und (98) nicht identisch sind, durch Koeffizientenvergleich von =, mit der ent- 


sprechenden Reihenentwicklung von (T7,)„-o bestimmt werden, siehe (73), (76), (77) und (82). 
&>0 


Aus den obigen Lösungen für die Halbebene kann mit Hilfe von (64) das Temperaturfeld 
des Heizbandes bestimmt werden. 


4. Lösungen für Temperaturabfall auf endlicher Strecke 
4.1.Die Integraldarstellung für das Temperaturfeld 


In diesem Abschnitt sollen Lösungen für ein Temperaturfeld aufgestellt werden, bei 
denen T am Rande auf endlicher Strecke 2 D — Const. auf Null herabfällt. Für das Heizband 
soll dabei gelten 


T=-1 für @<ı/a, n=0, 
T=0 für 22 (1/2 -+2D)%, 7=0; 
—D—]/: D+1/2 
RE EN u Pre 23600). 
7 D D 


für B/AS &<(2D+1/2)%, n=0, 
T=0, für £ beliebig, Nn>® 
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Die entsprechenden Bedingungen für die geheizte Halbebene, die hier, um eine möglichst ein- 
fache Schreibweise zu erhalten, bei & = D beginnen soll, lauten: 


T=1 für ED, n=0, | 
1 

ea) : ? 2 2 an 
1 = arecoge/ D arme 770 =, en 
Le) für ED, Wal | 
T=0 für Z£ beliebig, Nn>&© R 

Die Integraldarstellung für die Halbebene hat dann die Form: 
nut et, D)h(«,n)da _h(a, n) siehe (21) (102) 
a onen, o(& a,n)da (a, n) siehe ae Ei j 


08 


J,ist die 1. Bessel-Funktion der Ordnung Null. (102) geht für D—0 in die Lösung mit 
einem Temperatursprung über, d.h. je kleiner D, um so steiler der Temperaturabfall. 
Für das Heizband gilt die Superpositionsformel 


= TE FW + D, N) TE V2ZD vr Bar ea 
Physikalisch könnte man sich den Temperaturabfall auf endlicher Strecke durch ein Heiz- 


band realisiert denken, das auf einem die Wärme schlecht leitenden Einsatzstück in der Wand 
eines Kanals mit Metallwänden aufgebracht ist. 


42. Der Temperaturgradient (7,)=o 
Durch Differenzieren von (102) erhält man 


nz | #6@D [Acosa&— Bsina&]lda . . :. . . . (104). 


Integriert wird nach analogen Methoden wie in 2.32 und 3.2 unter Benutzung der dort 
gefundenen Werte für A, und B, (81). Umjedoch den Rechenaufwand in erträglichen Grenzen 


zu lassen, wird im Integrationsbereich a>2: Ar und B=1 genommen. Damit 
wird für kleine & und D: 
sign & Bi 
0 = a pam Tsotr g reset, te 000) 
mit 
D: D4 D® 
h=ra+7 An+2 +4 er + 550g @n+6 Ha Sun. Le! BIER (106), 
wo 
a5 = — 0,3816 ai = — 0,0566 
a = — 0,05612 o&= + 0,0341 
af = + 0,02823 at = — 0,0586 
as = — 0,04533 ar—= + 0,144 
as = + 0,1061 = — 0,411 
alo = — 0,2932 aın=+ 1,27. 


Das 1. Glieh in (105), welches das Verhalten in der Umgebung von &? = D? für &2 > D? be- 
stimmt, entspricht dem Temperaturgradienten bei reiner Wärmeleitung. 


Für &2<D® ist ai und nn, von (105) abzuziehen und 
2\1/2 
ln) + are sin [EI ee ee Rn. Are (107) 


hinzuzufügen. 

Für D— 0 geht (105) prinzipiell in (33) über. Die hier vorhandenen kleinen Differenzen 
beruhen darauf, daß bei (33) in der asymptotischen Entwicklung für A und Bein weiteres _ 
Glied berücksichtigt wurde. Die Darstellung (105) für (7,)„-o ist praktisch bis DS 2 brauch- 
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bar. Ein besonderer Ansatz für große D wurde nicht gemacht, da dort am besten mit der 
Simpso.nschen Regel gearbeitet wird. 


Die asymptotische Darstellung des Temperaturgradienten für &>D-+1und&>2D 
(diese Bedingung folgt aus der hier nicht gebrachten Restabschätzung!) lautet mit Hilfe des 
Integrals vonSonine-Schafheitlin, siehe [5], S. 49, nach ähnlichen Rechnungen wie 
in 2.32: 


(T,)n-0= — 0,538 &13 . „P (2/3, 1/6; 1; De/£2) 
— 0,1328 5% - „F(1/8, 5/6; 1; D/E®) HOLE]... 2.2.2... (108). 


Für den Temperaturgradienten des Heizbandes gilt die (103) entsprechende Superposi- 
tionsformel. 


4.3. Die Nusseltsche ZahlNu. 
Die Integraldarstellung für Nu lautet bei der Halbebene 


& 


Nu—n [ aD) | (eo (1 + .D) — cos En a(l+D)— sin aD)|da (109) 
: | 


und bei dem Heizband 


Nu= [6 D) Zleos al + D)—cos aD]dh az nll0) 
0 
Die Integration von (109) ergibt für kleine Werte von / bei der Halbebene 
N 1 s De n &2 4 En+l D+Ll 


(.f„ siehe (106)) und die von (110) bei dem Heizband 


2 I a 212 
Nut pH lmt) \HA@+2DHH lat DIDI 4 2) 


Für große / ist bei der Halbebene und dem Heizband 


- 11 D:\ ]? +1 


5D 55 D 1309 2 


a), 


+ 0,198 ef nr 72 & ! 2268 & 104976 8 


wobei &, > 1+Dbzw. 2D zu wählen ist. 


Nur ist für Nur einmal der der Halbebene, das andere Mal der dem Heizband entspre- 
chende Wert zu nehmen. Für sehr große !-Werte gelten damit wieder (95) und die daran an- 
schließenden Bemerkungen, wenn an Stelle von c hier 1/.D gesetzt wird. Bild 4 gibt ein Beispiel 
für den Verlauf von Nu. 


44.Berechnung des Temperaturfeldes 


Der Bereich des berechneten Temperaturfeldes unterliegt den gleichen Beschränkungen 
wie in 2.34 und 3.4. Für große positive &-Werte wende man wieder (50), (51) an und nehme für 
to, = 1 und für t, die Darstellung (108). 

Zur Gewinnung von Lösungen für kleine Werte von &? + n? werden elliptische Koordi- 
naten o, r eingeführt, die durch 


Ee=DCosocose und n=DSinosint...... 2... (114) 


definiert sind, da die Benutzung von Polarkoordinaten, wie in 2.34 und 3.4, wegen der Rand- 
bedingungen nicht zweckmäßig ist. Die Transformation von (3) ergibt 


2 
Tao-tilr: -7 [sin 27 (Cos 20° —1) T,+ Sin2o0(cos2r—1)T,] . . . (115). 
Zur Lösung wird der Ansatz gemacht: | 


uam (116) 


= . : : hi. e h. h. 
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mit 
T 
Kur RE 


der Lösung der entsprechenden reinen Wärmeleitungsgleichung, und 
Real. Dsin tl4,1e0 


E23 <D: 


HE D D Sal 
I — —t(t ———— O8 —, N 118). 
Din e g,snrt z cos T 4m 7 cos T- > 7 D cos? T (118) 


Cos 2 und Sin 2c in (115) werden nach Potenzen von o entwickelt. Durch Potenzenvergleich 
sind dann die weiteren X, eindeutig bestimmt. 
K,=0 
= Ta7.D8 . = 
R,= 5 1550 — cos2r) — Ki 


D% SERELLLN 
K,= 57 [sin 27K,+ (cos27—1) Ki] 


Ki D: 
= 96 60 


(2 cos2r—1) 


Mit Hilfe der obigen Berechnungsformeln für die Halbebene und der Superpositionsformel (103) 
kann dann ebenfalls das Temperaturfeld des Heizbandes bestimmt werden. 


5. Lösungen für quadratisch exponentiellen Temperaturabfall am Rande 


5.1.Die Integraldarstellung für das Temperaturfeld 
Zum Abschluß sollen noch Lösungen berechnet werden, bei denen im Gegensatz zu 2.—4. 


oT(£,0) 


keine Sprungstelle hat. Hierdurch wird erreicht, daß (7,)„-o für jeden &- Wert end- 


0 
ee Als spezielle Randbedingungen für das Heizband werden angenommen: 
It für &2< 12/4, n=L, 
T= e-#(e—49' für & > 12/4, U a it 
T=0 für beliebige & no 


Um das Temperaturfeld für verschiedene Heizbandbreiten aus Lösungen der Halbebene zu 
berechnen, ist hier, wie in 3. von zwei verschiedenen Ausgangslösungen der Halbebene auszu- 
gehen, die erste mit den Randbedingungen. 


ya für &>0, n—0, 
T=e-t#& für E<0, ne . iur um Anz Se 
I für £ beliebig, n> © 

und die zweite mit den Randbedingungen 
T=1—-e# fur 0° nV, | 
Tl für &<S0, nl, le ER 
T,=0 für beliebige &£ n—m | 


Es bedeutet wie in 3. Index d Temperaturabfall und Index « Temperaturanstieg. Die Fou- 
rier- Integraldarstellung für das Temperaturfeld zu (121) und (122) lautet: 


oo 


Te S Auer Kae ic) h(o, n) da WE on (123), 

h(&, n)siehe (21) und R ur. | 1 | 
let 3 Kluge | 

Fl) 4% oena ER T ae a 


wobei in (124), wie in den folgenden entsprechenden Formeln, sich das obere Vorzeichen auf die 
Randbedingungen (121) und das untere auf die von (122) bezieht. Für d— erhält man (20)! 

Die Superpositionsformeln für das Temperaturfeld und den Temperaturgradienten des 
Heizbandes (64) und (86) gelten hier entsprechend. 
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9.2. Der Temperaturgradient (T,),-o 
Die Integraldarstellung für den Temperaturgradienten lautet 


00 


1 
TE E e-"td#(Asina&+BcosaE£)da 


0 


(Ty)n=o =7T 


. (125). 


1 ä l 2 3 a? 
a | (Acosa&£—Bsina£) " ll DE a) da 
Ö ; 
Die Auswertung von (125) wird nur für d> 1 angegeben, da für d < 1 eine weitgehende An- 
wendung der Simpsonschen Regel zweckmäßiger ist. 
Man erhält für die Randbedingungen (121), wenn für a>2 As=--—, Bz1 ge- 
nommen wird: 


d2 mM 


iR 
Für negative £ ist Er; R_ g(8) ($ Fehlerintegral) hinzuzufügen. b,, siehe (128). 
Für die Randbedingungen (122) gilt: 


d 1 
(Ta)ız rer Bild $). m u Ups 2-5 1 ; $as)- Haren 


m! 


Für negative & ist — & $(dE) zu streichen! 


Die beiden Ba Pe kichn in (126) und (127) entsprechen dem Temperaturgradienten bei 
reiner Wärmeleitung. Im Gegensatz zu 2.32, 3.2 und 4.2 sind diese Glieder aber nur für sehr 
große d-Werte für das Verhalten in der Umgebung von & = O allein maßgebend. 


Die Faktoren b,, in (126) und (127) berechnen sich zu 


= Im Am+ı |, Imt3 | Amts , Imtz | ) * , Imt2 Amts, Amis, 9, 
er 2 d n 4.d? 1 32 d3jih 3841) + mA 9.d2.:1,12 2 Do Tr kuadı 


Die Größen a}, siehe (106). 

Selbst wenn in I und (35) nur das 1. Glied Ban wird, haben die asympto- 
tischen Formeln für (7,)„-o eine ziemlich komplizierte Struktur: Summe und Produkte von 
Whittaker- -Funktionen, die für große &- und d-Werte in die Lösung für Temperatursprung 
(39) übergehen. Da der Aufwand bei der Auswertung der asymptotischen Formeln für d>1 
in keinem Verhältnis zum Gewinn an Genau keit gegenüber (39) steht, wird auf diese Formel 
hingewiesen. 

Wie eingangs betont, unterscheidet sich der Temperaturgradient hier sehr wesentlich von 
dem Verlauf bei den anderen Randbedingungen, da in (126) und (127) der Temperaturgradient 
für endliches d immer endlich bleibt im Gegensatz zu (33), (71), (72) und (105), siehe Bild 3. Es 
ist anzunehmen, daß dies dentatsächlichen Verhältnissen am ehesten entspricht. 


5.3. Die Nusseltsche Zahl Nu 
Die Integraldarstellung für die Nusseltsche Zahl der Halbebene lautet: 


BR; 


A(cos «1 — 1) — Bsin al] e= 144 da 
Vz f ( hab l 
(129). 


3 2 
4a . [ewsi-ntamen]; al errll 


Durch Integration von a und (127) über £in den Grenzen 0 his ‚laßt sich Nu mittelbar 
gewinnen! 
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Die a le für die Nusseltsche Zahl des Heizbandes lautet: 


1 2 1 a 3 a? 
U NUN ingle-ajı@ = ee ae Re 
Nu= am. Bsin al.e da+ 2 [nie al ] r |" 5° “)] da (130). 
y k 
Für kleine !-Werte ergibt die Auswertung: 
1 
2d 2 72 2 £2 1 2 78 = > im 
Nuz — —er#P | e®®de+— [eT®?. Bild?) —Inyd2l2] + 2[c0l4+ "+ cmı — + °°° 
Ya 7 m! 
. (131), 


wobei die c„, aus den d,, der Formel (128) durch Streichen der Glieder mit a$, gewonnen werden. 


Zur näherungsweisen Berechnung von Nu für große I-Werte, siehe Bemerkung an- 
schließend an (128). 


5.4. Die Berechnung des Temperaturfeldes 


Der Bereich des berechneten Temperaturfeldes ist hier auf kleine Werte von &+ n? be- 
schränkt, da der Temperaturgradient nur für kleine &-Werte bestimmt wurde. 


Nach Transformation auf Polarkoordinaten r, ® erhält man als Lösung für die Halbebene 
mit den Randbedingungen (121) nach der in 2.34 angegebenen Methode: 


14 Dar sinnd+ (—1) 


Nn= 


n 2n 
a dan 


——— (d cos2nd +Inrsin2nd) 


2 
+ Ai sin? ne 294 sin 49 —4sin 20); ,mrsin| 


+75 = cos #sind + --- (132). 


(133) Für die Randbedingungen (122) erhält man formal eine analoge Formel, bei der zu (132) 
noch 


4 
— 1-+ r? d2 cos 29 -+ gan 292 cos D + ..- 


addiert werden muß. Die unbekannten Koeffizienten @,, die in (132) und (133) nicht identisch 
sind, erhält man wieder durch Vergleich von ( mit (T,)n-o (126), (127). 
&>0 


Für große positive &, n & & kann die Lösung für Temperatursprung (52) angeschlossen 
werden. 

Das Temperaturfeld des Heizbandes folgt aus (132) und (133), bei breiten Bändern unter 
“Zuhilfenahme von (52), durch Anwendung der Superpositionsformel (64). 

Hiermit sei der rechnerische Teil der Arbeit beendet und abschließend auf die Messung 
der Nusseltschen Zahl eingegangen. 


6. Vergleich mit experimentellen Daten 


Messungen der Nusseltschen Zahl im Bereich 10 < 1 < 150 liegen vor [2]. Sie be- 
stätigen das asymptotische Gesetz Nu Const — 0,807 12/3. Diese Untersuchungen wurden 
durch eigene Messungen für kleine l-Werte (1 < 10) ergänzt. Zu diesem Zwecke wurden elek- 
trisch geheizte Pt-Ir-Bänder von 0,007 mm Dicke, 12—15 mm Länge und 0,09—1,2 mm Breite 
auf verschiedenen Unterlagen (u.a. Trolitul, Bakelit, Glas, lackiertes Aluminium) in eine der 
senkrechten Wände eines Luftkanals von 40cm Höhe und 5cm Breite eingebaut. Mittels 
einer Wheatstone-Brücke wurde bei konstanter Heiztemperatur, 30—50°C über Zim- 
mertemperatur, ihre Wärmeabgabe als Funktion der Wandschubspannung bis r=1CGS, 
hauptsächlich bei turbulenter Strömung, gemessen (siehe Bild 5). 

Bei der Bestimmung der an die Strömung abgegebenen Wärmemenge wurde die bei der 
Windgeschwindigkeit Null abgegebene Wärmemenge, die im wesentlichen durch die freie Kon- 
vektion bestimmt ist und von der Theorie nicht erfaßt ist, von der Wärmemenge beierzwungener 
Konvektion abgerechnet und diese Werte als Funktion von 12/? aufgetragen. 

Aus dem eben Gesagten folgt zwar, daß zusätzlich zu den eingangs erwähnten Vernach- 
lässigungen zur Vereinfachung der Theorie in einigen Punkten weitere Abweichungen hiervon 
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auftreten (u.a. endliche Bandlänge in z-Richtung, geringe Wärmeabgabe von den Seiten- 
flächen, Gültigkeitsbereich des Geschwindigkeitsgesetzes (1) nur in Wandnähe), jedoch zeigen 
Kontrollrechnungen und leichte Abwandlung der Versuchsbedingungen, daß hiervon nur bei 
sehr schmalen Bändern die Wärmeabgabe von den Seitenflächen ins Gewicht fallen kann. Die 
wesentlichen Ursachen dafür, daß es nicht möglich ist, die Wärmeabgabe eines Heizbandes bei 


2 


0,807 0% 


Bild 5. Gemessene Nu-Werte als Funktion von 0,807 12/3 


a) 0,09 mm-Band auf Cellon e) Berechnete Werte nach L&v@que 
b) 0,09 mm-Band auf Bakelit a*) 1,2 mm-Band auf Trolitul 
c) 0,09 mm-Band auf Glas b*) 1,2mm-Band auf Bakelit 


d) 0,09 mm-Band auf Aluminium, eloxiert 


kleinen Bandbreiten und Windgeschwindigkeiten theoretisch vollständig zu erfassen, liegen 
darin, daß 1. der Temperaturabfallam Rande nicht bekannt ist, und 2. die Wärmeabgabe durch 
freie Konvektion, die hier bei der Windgeschwindigkeit Null selbst bei der größten Wand- 
schubspannung von 1 CGS größer als die zusätzliche Wärmeabgabe durch die Strömung ist 
berücksichtigt werden müßte, da die Annahme einer einfachen Addition nur eine Näherung 
darstellen kann. 


Aus diesem Grunde sind auch hier nur Vergleiche mit der asymptotischen Formel (95) 
für Nu möglich. Dabei kann jedoch in Ergänzung zu den früheren Messungen [2] der Einfluß 
der Randbedingungen durch die Wahl verschiedener Heizbandunterlagen geprüft werden. 

1/3 

Bild 5 zeigt die gemessenen Nu-Werte als Funktion von 0,807 123 — 0,807 nn “2/8, 
wobei für die Stoffwerte c,, A, » die zu dem Temperaturintervall gehörenden Mittelwerte ge- 
nommen sind. Während bei dem 1,2mm-Band, dessen Breite bin «-Richtung mehr als 100 mal 
größer ist als die Dicke in y-Richtung, Gleichung (95) gut erfüllt ist, macht sich bei den 0,09 mm- 
Bändern die Wärmeabgabe an den Seitenflächen durch steileren Anstieg, der einer scheinbaren 
Bandverbreiterung entspricht, bemerkbar. Trotzdem wurde der Einfluß der Unterlage an den 
schmalen Bändern untersucht, da sie wegen der kleinen Temperaturgrenzschicht gewisse Vor- 
teile haben und eine Eichung mit bekannten Wandschubspannungen jederzeit möglich ist. Er- 
wartungsgemäß zeigen die Unterlagen mit größerer Wärmeleitfähigkeit, d. h., langsameren 
Temperaturabfalls am Rande, eloxiertes Aluminium und Glas die kleinste Wärmeabgabe 
durch erzwungene Konvektion. Dies gilt jedoch nicht für die gesamte Wärmeabgabe, die bei 
der Windgeschwindigkeit Null über der anderen liegt. Es besteht also die Möglichkeit, durch 
die Wahl der Unterlage bei kleinen Wandschubspannungen die Kennlinien weitgehend zu 
beeinflussen. 


Bei Annahme eines exponentiellen Temperaturabfalles am Rande hat bei den Messungen, 
Bild 4, die Größe ce Werte unter 1. Aus den geschätzten c-Werten kann ein Überschlagswert 
für die Temperaturgrenzschichtdicke ö mit Hilfe der Formeln für das Temperaturfeld gewonnen 
werden. Da der Temperaturabfall am Rande aus dem Verlauf von Nu nur geschätzt werden 
kann, genügt es, den einfacher zu berechnenden Wert der Halbebene zu nehmen, der etwas 
größer als der des Bandes ist. Der in Formel (53) gegebene Wert für ö stellt dabei nur eine 
untere Grenze dar, denn je langsamer der Temperaturabfall am Rande, um so größer die 
Grenzschichtdicke! 

Die Übereinstimmung zwischen Theorie und Experiment entspricht den Erwartungen. 
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Die Randwertaufgabe der Geoelektrik für Kugel und Zylinder 
Von A. Huber in Wien 


Die theoretische Grundlage der geoelektrischen Widerstandsmethoden bildet das im Unendlichen ver- 
schwindende Potential des von einer Punktquelle in einem inhomogenen Medium erzeugten Feldes. Es 
ist bestimmt durch die Ergiebigkeit der Quelle und durch zwei homogene lineare Randbedingungen entlang 
der Trennungsflächen elektrisch verschiedener Medien. Dieses Potential wird für eine kugel- und dreh- 
zylinderförmige Begrenzungsfläche konstruiert und für den durch einen Stollen mit kreisförmigem 
Querprofil realisierten Sonderfall eingehender diskutiert. 


The so-called geo-electrical resistance-methods are theoretically based upon the potential of the field 
that is generated by a point source in a non-homogenous medium. The potential — vanishing at the 
infinite distance — is determined by the sirengih of the source and by two homogenous linear boundary 
conditions at the contact surface of two media differing by their electric constants. The potential is 
constructed for surfaces of eylindrical and spherical symmetry, and is discussed in detail for the special 
case that is realized by a gallery with circular cross-section. 


La base theoretique des methodes de resistance geo-Electrique est formee par le potentiel du champ 
produit par une source ponciuelle dans un medium inhomogene. Ce potentiel, disparaissant dans l’infini, 
est determine par le rendement de la source et par deux conditions aux limites homogenes lineaires le 
long les surfaces de separation de medium de constants electriques differents. Ce potentiel est construil 
pour une surface de limitation de la forme d’une boule et d’un cylindre tournant et est discut& en detail pour 
le cas special r&aliseE par une galerie avec un profil transversal en forme de cercle. 


Teoperuyeckoi OCHOoBOÜ TEO3IERTPHYECKUX METONOB CONPOTHUBJIEHHA ABIAETCH HCYE3AHIIHH 
Ha ÖeCKOHEYHOCTH HOTEHIHAI MOJIA, MHOPOKAAEMOTO TOYEYHBIM HCTOYHHKOM B HEOAHOPOAHOH 
cpere. OH onpezenaertca ODHABHOCTLBD® UCTOYHHKA H ABYMA ONAHOPOAHBIMH JIHHeÜHBIMH KPAeBbIMH 
VCJIOBHAMH BOB IOBePXHOCTeÜu pasgeıa BNeKTpHyecKu pasılmyHsIXx cpen. Ilorenmnan 3ToT 
KOHCTPyHpyerTca Aa et HU BPAImaTeJIBHO UHJIHHAPHYeCKOÜ OTPAHHYHBAMIIHX IOBEP- 
XHOCTeH; TOAPOÖHee OÖCYKNaeTCA YAcTHBIa cAyyal, OCYINecTBIdeMEIH INMIOM KPYTIIOTO 
NOMepeyHOTO ceyeHuf. 


Die sogenannten Widerstandsverfahren der Geoelektrik zielen darauf ab, aus den an 
einem der Messung zugänglichen Ort festgestellten Spannungen in einem von einer Punktquelle 
erzeugten Feld Schlüsse zu ziehen auf die Lage der unzugänglichen Trennungsflächen zwischen 
elektrisch verschiedenen Medien und auf deren spezifischen Widerstand. Dazu muß natürlich 
die Funktion bekannt sein, welche die diese Trennungsflächen bestimmenden Größen mit den 
gemessenen verknüpft. Bisher kennt man sie nur in dem einfachen Fall, wo die Trennungs- 
flächen parallele Ebenen sind. Die Auffindung dieser Funktion auch für andere Begrenzungs- 
flächen bedeutet daher einen wesentlichen Fortschritt in der Theorie der geoelektrischen Wider- 
standsmethoden. In dieser Arbeit soll ihre Konstruktion für kugel- und drehzylinderförmige 
Trennungsflächen gezeigt werden, wobei die in einem Stollen vorliegenden Verhältnisse be- 
sonders berücksichtigt werden mögen. 


1. Die Randwertaufgabe 


Es sei der ganze Raum von zwei Medien mit den spezifischen Widerständen og, und p, er- 
füllt, die durch die Fläche 7 voneinander getrennt werden. Die in diesen beiden Raumteilen 
existierenden Potentiale U, und U, eines elektrischen Feldes genügen der Gleichung AU; =0 
(* =1, 2)und erfüllen entlang der Fläche r die folgenden Bedingungen: 


DEU, anereiea  e 


OU 0Ur 
02 Fr u! On" 


\ 


(@ 
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Da das uns hier interessierende Feld von einer Punktquelle @ erzeugt werden soll, deren Er- 
giebigkeit J in der Geoelektrik eine maßgebende Rolle spielt, liegt der folgende Ansatz nahe, 
wenn Qim Medium (g,) liegt: 


rose 1 habt; 
De tu), Dee IE LORD NSL EEE INN 79) 


wobei R die Entfernung des Feldpunktes P vom Quellpunkt @ bedeutet und %, ein Integral 
von AU = ist, das im allgemeinen in der Umgebung von Q beschränkt bleibt und den Be- 
dingungen (1) und (2) entsprechend zu bestimmen ist. Es ist die aus U, und U, zusammenge- 
stellte Funktion also einfach die in gewisser Weise normierte Greensche Funktion dieses 
Randwertproblemes. 

Liegt jedoch der Quellpunkt Q auf der Trennungsfläche 7, was in der Geoelektrik meistens 


J 
der Fall ist, dann ist der Faktor —* in (3) folgendermaßen abzuändern. Betrachten wir dazu 


4 
gleich den auch praktisch vorkommenden Fall, — man denke etwa an das Ausstreichen einer 
Schichtgrenze an der Erdoberfläche! — daß drei Medien mit den spezifischen Widerständen 


q(l=1, 2,3) entlang einer Raumkurve ( aneinanderstoßen, durch die also ihre Trennungs- 
flächen 7, hindurchgehen, und es liege der Quellpunkt @ auf der Kurve ©. Die Tangential- 
ebenen der 7; mögen bei Q die Winkel 277g, mit I q; = 1 untereinander einschließen. Da ent- 
lang C insbesonders (1) bestehen muß, so gilt mit der in allen drei Medien gleichen Spannung U 
in einer hinlänglich kleinen Entfernung r von Q entlang (©: 


RETTEN 
4nqg r 


(=1,2,3), 


wobei J; der auf das Medium (o,) entfallende Anteil des aus Q austretenden Stromes J ist, also 


J= N J,. Daher ist 
SEE 2 eh 
el og An-r ee 


so daß statt (3) zu setzen ist 


4 tl > 
SAUREN 
An.yg \& 

01 
Ist insbesonders @ ein regulärer Punkt der Trennungsfläche bloß zweier Medien, so hat man 
daher wegen g, =g, — 1/2 statt (3) zu setzen: 


JÖ4" 08 (4 “) 7 
N 6 Lie lach iiefegngig afiie lt) gar 
3 (oo AR Er 


Bildet aber die Trennungsfläche einen Knick mit etwa aufeinander senkrechten Tangential- 
ebenen, — man denke z.B. an die in einem Stollen zwischen Sohle und Ulm vorhandene Ver- 
schneidung! — dann ist , =4, 9, =? und an Stelle von (3) tritt nun: 


I’ 01°0 1 Eee 
U, =. 122. e +). 

a (4 +3) \R 
In derselben Weise ließe sich auch der Fall erledigen, daß Q in einem Schnittpunkt der Tren- 
nungsflächen mehrerer Medien liegt, doch wird man dazu praktisch kaum einen Anlaß haben. 
J 0% 
4 
analoger für das Zusatzpotential %; belanglos, so daß wir statt U; die Funktion 
Ar: U, En 
ET 


Da die Randbedingungen (1) und (2) homogen sind, ist der Faktor oder ein ihm 


ee 
“=, tM 


betrachten können, die das reduzierte Potential heißen soll. 
Wir wenden uns nun zur 


2. Lösung der Randwertaufgabe für die Kugel 
Sie liege in einem unendlichen Medium mit dem spezifischen Widerstand go, und sie selber 
von einem Medium mit dem spezifischen Widerstand o, erfüllt. Ihr Radius sei gleich der Längen- 
einheit und ihr Mittelpunkt ist der Ursprung eines räumlichen Polarkoordinatensystemes. 
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Legen wir den Quellpunkt Qin die z-Achse, dann sind die Niveauflächen offenbar Dreh- 
flächen mit der z-Achse als Achse, das Potential im Feldpunkt P hängt also dann nur ab von 
MPR=r und = 9MR=%: 

Wir behandeln zunächst den Fall, daß 

a)Q@ außerhalb der Kugelliegt, also MQ=a> 1, und machen mit der zonalen 
Kugelfunktion n-ter Ordnung P,(®) den folgenden Ansatz: 


1 oo 
ur d)=t DB BORN (r>1) 
n—=0 
1 Sr 
u(,d)=-n + Don. Pro). .... (r<1). 
n=0 
Wegen (1) folgt daraus sofort fürr =1: 
On — An; 
und dafürr<a: 
Sie 
RE 


ergibt (2) die folgende Gleichung für a,: 


rate +l)al=nr: nn 


oder mit der in der Geoelektrik üblichen Bezeichnung 


Drehen 
0% +01 
> ax n 
"Ton +IiFr)-art 
Somit wird: 
> PER —n—1. 
ud) + 2 nee (ar) A nr 3) 
1 DIE > y\n , 
u,(r, Es ae dere 23 en ml) Pa) @<l)e 2.2... A) 
Aus (4’) sieht man sofort, daß für 0, = 0, also x = — 1, also wenn die Kugel ein vollkommener 
Leiter ist, %, = 0 wird, wie es ja sein muß. a 
Um einen 


b)Vergleich mit einem homogenen Feld zu ermöglichen, darf man in 
(4) und (4°) nicht so ohne weiteres einfach a über alle Grenzen wachsen lassen, sondern man muß 
vorher die U, anders normieren, indem man sie ersetzt durch 


An a®U, 5 .( 2) eo 
ur (2,0) = r a en = 12): 


1 
Man erhält sodann zunächst für r>1 TER 


a? 


S-arreosd+Z TED TE EETzDE 


uN) = im ut, 0) = Im (ga +20 [TE ante] 
1 2,# 2 ® 4 
trage 4258. >): 


also die bekannte Darstellung des Potentiales eines durch eine Kugel gestörten homogenen 
Feldest). In entsprechender Weise ergibt sich auch für r< 1 der bekannte Ausdruck für das 
Potential eines homogenen Feldes im Inneren der Kugel. Die Abb.1 zeigt den Verlauf der 
Niveaulinien für den Falla =2 und o,/0 = 10, also x = — -%. 


1) J.N.Hummel, Beiträge zur geoelektr. Methode, Ztschr. f. Geophys., 4., 1928, 181. 
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c) Q innerhalb der Kugel,also«@ <1. Mit demselben Ansatz wie vorher 


ergibt sich nun wegen 
1 Shi 
A (a<r): 


R n=9 
nat +1) 2])=(e 9)’ (n-+ 1): a” 
oder: 
ee B 2.(n-+1):x-a% 
bs an+1l-+x 
und daher: 
1 x n-+]1 a" 
REDE AL RE et are er, sr 5 
ET Drnlahe Areeen: u ee (5) 
1 < n-+1 
WR ee ; Br iin, ee A! (57 
ur d)=,— 2% re (ar)"- P,(9) r<1 (8) 


Bild 1. Verlauf des von einer Punktquelle erzeugten Potentiales in der Umgebung einer 
besser leitenden Kugel 


Für 0, =0 liefert eine einfache Rechnung: 


J. 
I 
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Ferner ist leicht zu sehen, daß die entsprechenden Potentiale U, in (4) und (5) ineinander 
übergehen, wenn 

d) Qauf der Kugel liegt. Istnämlich « =1so kann man zunächst in (4) das 
YR mit der Summe zusammenziehen und erhält: 


De 2xn \ P,@ 
ul. 0) =. Sili4, 2). 
. n=0 


a2n+1il+x 
oder wegen 
1 Zune RENT Et . 
2n+l+“ an ıer Tat 2 
eo 0s al. Bit 
De Bam 


TR 


Ebenso ergibt sich aus (5) : 
| 2%x.(n-+1) 2.01 2n-+1 
2n+l+x Ga+a% 2n+1l+x 
und damit dieselbe Darstellung für U, in Übereinstimmung mit (3’). In entsprechender Weise 
zeigt man, daß auch die beiden Ausdrücke für U, ineinander übergehen, wenn a—1. 

Wenn sich der Feldpunkt dem Quellpunkt nähert, wachsen wegen R— 0 natürlich die 
Potentiale über alle Grenzen, aber die Zusatzpotentiale %,,bzw.%, bleiben, wie aus (4) und (5) 
sofort zu sehen ist, im allgemeinen endlich. Nur wenn Q auf der Kugel liegt, wachsen bei An- 
näherung des Feldpunktes an ihn auch die %, und %, über alle Grenzen, da ja ?„(0) =1. Dies 
leuchtet auch aus physikalischen Gründen ein, denn die Zusatzpotentiale rühren ja von der 
auf der Kugel befindlichen ‚‚freien‘ Ladung her, deren Flächendichte in der Umgebung des auf 
der Kugel liegenden Quellpunktes selber über alle Grenzen wächst. 


Bevor wir die 


3. Lösung der Randwertaufgabe für den Zylinder 

in Angriff nehmen, ist es zur Vermeidung von Unterbrechungen der betreffenden Erörterungen 
zweckmäßig, vorher einige 

a) Hilfsätze über Zylinderfunktionen zusammenzustellen. Sie 
betreffen besonders den Fall des Additionstheoremes, wo Feld- und Quellpunkt von der Zylin- 
derachse gleich weit entfernt sind, der 
aber bei der Ableitung desselben aus- 
drücklich ausgeschlossen wird. Mit die- 
sem Sonderfall hat man es aber gerade 
zu tun, wenn man die theoretischen 
Grundlagen für geoelektrische Messun- 
gen in einem Stollen mit wenigstens 
nahezu kreisförmigem Querprofil ent- 
wickeln will. 

Es sei: R?(a, r,9,z) = Rö(a, r, p) 
Bild 2. Zum Additionstheorem der Zylinderfunktionen +.2—r? BE a —2arcos p 4 2? und 

ar, dann ist bekanntlich: 


Bild 2a Bild2b 


ge 2 (I Sa La 

m = R R . t » RN ° = — N n . . . en 

R | Ko(Ro tt) costz- dt | Bei SER: | cos n | cos tz’ dt asr 
: E (6) 


wobei I,(z) =i*:-J,(i x) und K,(z) =3% "+1. H/®(ix) die sogenannten modifizierten 
N bedeuten. Die Frage, ob man in dieser Darstellung zur 
u unter dem S-Zeichen nach r differenzieren darf, wollen wir durch die 


folgende Überlegung entscheiden. Es sei zunächst «> r (Abb. 2a), dann ist: 


Bildung von 


OR, r—a:cosp Pi: 
RT Tee R, rn008(Rirugire Sr (7) 


wobei PQ@’ in der Richtung des wachsenden r als positiv gilt und + spitz oder stumpf sein 
kann. Ist aber a<r, dann hat man offenbar (Abb. 2b): 
oh, Q9'P 5 » 
ee ei a er 79) 


N 
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zu setzen. Nun gilt aber für « > r das folgende Additionstheorem?): 
Er Re) I Ile Dokarla:i)-err, 
Nn=—o 


woraus nach Multiplikation mit €?, da I„=I,und21, =1,1+ T,+ı, wegen (7’) folgt: 


oR = 
cs +0) K, (RR Ü)=— Fi KR t)= 3 Inı(rt): Ku(a:t)cos np 
= 19:6, + sd 1 + Inrı) En cosno= 3 In(r-t): K.(a-t)- cos NY. 


Da aber in dem mittleren Ausdruck von (6) die Differentiation nach r unter dem [-Zeichen 
sicherlich erlaubt ist, so hat man tatsächlich: 


al/R) 2 or, | | 
a t K,(R, t) -cost2' dt 

A (8) 

2 Se 
ns | n(r t)- Kn(a*t) cosngo|-costz-dt---a>r | 

Ö =—o 
Auf eine analoge Weise folgt aus dem für « <r gültigen Additionstheorem: 
Re Der la teK..(rtjvene, 

R=— oo 


daK„=K, und —2K,„=K,„ı+t Kyarı, wegen (7): 


o(1/R 2 s S 
922. 2 Ina) Kufr-t)>cosmg|: cos tz: de: Er. il 
r 7 n=— © 
Für « =r schließlich sei zunächst bemerkt, daß: 
oR,(a, r) _ 1 oRser, n)| TERN) 
or ar OT VOaRr- VER AR — 


Indem wir nun im sogenannten Watsonschen Integral °): 


2 2, hr 
Jar ft): K,(r ir Jn@r'®), dx 
Vi: + 22 
0 
nach t£ differenzieren, kommt: 
fe Inl2r: a) 
———— — de. 

2 213/2 
| 2) 
Multipliziert man beide Seiten mit cos n 9, so liefert die Summation über n wegen der Jaco- 
bischen Gleichung?®): 


cos (27 2: sin 2) - > Jonl27°8%)-cosno 


n=— 


d I; 
mn’ En) EN 


und der bekannten Darstellung: 


2] 


E41 


cos(2r- sin®.o) 
- de: 


Kular: sin® 1) = 


Vi? + a2 
0 22 
King cos(2 7: sin 2...) 
r 
EDS Er): ana t  (? + a2)R : da 


0 
BRRTAN | uud ) 
=—2r:sin--K,l2r sin *t]. 
rsinn'Kı 2 9 
2) Vgl.z.B.R.Weyrich, Die Zylinderfunktionen und ihre Anwendungen. Leipzig u. Berlin 1937, S. 90ff. 
»3)G.N. Watson, Atreatiseonthe Theory of Bessel Functions, Cambridge 1922,S. 435, (3) ;W. Magnus 
u.Fr. Oberhettinger, Formeln und Sätze für die spez. Funkt. der math. Physik. Berlin 1943, S. 30, 
%W.Magnus u.Fr.Oberhettinger,l.c. S.18. 


>3 [Ile Ku(r-leos ” "cos I d=—2sin®. fe Kr (27: sin 07 ) 
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Daraus folgt aber sofort die gesuchte Darstellung von 


o/R)| 
or Ir=a' 


oo oo 


en U J 


0 6 


oK,[R,la, r)-t 
Xcost2:.di==2 Lies “cost 2: dı,— 
Or r—a 
ö ) 
Damit können wir nun unser Randwertproblem lösen, wenn sich der 


f .. (8). 


„s(/R) 
S or 


b) Quellpunkt Qin beliebiger Lage gegenüber dem ihn 
befindet. Der Radius des Zylinders sei gleich der Längeneinheit, er sei mit einem Medium vom 
spezifischen Widerstand o; erfüllt, in einem Medium mit dem spezifischen Widerstand o, ein- 
gebettet und seine Achse falle mit der z-Achse eines Zylinderkoordinatensystemes zusammen. 
Bedeutet R (a, r,9,z) den Abstand des Feldpunktes P (r, 9, z) vom Quellpunkt @ (a, o, o), 
wobei zunächst « #1 sein soll, dann kann man für das Potential im Äußeren und im Inneren 
des Zylinders den Ansatz machen: 


| SS 
UN (r, ©, 2) al | = [| SS K„(r ti) x@(t; a) - cos ..). 21 er eig 
| : n=—© 
(i) J "o®; | 1 ı 2 () Ä | | ” 
U (r, 9,2) = im | nr | 2,%0 Di WPto)sconng|-esterdilr<t| 
J 


Dabei ist J die Stärke des aus dem Quellpunkt in das umgebende Medium fließenden Stromes 
und? =1 oder i—=2, je nachdem aZ 1. Die in den Zusatzpotentialen auftretenden Funk- 
tionen %, und %,„, die so beschaffen sein sollen, daß unter dem &-Zeichen nach r differenziert 
werden darf, sind folgendermaßen aus den Randbedingungen zu bestimmen. Zunächst folgt 
aus (1): 


I,(t)-yP (t; a) — Ku(t)  KOlt; a) =0 G=e1,2.t nel) EI E+2 2ER 
und aus (2) ergibt sich wegen (8) und (8°): 
Isle Klaren) 
Lt) vet; — a Kr ()- xl; a) =(, — 8 3 
) vn (ta) 2 )xp(l;a) =(e 01) Kerken (a<1) 
(=, 31,5 2,50 Wal109. 
Setzt man noch 
An (t; 01 02) = 01° Kult) Int) — Ent) Kr) >0, r=0, +1, 42% °°'), 
dann liefern (10) und (10) die gesuchten Funktionen: 
s I.) 1I.()- K,(l@ de >T1T;,: =1; 
xt; a) = (0, — 91)' Z meh j © ( 3 
KA) - In(at)- a<1;i=2; 
Ky(t) (IC) HERnlarHt) a > 1er: 
Ant) KA) Inlart)- -a<1;i=2; 
Für a =1fällt Er obere Index : weg und man hat zunächst im Ansatz (9) nach (3’) das 


2 #01 ' 0 
01 02 


wegen (8°) die rechten Seiten von (10’) durch nz . U)  K„(t)] zu ersetzen, so daß an 
die Stelle der Gleichungen (11)und (11’) jetzt die folgenden treten: 


»—=0, EL, E22. EU 


vet; a) = (0 — 91) n=0, +1, +2 .- . (11). 


o; zu ersetzen durch . Während aus (1) wieder die Gleichungen (10) folgen, sind nun 


=) +1] 


* ni == Oy ib, 1% 2 ah 2 LE 
md EN) HL; vl 


SS 
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Wegen /,(0) =0 für n > 0 verschwinden für a = 0, also wenn der Quellpunkt auf der 
Zylinderachse liegt, sämtliche x, und y, (n> 0) identisch, die Summen in (9) reduzieren sich 
bloß auf ein einziges, von p unabhängiges Glied, das Feld ist also zylindersymmetrisch >). 

‘Mit Rücksicht auf die Anwendung dieser Ergebnisse auf geoelektrische Messungen in 
einem Stollen möge der 


4. Sonderfall 


02=% etwas ausführlicher erörtert werden, wobei wir annehmen wollen, daß zunächst 
a)Feldpunkt und Quellpunkt außerhalb des Zylinders liegen 
sollen, also a >1lundr >1 seien. Aus (9) und (11) ergibt sich sodann für das Potential: 


Jo 1 2 BIETE REN ER (20) ° I,tt) e | 
u) for re a I { N } ZEnNE on a nos Na 
&% (r, ®, 2) nn I? : | | > Kl) cos N 2 costz “| - 


und wir wollen zeigen, daß as ddurch das Integral dargestellte Zusatzpotential auch dann end- 
lich bleibt, wenn sich der Feldpunkt dem Quellpunkt nähert, wie es ja sein muß. Für t- 0 und 
n>]1 erhält man zunächst aus den entsprechenden Entwicklungen der Funktionen I,„und K,: 
Kr" tt). K,(@a®t) °I,(t) 1 

K,(t) 2nlar)"’ 


Nn=—o 


die Reihe konvergiert also rascher als eine geometrische mit dem Quotienten <ır Kür 
t> O0 aber hat man, wenn etwar>a: SR 


| K.tr t)° Kn(a t)° In (t) | a (a "i) 12-1) 

Kult) K,-ı(t) 
und daher mit Benützung der asymptotischen Entwicklungen von Debye®*) mit den Ab- 
kürzungen 


u (t) = ye—tn: und Ft, n) = w,(t) — n-Yr Tg te 
N 
für genügend großes n: 
Ent: t) Intl) ee . e- 2. [Fü -t,n) — Flt,n =1)7 j 
K,-ı(t) 2 w,(t) 


Setzt mana-t=t+(a— 1) t=t-+7r,so wird 

n 
wobei F, und riR Mittelwerte der Ableitungen von F nach t bzw. n bedeuten. Eine einfache 
Rechnung liefert sodann die für t Ss n gültige obere Schranke: 


Fi+v,n— Fi,n—-1)=F,- +, 1m ml): (a— 1)— ArTg 


ll) ISBN. we) +n un une 
Ka-ı(t) 2 un(at) w(t) —n t 
Nun ist aber 
Sn Le i 
u (er — 1)? 


womit wegen & =2 (a — 1)>0 nicht nur die Konvergenz der Reihe sondern auch die Existenz 
des Integrals selbst für 9 =z =0 gesichert ist. 

Wir wollen nun noch den praktisch besonders wichtigen Fall, daß 

b) Feldpunkt und Quellpunkt auf dem Zylinder liegen, auch hin- 
sichtlich der numerischen Durchführung eingehender behandeln. Aus (9) ergibt sich zunächst 
wegen (3’) und (11’’) für das Potential: 


U (9, 2) = are Er | | > Fr N. - K„(t)) cos 2 cost? 2 (12). 
/ on fen 


5) Vel.z.B. Fr. Ollendorff, Erdströme. Berlin 1928, S. 55 ff. 
6) Vgl.z.B. Jahnke-Emde.Funktionentafeln. 3. Aufl., Leipzig u. Berlin 1938, S. 140. 
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Man kann nun auf eine der Ableitung von (8) analoge Weise zeigen, daß (6) auch für a =r 
gilt, so daß man 1/R mit dem Integral vereinigen kann, wodurch man wegen 


Ink K,—Kı h=- 
für das Potential den eleganten Ausdruck bekommt: 
J: 01 f 5 K,(t) B 7 
T ER i Ei ET, -costz' VE REN Br.- 
U(9,2) = Im P3 1. Kuh) cos np|-costz-di (12). 
ö 
Obwohl hier wegen | En Be (?+ n?)1? für 0O<p <a die Reihe gleichmäßig konvergiert, 


RE 
eignet sich (12°) doch nicht so gut für die numerische Auswertung wie (12), wie wir gleich sehen 
werden. Man kann nämlich in (12) für den numerisch ausschlaggebenden Teil wenigstens die 
Summation ausführen, indem man zunächst die folgende Umformung vornimmt: 


K K, 
ar WR K, —21, Aa un (t). 
ae RI 
Neben einer naheliegenden Rekursionsformel für Ki erweist sich zur Berechnung der 


Produkte L, (t) =I,„(t) K,(t) bei dem beschränkten Umfang der üblichen Tafeln die folgende 
als nützlich: 


k 3 
TER = Ln-2 Sr = Lin z 


k 1 
x f —2: 1’ K-4 |, —2L.+3L+ + 1m.” al): 
Mit den Abkürzungen 
wann ni Bee it — 300 1? n? +-A nt 
TE n 576 u 
findet man ferner für große n mittels der Deb yeschen asymptotischen Entwicklungen: 
3 1 1 2. (®?— An?) 
9 ä AL 9 BEN 
DER; Fa a 21er >) De 8%7 
ne 1 Dr 2-31 —4n?) . i au% ne 
LsCKAT: ao 8 w, Tg er En 
und daher: f 
6 a LEE (1% 9.10‘ “) | at 13 
gta IT — pn) TE De 


Im Gegensatz zu (12’) beginnt also hier die Entwicklung mit dem Glied von der Ordnung (— 2), 
and für dieses kann man auch die Summation ausführen. Aus der bekannten Fourier- 
Reihe für &of it (nr — @)’): 
ED ae Ri. a  &ftr — 9) 
Fe 21002 Sintn 
ergibt sich nämlich durch Differentiation nach it und einige einfache Umformungen: 
| Pil,o)=2: BACH cn I (a) np = a 
ne R (nt —- 9) Sin t(n — Cojtp on 
4 t Sintn Eintn Sin®tn 


Da diese Entwicklungen schon für n > 4 für den hier verfolgten Zweck hinlänglich genau sind, 
so können wir für die Summe in (12) schreiben: 


B> G„(t)-cosnp — Gyl(t) +2 3|@ (t) — la) | sr» | 


n=— © 


. (14) 


Ph arg 


(15) 
= Pu, 9) +2 [0 (t) + m(t)] cos np | 


”) Vgl.z. BBK.Knopp, Unendliche Reihen. 3. Aufl., Berlin 1938, S. 391. 


ver 
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und damit schon die Bedingungen feststellen, unter denen die 
c) Konvergenz des Integrales in (12) gesichert ist. Aus den mit 


ö=1g2 — y =0,6931 .. .— 0,5772... =0,116... 

bestehenden Entwicklungen bei t =: 

Gl) = gt Zr) eig + gl +26--202)- +--- 

Gt) — let o+4)+; AU 2 te... n>1 

4 Weges Aa 4 ’ "N 2n2-(m— 1) (n>]) 

2 
Pl, 9) =— 5, 17 m 3022 (ag)? +15 (a pi] 4 I 
folgt zunächst die Existenz des Integrales an der unteren Grenze. Da ferner für 0 <n<t 
1 1 
G,(t) Saar gen > 
und für große { 
P(t,g) m le Hp) er + 2m ertanmn un 2. /00.22...(16) 


so existiert das über die ersten Glieder 
der Summe in (12) erstreckte Integral an 
der oberen Grenze ebenfalls, u. zw. auch für 


o=2 U, 


Wie ferner aus (16) ersichtlich ist, existiert 
auch das über P(t,@) erstreckte Integral an 
der oberen Grenze für 


+20. 


Da aber P(t,0) für {> oo nur von der 
ersten Ordnung verschwindet, so existiert das 
Integral nicht mehr, wenn zugleich 9— 0 
und z— 0, also wenn der Feldpunkt in den 
auf dem Zylinder liegenden Quellpunkt 
hineinrückt, was ja auch hier aus den be- ) 
reits am Schluß von 2d erwähnten Grün- 0 30 60 90 he 70 180 
den klar ist. 


Die beiden Bilder 3a und 3b, die den 
Verlauf des Potentiales ent- 
lang des durch den Quell- 
punkt gehenden Kreises und 
der durch ihn gehenden Er- 
zeugenden des Zylinders dar- 
stellen, zeigen unmittelbar, 
daß der Betrag des Zusatz- 
potentials % bei Annäherung 
an den Quellpunkt erheb- 
lich langsamer über alle 
Grenzen wächst als das redu- 
zierte Potential . 

Für die numerische 
Auswertung des Integrale 
in (12) ist noch eine 

d) Abschätzung 
dervong(t)und r„(i) 
in (15) gelieferten Bei- 


träge erforderlich. Wir dür- 
fen uns dabei offenbar auf Bild 3. Verlauf des von einer Punktquelle auf der Wand eines zylindrischen Hohlraumes 


PR . erzeugten Potentiales Yund des Zusatzpotentiales U entlang des durch sie gehenden Paral- 
z=0 beschränken, ersetzen im ec Erenzenden 


Bild 3a 


Bild 3b 
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f B Mt) cosn e\ - dt zunächst die 5 durch hr dn und können dann das über den Bereich 
o In= n=5 
0 Ss ß 4 < nerstreckte Doppelintegral durch Erntarune von Polarkoordinaten n =r cos Ö, 
t=r:sin d in der folgenden Weise approximieren: 

Lee) ya "COSNY 

e| | ar apa di 

0 4 
| far R Dr u| 
dr : a 
Be 7 | 005 (w°r 005 0): sine 9: an Deos NY | vera 
4 ö yio—m 


Der Fehler, den wir hier dadurch begehen, daß wir im zweiten Glied den Viertelkreis durch ein 
Stufenpolygon ersetzen, kann sicherlich vernachlässigt werden. Wegen der Poissonschen 
Integraldarstellung ®) hat man aber 


n/2 
ER 
cos(p*r cos d)- sind dd = 57 BAER 
0 
und daher: 
ne 3% Je(r - 2). 1 : f tt di 
PU, er]; den Sp [ r& 7 A 8 ND° j (1? + n?y’R’ 
u 4 Vıs—n: 


wo nun das erste Integral numerisch und das zweite sogar elementar ausgewertet werden kann. 
Damit findet man die folgenden Näherungswerte für die von q, (t) herrührenden Beiträge: 


N BE ee 2 I nn 
2 BEREITEN: a 
».f[2 Mm(t) cosng|- dt: 7 —7 1 —5 4 ed ed 
0 In=5 


Danach dürfen die aus dem Glied (—4)-ter Ordnung r, (t) entstehenden Beiträge umso eher 
vernachlässigt werden, als r, (t) für £ > 0 zweimal das Vorzeichen wechselt. 


Die im Vorangehenden gemachten Ansätze sind einer naheliegenden 


5. Verallgemeinerung 


fähig, die eine Lösung der entsprechenden Randwertaufgabe für geschichtete Kugeln und 
«Zylinder ermöglicht. Wir wollen uns auf den Fall beschränken, daß der Quellpunkt außerhalb 
der Kugel, bzw. des Zylinders liegt, und nehmen an, es sei: 


e=AM fü r>ny pp = fin n>r Dr; Oo 9 EUNR RAU I 
Für eine derartige 


a) geschichtete Kugel führen mit dem Ansatz: 


n<riu(nd)= + Dan: 7771: P,(0) 
n—=(0 
n<r<n:ald)=g 2 (Ban + Baer): Pu(6) 


T<1a:Us(r, d) = + Diner: P„(9) 
n=0 


die Randbedingungen (1)und (2)für jedes nzu einem inhomogenen System von je vier Glei- 
chungen für die Koeffizienten &,, Pa, Pr, Y„ mit nicht verschwindender Determinante. 


8) W.Magnusu.Fr.Oberhettinger,l.c. S. 26. 
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In analoger Weise ist auch ein 
b) geschichteter Zylinder zu behandeln, wobei wir uns wieder auf den Fall 
a >r, beschränken wollen. Für das reduzierte Potential liefern mit dem Ansatz: 


oo 


1 oo 
n<r:ul®, )=4+- B03 K„(r *t)*yP (t) - cos nal costzc 


N=—w 


oo 


NH<Tr<rN:ulr,®, )=5 +—: | > [Ku(r-t) xt) 4+1,(r ty (t)]- cos er) -costz di 
0 


Bee 


1 
ER A LAK) er Ya 


LES EETTERRONFM a 
RB | 2 ht t) wu (t) cos 5 costz:dt 


ö [8 


die Randbedingungen (1) und (2) 
für jedes nein inhomogenes System 
von je vier Gleichungen für die 
Funktionen yP, y2, y2, „EP mit 
nicht verschwindender Determi- 
nante. Damit kann auch das 
Problem des Stromübertrittes zwi- 
schen der Erde und einem Metall- 
rohr?) vollkommen streng erledigt 
werden, ohne daß dabei die radiale 
Strömungskomponente und der 
innere Radius des Rohres außer- 
acht zu bleiben brauchen. 

Man sieht nun auch sofort, 
wie das in Frage stehende Rand- Z 
wertproblem für eine Kugel und 
einen Zylinder mit mehreren 
konzentrischen, bzw. koaxialen 
Schichten zu lösen ist. Für die 
„Belegungsdichten“ x, Yu --- 
erhält man dabei rasch kompli- 
zierter werdende Ausdrücke, doch 
genügt für die Geoelektrik meist 
schon die Kenntnis einer einzigen 
von ihnen zur Konstruktion des 
Potentiales in jenem Medium, auf 
dessen Begrenzung die Elektroden 
und die Sonden gesetzt werden, 
weil für das angrenzende, nämlich 
die Luft, U =const ist. 

Diese Ergebnisse ermöglichen 
nun zwar — natürlich nach ent- 


sprechender Idealisierung — die 

Behandlung zahlreicher in der ° : e = Y 
praktischen Geoelektrik auftre- 30° 60° 90 120 250 23 
tender Fälle, wie an anderer Stelle en 

gezeigt werden soll, aber noch viel Bild 4. Verlauf des von einer Punktquelle auf der Wand eines 


P B . . ‚ zylindrischen Hohlraumes erzeugten Potentiales 
interessanter sind für sie solche, bei ; s a 


denen die Trennungsflächen verschiedener Medien elliptische und parabolische Zylinder oder 
Ellipsoide sowie deren Ausartungen sind. Umaber mit den hier benützten analogen Ansätzen 
auch für solche Fälle die Randwertaufgabe lösen zu können, müßte man insbesonders der Gl. (6) 
entsprechende Darstellungen von 1/Rmittels Mathieuscher, Weberscher, bzw. Lame&- 
scher Funktionen verwenden, doch scheinen solche bisher noch nicht bekannt zu sein. 


®) P.Franku.R.v.Mises, Die Differential- und Integralgleichung der Mechanik und Physik, II, 
2. Aufl., Braunschweig 1935, S. 687 ff. 


Eingegangen am 15. September 1952. 
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BUCHBESPRECHUNGEN 


L. Bieberbach (ehem. o. Prof. d. Univers. Berlin), 
Theorie der geometrischen Kon- 
struktionen (Lehrbücher und Monographien 
aus dem Gebiete der exakten Wissenschaften. Mathem. 
Reihe Band 13.) 102S. m. 102 Abb. Basel 1952. 
Verlag Birkhäuser. Preis brosch. 15,60 Fr., geb. 
18,70 Fr. 


Verf. gliedert sein Buch nach den verschiedenen In- 
strumenten und Hilfsmitteln, die bei den Konstruk- 
tionen zu verwenden sind, und nach deren Anwendungs- 
vorschrift. Auf Konstruktionen mit dem Lineal 
allein in unbegrenzter und begrenzter Ebene folgen 
Konstruktionen mit dem Lineal und einem, zwei, drei 
festen Kreisen, mit Zirkel und Lineal, mit Zirkel allein, 
mit Lineal und Eichmaß. Es werden Instrumente wie 
Rechtwinkellineal, Winkelhaken, Einschieblineal, Zim- 
mermannshaken, Steckzirkel, Ellipsen- und Kissoiden- 
zirkel besprochen. Der innere Gang ist der von line- 
aren Konstruktionen über quadratische zu Konstruk- 
tionen 3. und 4. Grades im algebraischen Sinne. Verf. 
verwendet viel Sorgfalt darauf, die erlaubten Ver- 
wendungsvorschriften für die Instrumente stets genau 
zu formulieren. Er gibt dann jeweils mit geometrischen 
Mitteln die Lösung der Grundaufgaben und mit 
algebraischen Mitteln einen Überblick über den Be- 
reich der mit dem betreffenden Instrument behandel- 
baren Aufgaben und dessen Grenzen. Die Unmöglich- 
keitsbeweise für die Dreiteilung des Winkels und die 
Quadratur des Kreises werden ausführlich in aller 
Strenge durchgeführt. 

Verf. vereint in seinem Buch die Beherrschung der 
weitest reichenden modernen algebraischen und zah- 
lentheoretischen Hilfsmittel mit sorgfältiger Kon- 
struktion des einzelnen Falls bis zum letzten Zirkel- 
schlag und mit ansprechender Eleganz der Darstellung. 
Leider wird die Freude an seinem Buch getrübt durch 
einen polemischen Angriff, den er bei Gelegenheit der 
Dreiteilung des Winkels auf 8.53 —56 gegen van der 
Waerdenführt: VanderWaerdenmeintin 
seinem Buch über Moderne Algebra Bd. IS. 194, daß 
die (historische, von den Griechen gestellte) Frage ‚die 
nach einer einheitlichen Konstruktion sei, die für alle 
Winkel gelten soll. Auf die Betrachtung von Drei- 
teilungskonstruktionen spezieller Winkel geht er in 
diesem Buch nicht ein. Das besagt aber nicht, daß sie 
für ihn kein Interesse hätten. Daß er ‚autoritäre An- 
gaben oder Behauptungen‘ mache, die ‚uns Geome- 
ter belehren‘‘, kann der Ref. mit dem besten Willerrin 
dem erwähnten Buch nicht entdecken. ‚‚Der bewährte 
Standpunkt, daß jede Fragestellung und Meinung 
“erlaubt sei‘, wird von vander Waerden nir- 
gendwo angetastet. 


Halle/Saale. Ott-HeinrichKkeller. 


Dr.-Ing. W. Matz, Berechnung der Aus- 
mauerungstählerner Gefäße. (Verfah- 
renstechnik in Einzeldarstellungen Heft 1). VIII + 
758. m. 15 Abb. Berlin/Göttingen/Heidelberg 1953. 
Springer-Verlag. Preis brosch. 10,50 DM. 


Die Ausmauerung stählerner Behälter der chemi* 
schen Großindustrie ist schon seit Jahrzehnten in Ge- 
brauch. Trotzdem fehlte ein theoretisch begründetes 
Berechnungsverfabren. Die Schäden an den Kochern 
in der Zellstoffindustrie (Ablösung der Ausmauerung 
einerseits oder Explosion ganzer Behälter andererseits) 
führten allmählich zu der Erkenntnis, daß es sich um 
einen Verbundkörper handelt, bei dem der Wärme- 
durchgang und die verschiedene Wärmedehnung des 
Stahles und des Mauerwerks von grundlegender Be- 
deutung sind. Versuche und theoretische Betrachtun- 
gen ließen erkennen, daß der Verbundkörper (Stahl- 
mantel und Ausmauerung) nur dann von Bestand sein 
kann, wenn man ihn durch Verwendung von quell- 
fähigem Kitt eine Vorspannung gibt. Mat z hat nun 


die beiden. Grundprobleme des Wärmedurchgangs und 
der Wärmedehnung für den Verbundkörper Stahl und 
Mauerwerk in Verbindung mit einem quellfähigem 
Kitt tbeoretisch klar aufgezeigt und sie für die Praxis 
bis zur gebrauchsfertigen Formel und Rechenvorschrift 
entwickelt. 

Die Abstimmung der verschiedenen Faktoren 
(Stärke des Behälters, Stärke der Ausmauerung, Wahl 
des Kittes), setzt trotz der übersichtlichen und im 
Grunde genommen einfachen Theorie sehr große; Er- 
fahrungen voraus, weil Fragen der Wärmedehnung 
und der Festigkeitslehre aufeinander abzustimmen 
sind. Zahlenbeispieleam Schluß jedes einzelnen Kapi- 
tels tragen sehr viel zum Verständnis der theoretischen 
Ausführungen bei und zeigen die praktische Bedeu- 
tung der Theorie für die Bemessung der Einzelteile. 

In dem Buch ist nur der stationäre Wärmefluß 
erfaßt, eine Arbeit über die Vorgänge beim Anheizen 
und Abkühlen ist aber vom Verfasser in Aussicht 
gestellt. 

Das Buch ist ein sehr wertvoller Beitrag für die Wei- 
terentwicklung des Behälterbaues und wird von allen 
an diesen Aufgaben Beteiligten begrüßt werden. 


Dresden. Henn. 


W.Matz, Aufgabensammlung zur 
Thermodynamik des Wärme- und 
Stoffaustauschesin der Verfahrens- 
technik. XII+138S. m. 100 gelösten Auf- 
gaben und 29 Abb. Darmstadt 1953. Verlag Dietrich 
Steinkopff. Preis brosch. 16,— DM, geb. 18,— DM. 

Das vorliegende Werk ist eine Ergänzung zu dem 
gleichnamigen Buch des Verfassers aus dem Jahre 
1949. Es ist ohne dasselbe nicht für sich allein zu 
verwenden, da auf die Entwicklungen, Gleichungen 
und Abbildungen der zugeordneten Abschnitte 
laufend Bezug genommen wird. Nach dem Vorwort 
war es die Absicht des Verfassers, die Anwendung 
der in der Thermodynamik vorgetragenen Theorie an 
solchen praktischen Beispielen zu erläutern, wie sie 
„in der Industrie häufig an den Betriebsingenieur 
und Betriebschemiker herantreten‘“. 

Der praktische Betriebsmann wird jedoch die 
Aufgabensammlung enttäuscht aus der Hand legen, 
denn er findet unter den hundert Aufgaben allenfalls 
zwei Dutzend, die ein konkretes praktisches Problem 
zum Gegenstand haben und zahlenmäßig durch- 
rechnen. Alle anderen behandelten Aufgaben sind rein 
theoretischer und spekulativer Natur; sie bilden als 
solche allerdings eine vorzügliche und dankenswerte 
Ergänzung und Vertiefung des Hauptbuches und 
werden auch dem Fachlehrer als Prüfungsfragen 
willkommen sein. 

Auf zwei unterlaufene Fehler sei hingewiesen: In 
Aufgabe 4 ist die Berechnung von AW aus den 
spezifischen Dampfgewichten nicht exakt, da auch 
die bezogene trockene Luftmenge ihren Zustand 
ändert; in Aufgabe 15 ist die Integration von d@ zu 
beanstanden. 


Dresden. HrBa tn 


Prof. Dr. Gerhard Kowalewskiy. Einführung 
in die Analytische Geometrie. Vierte 
verb. Aufl. 364 S. m. 112 Abb. Berlin 1953. Walter 
de Gruyter & Co. Preis geb. 26.— DM. 


Das Buch ist aus Vorlesungen des Verf. an den 
Universitäten Leipzig, Greifswald und Bonn ent- 
standen und stellt eine erste Einführung in die 
Analytische Geometrie der Ebene und des Raumes 
vor. Die vierte Auflage ist einer gründlichen Revision 
unterzogen und teilweise ergänzt worden. Verf. be- 
handelt die Analytische Geometrie mit Verwendung 
der Vektoren, die hier „Strecken‘‘ genannt werden. 
Die Grundlagen der Streckenrechnung werden zu- 
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nächst abgeleitet und hierauf die schief- und recht- 
winkligen kartesischen Koordinaten, die Polar- 
koordinaten in der Ebene und im Raum und die 
Zylinderkoordinaten eingeführt. Homogene Koordi- 
naten dienen zur Erklärung der uneigentlichen 
Elemente des Raumes. Mit Hilfe projektiver Koordi- 
naten trägt Verf. die projektive Geometrie auf der 
Geraden und in der Ebene und die Lehre von den 
Kollineationen und Korrelationen, insbesondere der 
Polaritäten vor. Es folgen die allgemeine Theorie der 
ebenen algebraischen Kurven und Beispiele von be- 
sonderen höheren algebraischen und transzendenten 
Kurven. Die Theorie der quadratischen Formen 
dient zur Behandlung der Kurven 2. Ordnung und 
2. Klasse. Es werden auch Büschel von Kurven 
2. Ordnung und Scharen von Kurven 2. Klasse be- 
trachtet. Im Zusammenhang mit letzteren stehen die 
Fokaleigenschaften der Kegelschnitte, die zur Ein- 
führung elliptischer Koordinaten in der Ebene ver- 
anlassen. Die Kreisbüschel wieder führen zur Inver- 
sion und diese zur stereographischen Projektion der 
Kugel. Die projektive Geometrie des Raumes wird 
mit Hilfe projektiver Punkt- und Ebenenkoordinaten 
undder Plücker schen Linienkoordinaten behandelt. 
Sie wird einschließlich der Elemente der Linien- 
geometrie bis zu den Flächen 2. Ordnung und 2. Klasse 
vorgetragen. — Kowalewskis Lehrbuch ist aus- 
gezeichnet geeignet, die Studierenden der Anfangs- 
semester mit den wesentlichen Begriffen und Ver- 
fahren der Analytischen Geometrie vertraut zu 
machen. Sein bekannt leichter und flüssiger Stil und 
die Klarheit und anschauliche Kraft seiner Sprache 
machen die Lektüre dieses Buches zu einem Genuß. 
Er versteht es meisterhaft, den Leser zur selb- 
ständigen Weiterverfolgung der im Lehrbuch ein- 
geleiteten Gedankengänge anzuregen. Auch die neue 
Auflage des Kowalewskischen Werkes wird dem- 
nach eine gern und mit vielem Nutzen gelesene Ein- 
führung in die Analytische Geometrie bleiben. 


Freiberg (Sa.). W. Schmid. 


H. Rutishauser, Automatische Rechenplan* 
fertigung bei programmgesteuerten Re- 
chenmaschinen. (Mitteilungen aus dem Insti- 
tut für angewandte Mathematik an der ETH. Zürich 
Nr.3.) 458. m.8Abb. und 3 Strukturdiagrammen. 
Basel 1952. Verlag Birkhäuser. Preis brosch. 5,70 
Schw. Fr. 

Neben der rein technischen Weiterentwicklung der 
programmgesteuerten Rechenautomaten tritt in letzter 
Zeit die Notwendigkeit nach Vereinfachung der oft 
recht zeitraubenden und fehlerempfindlichen Auf- 
stellung der Befehlsreihe (Rechenprogramm) immer 
mehr in den Vordergrund. Verschiedene Wege sind 
beschritten worden: Formelgemäße Operationsaus- 
führung im Rechenautomaten, Verwendung von 
Programmbibliotheken, Einsatz spezieller Hilfsgeräte 
zur Programmfertigung. Verf zeigt, wie der Rechen- 
automat selbst als Planfertigungsgerät benutzt werden 
kann. Das Ziel ist die Verkürzung der Vorbereitungs- 
zeit und insbesondere die Elimination manueller 
Arbeit und damit die Beseitigung von Fehlerquellen. 

Die vorgeschlagene Methode besteht darin, daß 
alle auszuführenden Formeln arithmetisiert werden, 
indem jedem ihrer Bausteine (Klammern, Operations- 
zeichen) charakterisierende Zahlen, den Operanden 
die vorgeschriebenen Speichernummern zugeordnet 
werden. Mit dieser Zahlenfolge wird die Aufgabe voll- 


ständig beschrieben und damit die automatische 
Auflösung in eine Befehlsreihe gemäß der für das 
spezielle Gerät geltenden ‚Spielregeln‘‘ möglich. 
Dabei können auch die Operanden mit laufenden 
Indices berücksichtigt werden — das Ergebnis sind 
dann entweder teilweise zyklische Rechenpläne 
(Minimum an Befehlen) oder die vom Verf. besonders 
für langsame Maschinen empfohlenen teilweise ge- 
streckten Rechenpläne (Elimination von Zähl- und 
Sprungoperationen). 

Im einzelnen werden vom Verf. für weite Problem- 
klassen (Beispiele: Matrixinversion, Eliminations- 
verfahren für lineare Gleichungssysteme nach Gauß- 
Banachiewicz) die zur automatischen Herstellung 
zyklischer und gestreckter Rechenpläne erforderlichen 
Überlegungen durchgeführt und die Struktur der be- 
nötigten Superpläne untersucht. 


Dresden. N. J. Lehmann. 


Noordhoff’s Wiskundige Tafels in 5 decimalen, 
Tafel H (Fünfstellige mathematische Tafeln). Heraus- 
gegeben vonP. Wijdenes, Amsterdam. 5. Aufl. 
VIII + 2698. Groningen (Holland) 1953. Verlag 
P. Noordhoff. Preis geb. 8.75 f. 

Eine kurze Angabe des Inhalts des in 5. Auflage 
vorliegenden Tafelwerks verdeutlicht wohl am besten 
seinen Wert. 

Tafel I: Briggsche Logarithmen, fünfstellig von 
1 bis 10 000, sechsstellig von 10 000 bis 
12 000; 

Fünfstellige Legarithmen der trigono- 
metrischen Funktionen: 0° (1’’) 20’ (10°) 
2°°(17) 45°; 

Umwandlungen zwischen Altgrad, Neu- 
grad und Bogenmaß; 

Fünfstellige Werte der trigonometrischen 
Funktionen von Minute zu Minute, dgl. 
von ctg&a für 0° <a< 3° von Sekunde 
zu Sekunde, dgl. von ctg a für 3’ <a<<10° 
von zehn zu zehn Sekunden; 

Natürliche Logarithmen, fünfstellig von 
1 bis 1000 sowie für die Primzahlen von 
1000 bis 10 000, achtstellig für die Prim- 
zahlen von 1 bis 1000; 

e? und e-%2, sechsstellig: O0 (0,001) 

1 (0,01) 4; 
sinh z’und cosh x, sechsstellig: 0 (0,01) 

44 (051)27; 

Primfaktorentabelle von 1 bis 12 000; 


n?, n®, Yn, Yn, E siebenstellig von 1 bis 
000; IL 


Tafel II: 


Tafel III: 
Tafel IV: 


Tafel V: 


kleinere Tabellen für die Fakultät, Ex- 
ponentialintegral, Integralsinus und -co- 
sinus, Fehlerintegral, Besselsche Funk- 
tionen Ju(®), Jı(®), Yo(®), Yy(e). 
Zusammenfassend ist zu sagen, daß die Tafel dem 
im Vorwort genannten Zweck — Studenten an Uni- 
versitäten und technischen Hochschulen eine wich- 
tige Hilfe zu sein — voll gerecht wird. Darüber 
hinaus wird sie auch für den in der Praxis stehenden 
Ingenieur sehr wertvoll sein. Hervorzuheben ist 
noch die außerordentlich übersichtliche Anordnung, 
die durch verschiedenfarbiges Papier für die einzel- 
nen Abschnitte noch beträchtlich unterstützt wird, 
sowie der klare Druck und die gute Gesamtaus- 
stattung. 


Dresden. A. Schubert. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen, ausführliche Besprechung bleibt vorbehalten. 


Dr. A.Müller, Nomographie für die 
technische Praxis 2688. m. 97 Abb. 
Leipzig 1952. Fachbuchverlag GmbH. Preis geb. 
9,80 DM. 


J. P6res, (Prof. & la Facult6 des Sciences de Paris), 
M&canique Gö&ne&rale VI-+407S. m. 


64 Abb. Paris 1953. Masson & Cie. Preis brosch. 


2110 frs., geb. 2545 frs. 
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Dr. Günther Pickert, Analytische Geo- 
metrie. (Mathematik und ihre Anwendungen in 
Physik und Technik Reihe A Bd. 24) X + 397 8. 
m. 76 Abb. Leipzig 1953. Akademische Verlags- 
gesellschaft Geest & Portig KG. Preis geb. 26,— DM. 


Dr. Ernst Schmidt (o. Prof. a. d. Technischen Hoch- 
schule München), Einführungin dieTech- 
nische Thermodynamik undiin die 
Grundlagen derchemischen Thermo- 
dynamik. 5. berichtigte Aufl. XVI +520S. m. 
244 Abb. u. 69 Tabellen. Berlin-Göttingen-Heidel- 
berg 1953. Springer-Verlag. Preis geb. 30,— DM. 


Dr. Ralph W.Detra, The Secondary Flow 
in Curved Pipes. (Mitteilungen aus dem 
Institut für Aeredynamik an der ETH. Zürich. 
Nr. 20) 508. m. 22Abb. Zürich 1953. Verlag 
Leemann. Preis 15,60 Schw. frs. 


ZeitschriftfürFlugwissenschaften(ZFW), 
Organ der WGL. Heft 1. Braunschweig 1953. 
Preis vierteljährlich (3 Hefte) 12,— DM. 


Prof. Dr. F. Böhm (Prof. a. d. Universität München) 
Versicherungsmathematik II. (Samm- 
lung Göschen Band 917/917a). Zweite, vermehrte und 
verbesserte Aufl. 205 S. Berlin 1953. W. de Gruyter 


& Co. Preis 4.80 DM. 
Jahrbuch 1952 der Wissenschaftlichen 
Gesellschaft für Luftfahrt e. V. (WGL). 


148 S. Braunschweig 1953. Friedrich Vieweg & Sohn. 
Preis geb. 24,— DM. 


Wissenschaftliche Zeitschrift der Uni- 
versität Leipzig, Mathematisch -naturwissen- 
schaftliche Reihe Nr.3. 1952/53. S 169—308. Selbst- 
verlag der Universität. 


H. v. Jürgensonn, 
Festigkeit im 


Elastizität und 
Rohrleitungsbau. 
Statische Bereehnung der Rohr- 
leitungen und ihrer Einzelteile. 
2. neu bearb. Aufl. VIII + 379 S. m. 236 Abb. u. 
38 Zahlentafeln. Berlin-Göttingen-Heidelberg 1953. 
Springer-Verlag. Preis geb. 39,— DM. 


NACHRICHTEN 


Ansprache des Vorsitzenden der Ges. angew. Math. 

Mech. A. Walther bei der Trauerfeier für Lud- 

wig Prandtl in der Universitätskirche Göttingen 
am 20. August 1953. 


In tiefer Trauer steht die Gesellschaft für angewandte 
Mathematik und Mechanik an der Bahre ihres Ehren- 
vorsitzenden und bringt ihm einen letzten Gruß dar. 

Als unsere Gesellschaft 1923 begründet wurde, konnte 
keinerlei Zweifel darüber bestehen, daß Ludwig 
Prandtlder Posten des Vorsitzenden gebührte. Da- 
zu hatten ihn nicht nur seine allgemein anerkannten 
wissenschaftlichen Leistungen qualifiziert, sondern auch 
sein hohes Ansehen als ein reiner und gütiger Mensch. 
Hinzutrat die Tatsache, daß in seiner Person das Stre- 
ben unserer Gesellschaft nach Einheit von Theorie und 
Praxis, von Wissenschaft und Leben in einzigartiger 
Weise verkörpert war. 

Er hat den Vorsitz über 20 Jahre lang bis zu dem für 
ihn besonders harten Zusammenbruch 1945 innegehabt. 
Durch seine Autorität, sein vermittelndes Wesen und 
seine Pflichttreue wußte er die Gesellschaft auch durch 
schwierige Situationen mit Glück hindurchzuführen. Als 
sie im Gefolge der Zoneneinteilung auseinanderfiel, ver- 
lor er den Mut nicht, sondern sammelte unter seinem 
Vorsitz bald wieder eine Teilgesellschaft für die bri- 
tische Zone um sich. Auf der Darmstädter Tagung im 
April 1950 erfolgte der Zusammenschluß mit der ent- 
sprechenden Gesellschaft in der amerikanischen Zone, 
und es kam wieder eine einheitliche deutsche Gesell- 
schaft für die ihm so sehr am Herzen liegenden Fach- 
gebiete zustande. Diese Neubildung war, wie ich selbst 
miterleben durfte, eine ganz große Freude für ihn. Eben- 
so sehr freute er sich, daß ihm einstimmig der Ehrenvor- 
sitz der Gesamt-Gesellschaft übertragen wurde. Er 
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nahm dann an unserer Arbeit rege und aktiv teil, so- 
lange ihm das möglich war. Wir bewunderten sein Inter- 
esse und seine Aufgeschlossenheit für alle wissenschaft- 
lichen Fragen, bis Krankheit und Siechtum ihn von uns 
trennten und nur noch unsere Gedanken und Wünsche 
zu ihm gehen konnten. 

Fragen wir nach den tiefsten Wurzeln der Persönlich- 
keit von Ludwig Prandtl, so scheinen sie mir 
darin zu liegen, daß er sich eine echte Kindhaftigkeit 
im biblischen Sinne bewahrt hatte. Sie befähigte ihn 
als Forscher, die Dinge völlig unvoreingenommen zu 
sehen, sich nicht an Vorhandenes zu binden, sondern 
neu und ursprünglich zu denken und zu arbeiten. Ihm 
war es — um mit Goethe zu sprechen — gegeben, zu 
den Müttern hinabzusteigen. Wie ein Kind vermochte 
er immer wieder nach dem Grunde beobachteter Er- 
scheinungen zu fragen, sich am Spiel mit dem rinnen- 
den Wasser, an der Betrachtung der Vorgänge im Wind- 
kanal zu freuen und daraus tiefe, aber doch einfache 
und anschaulich verständliche Erkenntnisse zu ge- 
winnen. Und als Menschen war LudwigPrandtl 
die Reinheit des Kindes eigen, das noch nicht in Kon- 
flikt mit der Welt gekommen ist, sein Gerechtigkeits- 
sinn, sein Wahrheitsgefühl, seine Aufgeschlossenheit 
und Hilfsbereitschaft. 


Was Ludwig Prandtlals Forscher geschaffen 
hat, wird Jahrhunderte überdauern. Er hat seine Spur 
tief in die Wissenschaft und ihre Anwendungen einge- 
prägt. Was er als Lehrer und Mensch war, werden seine 
Schüler und Freunde in Treue bewahren. Die Gesell- 
schaft für angewandte Mathematik und Mechanik neigt 
sich, indem ich diesen Kranz niederlege, in tiefer Ver- 
ehrung vor ihm. Habe Dank für alles, was du uns ge- 
geben hast LudwigPrandtl. 
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